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Resume : Etude et Moderation d’un refrigerateur solaire 

Notre travail de these a pour but d’etudier, dimensionner et de modeliser un refrigerateur 
solaire, l’objectif de la premiere partie est de realiser un refrigerateur domestique qui fonctionne 
en energie solaire photovoltaique avec une faible consommation et un moindre cout. Pour cela, 
trois prototypes de refrigeration a compression de vapeur ont ete realises, dont deux refrigerateurs 
de 160 litres et une chambre froide positive de 1.56 m 3 (1560 litres). Le premier refrigerateur 
(DC-AC) a ete directement alimente par energie solaire en utilisant un onduleur solaire 
(convertisseur DC-AC) et sans changement dans les equipements du cycle frigorifique 
conventionnel; le deuxieme (DC-DC) a ete adopte pour fonctionner directement en energie 
solaire photovoltaique sans utilisation du convertisseur en employant un compresseur DANFOSS 
BD35F qui marche en courant continu et un evaporateur eutectique AC 15 (stockage du froid par 
chaleur latente) pour reduire la capacite de stockage electrochimique; le dernier prototype (DC- 
AC), c’est une chambre froide positive alimentee aussi directement par energie solaire 
photovoltaique. Ces prototypes ont ete dimensionnes avec le logiciel PVSyst V5.59 pour une 
autonomie de 3 jours sans ensoleillement et testes sous les conditions climatiques sahariennes et 
mediterraneennes. La consommation journaliere du systeme frigorifique DC-AC est de 
0.865kWh/jour et un cout global du systeme estime a 462779.46 DA. La consommation 
journaliere du systeme frigorifique DC-DC sous les memes conditions est de 0.473kWh/jour et un 
cout global du systeme de 280436.35 DA, cela induit une reduction de 40%, c.a.d. une reduction 
de 182343.11 DA par rapport au premier systeme DC-AC. Cette reduction est due principalement 
au prix de F onduleur supprime et a la reduction de la capacite de stockage electrochimique. 
D’apres les tests menes sur la chambre froide con£ue (DC-AC), il s’est avere que cette demiere 
est capable de garder sa temperature interieure entre 2 et 8°C, avec une consommation moyenne 
de 1.2kWh/jour sous les conditions environnementales climatiques sahariennes (Ain-Guezzam : 
Tamanrasset). 

L’objectif de la deuxieme partie, est de dimensionner une machine frigorifique solaire a 
adsorption solide/gaz afin de concevoir un prototype capable de fonctionner sous les conditions 
meteorologique algeriennes. Mais les problemes rencontres au cours de soudage nous ont obliges 
d’arreter ce projet et de s’orienter vers la modelisation d’une machine frigorifique solaire a 
adsorption solide/gaz. C’est dans ce contexte que se situe notre travail dont le but est d’etablir un 
modele neuronale capable de predire le Coefficient de Performance Solaire (COP s ) d’une 
machine frigorifique a adsorption solaire. Cette modelisation s'articule sur l'elaboration d'un 

programme sur MatLab®. Le reseau de neurones statiques optimise (RNSO) a une architecture de 

• • • 
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(8-8-1) : {8 neurones dans la couche d’entree; 8 neurones dans la couche cachee et 1 neurone en 
sortie (COP s )}. Le RNSO reproduit fidelement les donnees experimentales tirees de la litterature 
avec une tres bonne precision: une erreur quadratique moyenne de 0.45 %o et un coefficient de 
correlation R 2 > 0.9985 dans le cas du test. Pour apporter encore plus de souplesse a l’utilisation 
du RNSO, un programme en langage MatLab® avec une interface conviviale a ete developpe. 
Une etude a ete aussi menee sur la sensibilite des parametres d’entree sur la sortie (COP s ), dans le 
but de determiner les parametres les plus influen 5 ants sur le COP solaire. Cette etude nous a 
permis aussi de determiner le COP solaire par la methode du reseau de neurones inverse optimise 
(RNIO) et par consequent estimer le temps de calcul de la simulation du RNIO, qui est de 10 
secondes pour le cas du test. 

Mots cles : Refrigeration solaire, adsorption, photovoltai'que, reseaux de neurones 
artificiels directe et inverse. 

Abstract: Study and modeling of solar refrigerator 

Our study is divided in two parts. The first part aims to design a domestic refrigerator 
which operates with photovoltaic solar energy with low power consumption, minimal cost and 
more environmentally friendly. For this, three prototypes with a vapor compression cycle have 
been constructed, two of them have 160 liters and a positive container with 1.5 m 3 (1500 liters). 
The first prototype (DC-AC) was fed directly by solar energy using a solar inverter without any 
modification in the conventional equipment. The second prototype (DC-DC) was adopted to be 
powered directly by solar energy without inverter. This was done by replacing the conventional 
compressor/condenser and the evaporator respectively by a set of DANFOSS BD35F 
compressor/condenser and an AC 15 eutectic evaporator, these modifications led to reduce the 
electrochemical storage capacity. The third prototype (DC-AC) is a container powered by solar 
energy using solar inverter. These prototypes were designed using the software PVSyst V5.59 to 
work for 03 days without sunshine. They were tested under the meteorological conditions of the 
Sahara and the Mediterranean. Results revealed that, the daily consumption of the first 
refrigeration system is 0.865kWh/day with an estimated cost of approximately 462779.46 DA. 
The daily consumption of the second refrigeration system under the same conditions is 
0.473kWh/day with an estimated cost of 280436.35 DA. This induces a reduction of 40%, i.e. a 
reduction of 182343.11 DA in comparison with the first system. This reduction is mainly due to 
the price of the inverter and also due to the reduction of the electrochemical storage. The tests 
conducted on the third system also revealed that the designed container is able to keep the indoor 
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temperature between 2 and 8°C with an average consumption of 1.2kWh/day and under 
environmental conditions of the Sahara (Ain-Guezzam: Tamanrasset). 

The objective of the second part of the study is to size a solar adsorption refrigerator 
machine to work under the Algerian climatic conditions, but problems met by wolding obliged us 
to stop the design of this prototype and work on modeling of solar adsorption refrigerator 
machine. In this content, this work aimed to develop a neuronal model to predict the performance 
(Solar Coefficient of Performance (COPs)) of a solar adsorption refrigerator machine. This model 
is based on the development of a MatLab® program. The static neural network optimized 
obtained (SNNO) permitted to reproduce the experimental data from the literature with very good 
accuracy. The architecture of the SNNO includes (8-8-1) {8 neurons in the input layer, 8 neurons 
in the hidden layer and one output neuron (the solar coefficient of performance) } with an average 
absolute relative error of 0.45 %c and a correlation coefficient R > 0.9985 in the case of the test. 
To provide more flexibility in using SNNO, a designed friendly interface was developed using 
MatLab®. The sensitivity analysis of the input parameters on the output was also conducted in 
order to determine the parameters that have the most influence on the solar COP. This study 
allowed us to determine the solar COP by the inverse artificial neural network (IANN) method 
and to estimate the elapsed computing time by the IANN which is about 10 seconds. 

Keywords: Solar energy; refrigeration; artificial neural network direct and inverse; 
photovoltaic. 
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Introduction generate et problematique 


Introduction generate et objectifs : 

Les besoins energetiques de l'humanite n'ont cesse de croitre avec son evolution, 
1'homme a utilise des ressources qui sont apparues avec le progres scientifique telle que le 
charbon, le petrole, le gaz ou encore le nucleaire. 

Actuellement les energies fossiles (petrole, gaz, charbon, ....) sont consommees plus 
rapidement qu'elles ne se torment dans la nature, et il est estime que les reserves mondiales 
seront epuisees vers 2030 si la consommation n'est pas radicalement modifiee (Miller et al. 
2001; Kassmi et al. 2007). L’utilisation de ces energies fossiles est manifestee principalement 
par la pollution et le rechauffement climatique de la terre par effet de serre. Face a ces 
problemes, et afin de limiter l’emploi des energies fossiles, certains pays, se sont toumes vers 
l’emploi des energies renouvelables faisant appel de facon directe ou indirecte a l’energie 
solaire, (Poitiers. 2003). 

L'energie solaire est une source d'energie inepuisable et propre car elle ne provoque pas 
d'emissions de gaz nocifs pour notre environnement, elle peut etre utilisee de deux manieres 
differentes (Royer et al. 1998), soit par conversion photovoltaique, permettant de generer 
l’electricite directement a partir du rayonnement solaire, soit par conversion thermique, en 
transformant le rayonnement en chaleur ou en froid. De part sa situation geographique, 
l’Algerie dispose d’un des gisements solaires importants, car la duree d’ insolation sur la 
quasi-totalite du territoire national depasse les 2000 heures annuellement et atteint les 3900 
heures dans les hauts plateaux et le Sahara. 

Ce travail est motive par le fait que le besoin de production de froid dans les zones 
rurales ou les bases de vie et les postes avances sont prives d’electricite, en vue de la 
conservation d’ aliments, fruits, legumes et viandes, ainsi que les produits pharmaceutiques, 
medicaments et vaccins. C’est dans ce contexte que s’inscrit ce travail de these dans le cadre 
d’un projet commun entre le laboratoire de Biomateriaux et Phenomenes de 
Transport (LBMPT- Universite de Medea) et le laboratoire de Froid et Climatisation par 
Energies Renouvelables (LFCER-UDES), la partie experimentale a ete realisee au sein de ce 
dernier labo. 

Objectifs : 

Les objectifs du present travail consistent a etudier et realiser un refrigerateur 

domestique a compression de vapeur economique fonctionnant en energie solaire 

photovoltaique. D’ autre part une strategic basee sur les reseaux de neurones artificiels (RNA) 
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pour une modelisation simple, fiable et precise a ete elaboree pour estimer le coefficient de 
performance solaire d’une machine frigorifique a adsorption. 

En effet, les modeles mathematiques sont tres complexes a etablir, car le phenomene 
d’ adsorption/desorption est gouverne par une serie d’ equations aux derivees partielles non 
lineaires dont la resolution analytique est difficile, (Daniala et al. 2005 ; Mansouri et Ouali. 
2009). Les methodes numeriques peuvent estimer des parametres physiques avec des 
precisions et des vitesses de calcul meilleures, mais elles ne sont pas fiables quant a la 
precision de la solution analytique, (Daniala et al. 2005; Mansouri et Ouali. 2009). Or 
l’approche experimentale, bien qu’elle soit precise et privilegiee, elle presente l’inconvenient 
d’etre souvent delicate, couteuse et limitee par les moyens technologiques. Pour pallier ces 
contraintes, la modelisation neuronale peut etre une alternative aux methodes analytiques ou 


numeriques, (Daniala et al. 2005; Mansouri et Ouali. 2009). 


Les reseaux de neurones formels ou artificiels (RNA) suscitent toujours un vif interet 
dans la plupart des domaines d’ingenierie. Ils trouvent leurs applications, non seulement en 
matiere de modelisation statique mais aussi en modelisation dynamique des processus 
evolutifs dans le temps (Bouzidi et al. 2007 ; Si-Moussa et al. 2008). Une caracteristique 
fondamentale pour les reseaux de neurones artificiels, qu’ils sont capables d’approximer des 
relations complexes et fortement non lineaires entre un vecteur des inputs et un vecteur des 
outputs. Ils peuvent aussi traiter plusieurs systemes de refrigeration a la fois malgre la 
diversite des coordonnees geographiques, comme ils peuvent etre incorpores dans une boucle 
de controle, afin de determiner la tra§abilite du systeme c.a.d. la detection des susceptibles 
anomalies experimentales des inputs (Si-Moussa, 2010). 

La presente these s’organise en quatre chapitres: 

Apres une introduction generate, une synthese bibliographique a ete presentee dans le 
premier chapitre englobant quelques travaux concernant les refrigerateurs solaires 
photovoltaiques realises et la modelisation neuronale des refrigerateurs solaires a adsorption. 

Le deuxieme chapitre est consacre a la presentation des concepts de base du 
rayonnement solaire avec les conversions possibles de cette forme d’energie, ainsi que des 
generalites sur la refrigeration solaire photovoltaique et thermique. 


Dans le troisieme chapitre, apres une breve introduction et un apergu historique sur les 
reseaux de neurones artificiels, nous rappelons la notion de neurone formel et les fonctions 
d'activations, ensuite nous decrivons le perceptron et ses regies d'apprentissages. Par ailleurs 
nous presentons les perceptrons multicouches, l'algorithme de retro-propagation du gradient et 
ses principales variantes ainsi nous explicitons les principaux modeles de reseaux de neurones 
et les regies heuristiques necessaires au developpement d'un reseau neuronal. 
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Le quatrieme chapitre est scinde sur deux contributions : 

• la premiere contribution peut etre resumee comme suit : 

o L’ etude du dimensionnement et de la realisation de trois systemes 
frigorifiques photovoltaiques sous diverses conditions operatoires seront 
presentes. 

o L'optimisation du dimensionnement de ces trois systemes photovoltaiques 
sera effectuee avec le logiciel de dimensionnement de systemes et 
d'installations solaires photovoltaique PVSyst V5.59. 

• la deuxieme contribution peut etre recapitulee dans les points suivants : 

o L’optimisation des parametres de 1’ architecture neuronale directe et 
inverse d’un reseau de neurone statique pour 1’evaluation du coefficient de 
performance d’un systeme de refrigeration solaire a adsorption sera 
presentee, un programme sous MatLab® dans le but d’apporter plus de 
souplesse au reseau de neurones statique optimise (RNSO) sera elabore. 

o La comparaison, entre les valeurs du reseau de neurones statique optimise 
(RNSO) et les donnees experimentales sera realisee en terme d'erreurs 
dans le but de valider le RNSO et de montrer sa robustesse. 

o Une optimisation du reseau de neurone statique inverse sera conduite afin 
de trouver la tra§abilite du systeme global (les inputs en fonction des 
outputs), et par consequent mieux cernes le phenomene 
d’ adsorption/desorption et faire face aux anomalies au cours du processus 
experimental. 

o Pour apporter encore plus de souplesse a l’utilisation du RNSO, un 
programme informatique en langage MatLab® avec une interface 
conviviale sera developpe. 

o Les (RNSO et RNIO) ainsi developpes permettront de determiner le COP 
solaire et d’estimer le temps de calcul de la simulation. 

Les conclusions tirees de ce travail ainsi que les etudes qui peuvent etre envisagees dans 
le futur sont presentees a la fin de cette these. 
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Chapitre 1: Revue critique de la litterature 


1.1. Introduction 

Du fait des changements importants qui se sont produits durant ces dernieres decennies, le 
contexte energetique mondial connait aujourd’hui de grandes mutations. II appelle a une 
redefinition des besoins et des modes de consommation, car les risques d’epuisement des 
ressources non renouvelables tels que le petrole, le gaz naturel et le charbon largement 
explodes sont en voie de disparition. Le rechauffement climatique et les catastrophes 
nucleaires font apparaitre clairement que le modele energetique actuel n’est pas durable. II est 
done necessaire d’ engager une transition vers un modele qui serait plus durable pour les 
besoins des generations presentes et futures. 

1.2. Revue critique sur les refrigerateurs solaires photo voltaiques realises 

Pour la plupart des regions du monde, la mise a disposition de sources de froid est une 
base considerable de progres. Le besoin du froid pour la conservation des fruits, legumes et 
viandes, pour la conservation des produits pharmaceutiques et surtout des vaccins, est une 
necessite de premiere envergure. Pour subvenir a la demande en production de froid, 
nombreux systemes frigorifiques ont ete mis en pratique selon divers cycles 
thermodynamiques. Parmi ces cycles, le cycle a compression de vapeur, le cycle d’ absorption 
et le cycle a adsorption sont les plus utilises dans le domaine de la production du froid et de la 
climatisation. 

Les systemes de refrigeration les plus usites sont des systemes a compression de vapeur 
et fonctionnant aux fluides frigorigenes R12, R22, R502, etc. Ils necessitent de l’energie 
electrique pour leur fonctionnement. Dans les regions ou l’energie electrique est disponible, 
ce systeme de refrigeration est habituellement adequat pour satisfaire la plupart des exigences 
en refrigeration. Cependant, plusieurs regions ne disposent pas d’electricite. Ceci inclut des 
villages ruraux et des sites isoles dans les pays en voie de developpement. La refrigeration 
dans ces regions, par consequent, exige une solution alternative a ces systemes de 
refrigeration. L'utilisation de l’energie solaire comme source d’energie electrique pour 
1’ alimentation des systemes frigorifiques est tres interessante puisqu’elle est inepuisable et 
abondante dans les regions ou le besoin en froid est enorme. 

Plusieurs systemes ont ete realises dans ce contexte. Parmi ces systemes, les systemes 
frigorifiques photovoltaiques sont les plus appropries. On trouve dans la litterature plusieurs 
travaux qui montrent le progres de realisation des refrigerateurs domestiques, dans le but 
principalement de reduire la consommation electrique et proteger l’environnement, et cela soit 
par le changement des equipements frigorifiques tels que l’emploi des compresseurs qui 
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marchent en courant continu, evaporateurs eutectiques ou d’autres equipements ou soit par 
l’emploi des refrigerants qui ne participent pas a la degradation de l’environnement. De plus, 
F integration des energies renouvelables comme source de production d’electricite a aider 
d’ avoir des refrigerateurs autonomes, propres et moins chers. 

La recherche bibliographique effectuee dans cette partie est relative aux recents 
refrigerateurs domestiques a compression de vapeur realises jusqu’a 2013 et leurs 
alimentation en energie solaire photovoltaique, ces travaux sont resumes comme suit : 

Les premiers travaux relatifs a F alimentation de petits refrigerateurs alimentes par des 
generateurs solaires photovoltaiques realises par Tomachan et al. 1996 ; Thomachan et al. 
2000 ont prouve que les refrigerateurs domestiques conventionnels peuvent etre alimentes en 
courant altematif sinusoidal ou non-sinusoidal sans aucune degradation thermique des 
compresseurs hermetiques. 

Dans un autre laboratoire, Ulrich et al. 1999 ont etudie la performance de deux 
refrigerateurs alimentes en energie solaire photovoltaique autonomes sous les conditions 
climatiques ethiopiennes. Les resultats montrent la rentabilite des deux prototypes sous les 
conditions climatiques ethiopiennes. 

L’ etude effectuee par Bias et al. 2003 sur une installation frigorifique concue pour 
refroidir une production quotidienne prevue de 150 litres de lait a permis de montrer que 
l'autonomie du systeme est 2.5 jours s’il y a une formation de glace de 80%. L'energie 
photovoltaique obtenue pendant la joumee est stockee sous la forme de chaleur sensible et 
latente de l'eau gelee dans une cuve entourant un reservoir de lait. 

De leur cote, Socrates et al. 2006 ont choisi un refrigerateur conventionnel pour 
l’alimenter en energie solaire, pour cela quelques modifications ont ete apportees afin de 
reduire la charge de refroidissement et par consequent la puissance requise. 

Les etapes de fabrication, l'experimentation et la simulation menees par Anish en 2009 
et ses collaborateurs sur un refrigerateur domestique conventionnel converti pour fonctionner 
en energie solaire ont ete presentees. Divers tests ont ete menes pour etudier la performance 
du systeme. Ils ont arrive a reduire le cout des equipements principalement le champ et les 
batteries, et par consequent le cout global du systeme. 

De leur part, Petros et al. 2009 ont realise un congelateur solaire photovoltaique 
fonctionne sans stockage d’ energie pour etre utilise dans des applications veritablement 
autonomes dans les regions eloignees. Un nouveau concept de commande a ete etabli pour 
permettre la gestion de l'alimentation et le demarrage du compresseur, en outre le suivi de 
puissance maximale. Les resultats de prototype ont montre une tres bonne capacite de 
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fabrication de glace et un fonctionnement fiable et une grande amelioration dans les 
caracteristiques de demarrage des compresseurs, qui restent operationnelles meme pendant les 
jours a faible ensoleillement et de fonctionner avec une meilleure utilisation de l'energie 
photovoltaique disponible. 

Concemant le travail realise par Zhang et al. 2010 sur F adaptation d’un refrigerateur 
solaire photovoltaique pour fonctionner dans la region de Shanghai. Les resultats trouves 
montrent que la plus basse temperature dans le congelateur peut atteindre pres de -15 °C avec 
une consommation raisonnable de l’energie. Vu les performances du refrigerateur il peut 
servir pour le stockage de denrees perissables, des vaccins, des medicaments. 

Dans ce meme contexte, Ekren et al. 2011 et Ekren. 2012 ont montre la faisabilite 
d’alimenter un refrigerateur domestique fonctionnant avec un compresseur qui marche en 
courant continu via vitesse variable par energie solaire photovoltaique sans utilisation de 
l’onduleur, dans le but de reduire le cout du systeme parce que le prix de l’onduleur est eleve. 
La comparaison a montre que le fonctionnement a vitesse variable du compresseur a courant 
continu peut etre beaucoup plus efficace que le fonctionnement a vitesse constante, en 
particulier a des vitesses plus elevees. 

Dans le meme axe de recherche, Mehmet et al. 2011 ont realise un refrigerateur 
photovoltaique polyvalent pour etudier experimentalement ses performances au cours de 
fonctionnement joumalier et saisonnier sous les conditions climatiques de la Turquie. Les 
resultats revelent que dans l'ensemble, ce refrigerateur peut etre utilise de maniere fiable en 
permanence independamment du reseau local. 

Dans le but de la protection de l’environnement, Nawaz et al. 2012 ont conduit une 
comparaison entre un systeme de refrigeration solaire et conventionnel. D’apres les resultats 
trouves, le refrigerateur solaire s'avere etre techniquement meilleur, et aussi economique 
(lorsqu'il est utilise dans les petites industries et des usages domestiques, ou encore lorsqu’il y 
a des crises de l'electricite ou dans les regions de climats chauds), la quantite d'emissions de 
CO 2 a ete completement supprimee alors qu'il y a des emissions produites par le systeme 
conventionnel. 

A travers les resultats obtenus par Liu et al. 2012, sur F etude experimentale 
comparative a entre un refrigerateur solaire photovoltaique qui fonctionne a courant continu 
avec et sans Materiau a Changement de Phase (MCP). Le refrigerateur avec MCP presente de 
nombreux avantages par rapport a celui sans MCP. 
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A travers les resultats de modelisation de Mba et al. 2012, ils ont aborde une 
modelisation suivie d’une simulation sur le systeme, notamment l’effet des differents 
panneaux solaires sur la relation existante entre les temperatures de sortie du systeme de 
congelation et le rayonnement solaire incident. II est a noter que lorsque la tension du reseau 
augmente, la taille de la batterie de fonctionnement du systeme diminue et le courant 
correspondant a la puissance maximale augmente avec l'augmentation de l'insolation. II a ete 
egalement observe que la puissance MPP (Maximum Power Point) augmente quand 
l'insolation solaire incidente augmente, alors que la tension MPP reste constante (14.2V) et 
independamment du niveau de l'insolation. 


1.3. Revue critique sur l’estimation neuronale des parametres de controles des 
refrigerateurs solaires a adsorption 
1.3.1. Introduction 

La conversion solaire thermique connait un fort regain d'interet notamment sur les 
systemes de refrigeration a adsorption. Ces systemes sont particulierement adaptes aux 
regions ou 1’ ensoleillement est eleve et ou il n’y a pas de possibilite d’etre raccorde a un 
reseau electrique. Ceci est essentiellement du a trois facteurs, (Lemmini et al. 2002) : 

- Le rendement des panneaux photovoltaiques commercialises n’ augmente que faiblement 
depuis de nombreuses annees 0/ = JO a 15%). 

- Le champ des applications des systemes solaires thermiques s'elargit. 

- Les problemes environnementaux lies a 1’ utilisation de fluides frigorigenes (CFC et 
HCFC) contenus dans les circuits des systemes frigorifiques a compression. Ces fluides 
frigorigenes sont interdits a cause de leurs effets nefastes sur la couche d’ ozone. 

L’ utilisation des reseaux de neurones artificiels (RNA) pour la modelisation de la 
performance des machines frigorifiques a adsorption solaire est une solution rapide et plus 
simple, car les RNA sont dotes de proprietes de non linearite, d’ approximateurs universels et 
parcimonieux de fonctions. L’objectif de l’utilisation des RNA est pour but de reduire le 
risque d'incertitudes experimentales et aussi supprime la necessite aux equations analytiques 
complexes qui necessitent le temps de calcul et des efforts. La modelisation par les RNA est 
basee sur des donnees experimentales et/ou semi-experimentales (deduites des modeles 
analytiques, empiriques et numeriques, ...), done la qualite et la quantite de donnees decident 
de l'exactitude et la robustesse du modele. 
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1.3.2. Travaux realises 

Les RNA peuvent etre consideres comme un outil d’ approximation universel avec une 
capacite inherente d’extraire a partir des donnees experimentales les relations hautement non 
lineaires et complexes entre les variables decrivant le systeme etudie. Ils constituent 
maintenant un outil de calcul efficace et robuste en energie renouvelable. 

L’ application des reseaux de neurones trouvent plusieurs applications dans le domaine 
de la refrigeration solaire en generate. Mais dans cette revue critique 1’ accent est mis sur les 
applications des RNA pour la modelisation, dans le but de prediction et/ou de correlation, des 
proprietes de performance des refrigerateurs solaires a adsorption. Dans ce contexte peu de 
travaux ont ete publie dans la litterature (voir figure 1.1). 
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Figure 1.1: Revue critique sur les travaux publies sur la modelisation neuronale des 

refrigerateurs solaire a adsorption. 

Le premier travail relatif a 1’ application des RNA pour 1’ estimation de differents 
parametres des systemes de refrigeration solaire a adsorption conduit par Baiju et al. 2012a ou 
deux RNA avec 04 entrees et six sorties ont ete optimises. Dans le premier ils ont opte pour 
ralgorithme de Levenberg-Marquardt avec 9 neurones caches et pour le deuxieme, 
ralgorithme de Gradient conjugue avec 17 neurones caches ont ete utilises. La performance 
de prediction par le reseau de neurones artificiel statique optimise RNSO est obtenue en 
utilisant l’erreur quadratique moyenne, le coefficient de correlation et la covariance. Les 
coefficients de correlation et les erreurs quadratiques moyennes (%) sont de l'ordre de 
{( 0.9843 a 0.9971 ) et ( 0,590 % a 2,375%)} respectivement. 
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Dans le meme laboratoire, Baiju et al. 2012b ont teste pour la modelisation neuronale 
du lit adsorbant dans un systeme de refrigeration solaire a adsorption, l'algori thine de retro- 
propagation avec trois differentes variantes (le gradient conjugue, Pola-Ribiere gradient 
conjugue et Levenberg-Marquardt (L-M)) et la fonction de transfert (d’ activation) la sigmoide 

logistiques afin d’evaluer les parametres de performance du systeme ( ). Apres 

l’apprentissage, il est constate que 1'algorithme de L-M avec 14 neurones caches est le plus 
approprie. Les predictions par le RNSO de parametres de performance est en bon accord avec 
les valeurs experimentales. Le coefficient de correlation (R 2 ) est proche de 1 et le pourcentage 
d'erreur quadratique maximale est inferieur a 5%. 

Dans la meme epoque et dans le meme laboratoire, Baiju. 2013 a mobilise les 
performances d’une machine frigorifique a adsorption a double etage en utilisant les RNA. 
Apres la phase d’apprentissage, il a ete constate que 1'algorithme de LM avec 9 neurones est 
le plus approprie pour la modelisation de systeme de refrigeration adsorption solaire. Les 
predictions ANN de parametres de performance sont en bon accord avec les valeurs 
experimentales avec les valeurs de R 2 et pres de l’unite et un pourcentage d'erreur maximale 
inferieur a 5%. 

En travaillant sur la modelisation faite par Hernandez et al.2012 d’une machine 
frigorifique a adsorption par les reseaux de neurones artificiels directs et inverses qui 
fonctionnent avec le couple ammoniac/nitrate de lithium avec des differentes concentrations. 
Le RNSO direct a ete trouve avec une erreur minimale et un R 2 egale a 0.986, apres une etude 
preliminaire de sensibilite des parametres d’ entrees sur la sortie, le RNSO comporte 06 
neurones d’ entrees : la concentration de la solution, la temperature de l'eau de 
refroidissement, la temperature de generation, la temperature ambiante, la pression, le 
rayonnement solaire, 06 neurones caches et 01 neurone de sortie : le coefficient de 
performance de cette machine (COP). Le RNSO inverse a ete trouve avec grande precision et 
un temps de calcul court 0.29s, le COP a egalement ete determine en inversant le RNA pour 
calculer le parametre d'entree inconnu d'un COP requis. 

1.4. Conclusions 

L’utilisation de l’energie solaire pour la refrigeration et la climatisation est une 
application a promouvoir, d’autant que les besoins en froid coincident a la plupart du temps 
avec la disponibilite du rayonnement solaire (fonctionnement au fil du soleil). Le passage en 
revue de la litterature fait sortir principalement les points suivants : 

• Les systemes de refrigeration conventionnels necessitent de l'energie electrique pour 
leur fonctionnement, les sources de cette energie sont responsables a remission de 
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CO 2 et d’ autres gaz a effet de serre tels que les CFC et les HFC qui sont consideres la 
cause principale de l'appauvrissement de la couche d'ozone. 

• Ces problemes poussaient les chercheurs a trouver des autres sources d’electricite, tel 
que la conversion du rayonnement solaire par les panneaux photovoltaiques. 

• Dans ce contexte, plusieurs prototypes de refrigerateur domestique solaire a 
compression de vapeur ont ete realises, mais ils restent jusqu’ a ce jour plus chers avec 
un faible rendement. Chose qui obligeaient une autre fois, les chercheurs a apporter 
d’ autres ameliorations soit dans le systeme photovoltaique ou dans le cycle 
frigorifique. 

• Autres types de refrigeration solaire ont ete developpes, tel que, les refrigerateurs a 
adsorption qui atteignent une attention considerable dans les annees 1970 en raison de 
la crise energetique et les problemes ecologiques lies a l'utilisation des CFC et des 
HFC. 

• Les analyses thermodynamiques des systemes a adsorption sont complexes en raison 
des equations differentielles complexes impliquees. Au lieu de resoudre ces equations 
differentielles complexes et appliquer le nombre limite de donnees experimentales, des 
solutions rapides et plus simples peuvent etre obtenues a l'aide de reseaux de neurones 
artificiels. 

• La modelisation neuronale constitue une alternative prometteuse, ou du moins un outil 
complementaire, a la modelisation parametrique, puisque les modeles neuronaux sont 
dotes de proprietes de non linearite, d’ approximateurs universels et parcimonieux de 
fonctions. D’apres l’etude bibliographique citee, on peut observer que peu de travaux 
ont ete consacres a la prediction de la performance des systemes de refrigeration a 
adsorption solaire. 

• La methode neuronale, adoptee pour le calcul s’avere efficace. En effet, les differents 
tests effectues sur le modele neuronal concu pour tous les travaux presentes, ont 
permis de trouver des erreurs acceptables et un coefficient de correlation ou 
determination (R ou R 2 ) pratiquement egal a 1 . 
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2.1. Conversion de l’energie solaire 

2.1.1. Introduction 

Le soleil est la principale source des differentes formes d’ energies renouvelables 
disponibles sur terre. Le rayonnement solaire global est compose de rayonnements de 
differentes longueurs d’ondes et intensites. Le maximum d’intensite a lieu dans le 
rayonnement visible, soit pour une longueur d’onde de 0.5 pm. La repartition de l’energie 
solaire a la surface de la terre est comme suit (Barde. 1991): 

• Rayons ultra- violets ~ 6% 

• Rayons visibles ~ 50 % 

• Rayons infra-rouges « 44 % 

Un panneau solaire est un dispositif destine a recuperer une partie de l'energie du 
rayonnement solaire pour la convertir en une forme d'energie utilisable par 1'homme. On 
distingue deux types de panneaux solaires : les panneaux (capteurs) solaires thermiques qui 
convertissent la lumiere en chaleur recuperee et utilisee sous forme d'eau chaude ; les 
panneaux (modules) solaires photovoltaiques qui convertissent la lumiere en electricite, 
(Bernard et al.1980). 

2.1.2. Procedes de conversion et de stockage photovoltaique 
2.I.2.I. Principe de fonctionnement 

II existe differentes techniques permettant la conversion directe de la lumiere solaire en 
electricite, la plus connue est la conversion photovoltaique effectuee a l'aide de materiaux 
semi-conducteurs tel que le silicium (Si), le germanium (Ge), le selenium (Se) ou les 
composes semi-conducteurs tel que l'arseniure de gallium ( GaAs ), le tellurure de cadmium 
(i CdTe ). La majorite des cellules photovoltaiques sont fabriquees a partir du silicium cristallin, 
car il possede la caracteristique d'etre non toxique contrairement au cadmium ou au selenium, 
en plus, il permet d'atteindre des efficacites de conversion remarquables, il constitue environ 
28% de l'ecorce terrestre sous forme de composes (silicates, silice), ce qui en fait une source 
quasi inepuisable. Les cellules de type GaAs sont tres couteuses dans leur fabrication, leur 
utilisation est aujourd'hui essentiellement limitee aux applications spatiales (Francis. 2007). 

La cellule solaire ou photopile a semi-conducteur permette de debiter un courant 
electrique dans une charge externe lorsque celui-ci est expose a la lumiere. Son principe de 
fonctionnement se resume comme suit (Luc. 1995): 

a) Les photons de lumieres heurtent la surface photovoltaique disposee en cellule ou en 
couche mince 
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b) Us transferent leur energie aux atomes de la jonction, Si cette energie est suffisamment 
elevee, elle peut faire passer les electrons de la bande de valence a la bande de 
conduction du materiau semi-conducteur et creer ainsi des paires "electron- trou". 

c) Les electrons charges negatives (-) et les trous charges positives (+), sont alors 
maintenus separes par un champ electrique qui constitue une barriere de potentiel. 

d) Si une charge est placee aux bornes de la cellule, les electrons de la zone N rejoignent 
les trous de la zone P via la connexion exterieure, donnant naissance a une difference 
de potentiel et un courant electrique circule. 

e) Des fils metalliques tres fins connectes les uns aux autres recueillent alors le courant 
electrique continu qui est cree. Le courant est alors additionne d’une cellule a L autre 
jusqu’aux fiches de connexion du panneau photovoltaique et va ensuite s’additionner 
aux panneaux suivants connectes en serie, (Francis. 2007 ; Luc. 1995). 


Photon 



Charge 

electrique 

Exterieur 


Electron (-) 


Jonction 
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Zone N 


Figure 2.1: Principe de fonctionnement d’une cellule photovoltaique (Luc. 1995 ; 

Francis. 2007). 


2.I.2.2. Generateur photovoltaique 

Chaque cellule photovoltaique ne peut pas produire qu'une tension continue nominale 
de 0.5V a 0.6V en circuit ouvert et une puissance nominale voisine de 1.5W C . Pour satisfaire 
les besoins des charges couramment utilisees, il faut envisager un assemblage de plusieurs 
cellules photovoltaiques soit en serie soit en parallele. Cet assemblage forme ce qu'on appelle 
"module photovoltaique"(Luc. 1995). Un module compose de 36 cellules en silicium cri stall in 
est adapte pour la charge d'une batterie de 12V. Le dimensionnement du systeme solaire 
determine le nombre de modules a mettre soit en parallele afin d'accroitre le courant en 
conservant la tension, soit en serie afin d'augmenter la tension en conservant le courant, et 
pour avoir une satisfaction en courant et en tension, un groupement mixte "serie- parallele" est 
obligatoire, (Luc. 1995). 
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2.1.2.3. Configuration des installations photovoltaiques 

La Configuration des installations photovoltaiques se fait selon ces trois manieres, 
(Labouret et al. 2000) : 

a) Alimentations electriques faibles : Des modules photovoltaiques foumissent du 
courant continu pour n'importe quels appareils electroniques mobiles, alimentes par 
des piles ou batterie avec une sortie de 5V-20Vet une puissance max de 40W. 

b) Installations autonomes : comme les maisons en sites isoles, elles necessitent le plus 
souvent un stockage de l'electricite a l'aide d'accumulateurs. Si elles marchent en 
courant alternatif, il faut aj outer un onduleur. 

c) Installations ou centrales photovoltaiques connectees au reseau : Un generateur 
photovoltaique connecte au reseau n'a pas besoin de stockage d'energie. C'est en fait le 
reseau dans son ensemble qui sert de reservoir d'energie, (Labouret et al. 2000). 

2.1.2.4. Systeme de stockage (Installations autonomes) 

Dans le cas ou l'energie photovoltaique doit assurer la totalite des besoins en electricite, 
il est necessaire de la stocker pour les periodes non ensoleillees. Ce stockage est generalement 
assure par des batteries. La plupart des systemes photovoltaiques comportent des batteries 
speciales (batteries stationnaires a alliages de Plomb) qui accumulent l'excedent d'electricite 
produit puis le restituent en cas de besoin. La presence d'un regulateur est necessaire pour 
proteger les batteries contre les surcharges et les decharges profondes. Ces batteries sont 
concues pour restituer un courant stable pendant de longues periodes en conservant leurs 
aptitudes a la recharge, et ceci a un grand nombre de reprises (cycles), on parle de batteries 
stationnaires ou a decharge profonde, (Labouret et al. 2000). 

2.1.3. Potentiel solaire en Algerie 

De par sa situation geographique, l’Algerie dispose d’un des gisements solaire les plus 
eleves au monde. La duree d’ insolation sur la quasi-totalite du territoire national depasse les 
2000 heures annuellement et peut atteindre les 3900 heures (hauts plateaux et Sahara). 
L’energie re 5 ue quotidiennement sur une surface horizontale de lm 2 est de l'ordre de 5 kWh 
sur la majeure partie du territoire national, soit pres de 1700KWh/m 2 /an au Nord et 2263 
kWh/m 2 /an au Sud du pays, (Amardjia. 2007). 
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Tableau 2.1: Potentiel solaire en Algerie, (Amardjia. 2007). 


Regions (ce gisement solaire depasse 
les 5 milliards de GWh) 

Region 

cotiere 

Hauts 

Plateaux 

Sahara 

Superficie (%) 

4 

10 

86 

Duree Moyenne d’Ensoleillement 
(h/an) 

2650 

3000 

3500 

Energie Moyenne Re 5 ue (kWh/rn /an) 

1700 

1900 

2650 


2.1.4. Conclusion 

A la lumiere ce cette partie theorique, on peut citer les conclusions suivantes : 

• L’energie solaire est disponible partout sur terre. Notre planete re£oit 15000 fois l’energie 
que l’humanite consomme. 

• La reduction des couts de fabrication des systemes photovoltaiques reste une priorite a 
court et mo yen termes. Elle concerne en particulier la production de silicium de qualite 
‘solaire’, moins onereux que celui de qualite ‘electronique’. A plus long terme, de 
nouveaux materiaux pourraient succeder au silicium cristallin : silicium amorphe, 
materiaux organiques, ..., en particulier sous forme de couches minces. Par ailleurs, 
P amelioration de la partie conversion-gestion peut permettre de reduire les pertes et 
d’ameliorer la fiabilite des systemes photovoltaiques. 

• La necessite des systemes autonomes avec stockage de 1’ electricite de nature 
electrochimique, sous forme de batteries d’accumulateurs surtout pour les zones 
enclavees. 

• Le Programme National Algerian de Developpement des Energies Nouvelles et 
Renouvelables et de l’Efficacite Energetique pour la periode 2011-2030 a ete adopte par le 
Gouvemement en date du 3 fevrier 2011. 

• L’ Algerie ambitionne de produire 40 % de son electricite a l’horizon 2030 a partir des 
energies renouvelables et egalement de se positionner comme fournisseur majeur 
d’ electricite verte en direction du marche europeen en se fixant un objectif d’ exportation 
de 10000 MW en partenariat a la meme echeance. 

• Ce programme constituerait aussi le vecteur de developpement d’une industrie nationale 
des energies renouvelables qui s’appuiera sur les competences existantes en mettant en 
valeur F effort de recherche et de developpement dans les differents domaines lies a ces 
industries. 


- 14 - 




Chapitre 2: Conversion de Venergie solaire - refrigeration solaire 


2.2. Refrigeration a compression de vapeur 
2.2.1. Introduction 

Le choix d'une installation de refrigeration solaire se fait generalement en fonction des 
ressources disponibles dans le lieu ou l'on se trouve. II existe plusieurs types de refrigerateurs 
adaptes a differentes sources d'energies disponibles (tableau 2.2) (Webasto. 2003). 

Tableau 2.2: Differents types de refrigerateur selon les sources energetiques disponibles, 


(Webasto. 2003). 



Mode de fonctionnement 

Source d' energie 

Type de refrigerateur 

Refrigerateur 

solaire 

Conversion du rayonnement 
solaire en energie thermique 

Rayonnement 

solaire 

Adsorption 

Conversion du rayonnement 
solaire en energie electrique 
(panneaux photovoltai'ques) 

Rayonnement 

solaire 

Compression 

Conversion du rayonnement 
solaire en energie thermique 

Rayonnement 

solaire 

Absorption 


2.2.2. Principe de fonctionnement d’un cycle frigorifique 

Dans un premier temps, le compresseur va aspirer le gaz frigorigene a basse pression et 
a basse temperature (1). L'energie mecanique apportee par la compression va permettre 
F elevation de la pression et de la temperature du gaz frigorigene. Une augmentation 
d'enthalpie en resultera (2). La deuxieme etape represente le passage du fluide caloporteur 
dans le condenseur. Le gaz chaud provenant du compresseur va ceder sa chaleur au fluide 
exterieur et ainsi se refroidir "Desurchauffe 2-3" et cela avant l'apparition de la premiere 
goutte de liquide (3) zone I du condenseur. La condensation s'effectue jusqu'a la disparition de 
la demiere bulle de vapeur (4) zone II du condenseur. Le fluide liquide peut alors se refroidir 
de quelques degres ‘’sous-refroidissement 4-5” avant de quitter le condenseur (5) zone III du 
condenseur. La phase de condensation est en rapport directe avec la quantite de chaleur 
evacuee au niveau du condenseur. La troisieme etape, necessite F insertion d’un dispositif de 
detente detendeur. Afin d’abaisser la pression dans le circuit (6), pour permettre le 
changement de phase dans F evaporateur, du fluide caloporteur, de l’etat liquide a l’etat gaz. 
Le fluide frigorigene se vaporise partiellement dans l’organe d’ etranglement ou detendeur. La 
quatrieme et demiere etape, represente la production frigorifique du systeme au niveau de 
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1’ evaporateur. Le fluide frigorigene liquide entre en ebullition et s'evapore en absorbant la 
chaleur du fluide exterieur (7). Enfin, le gaz forme est encore legerement rechauffe par le 
fluide exterieur, c'est ce qu'on appelle la phase de ‘’surchauffe 7-1”. Le cycle est ferme, le 
fluide frigorigene evolue sous faction du compresseur dans les quatre elements constituant la 
machine frigorifique. L'ensemble du cycle peut etre represente dans le diagramme enthalpie- 
pression. Sous la courbe en cloche se situent les etats de melange liquide-vapeur; a gauche de 
la cloche, le fluide est a l'etat liquide "sous-refroidi", a droite, le fluide est a l'etat vapeur 
"surchauffe", (Figure 2.2), (Recknagel et al. 2001). 


Pression 

(P) 


Haute 
pression 


Basse 

pression 


Sous refroidissement 
du liquide 




Quail rite de h chaleur a. evacuee 
an eondenseur 


Desurchauffe de la 

vapeur 


Detente 



Compression 




Evaporation 


Production frigorifique de 
revaporateur 


, compresiion 


lauffe 


Enthalpie (H) 


Figure 2.2: Diagramme enthalpique du cycle frigorifique : diagramme de Mollier, 


(Recknagel. 2001). 


2.2.3. La surchauffe 

Un compresseur de fluide frigorigene peut transporter, en raison de sa construction, que 
des gaz ou de la vapeur. Les liquides ne peuvent pas etre comprimes et ne seraient par 
consequent pas figurer dans la chambre a compression du compresseur, (Recknagel. 2001 ; 
Rapin. 1996). 

Si le niveau d’ aspiration du compresseur est situe directement sur la ligne de rosee, une 
diminution de charge de f evaporation peut entrainer une “aspiration“ humide. II peut en 
resulter deux consequences nefastes pour le compresseur. Le fluide frigorifere liquide lave le 
film lubrifiant entre le piston et les parois du cylindre; une lubrification insuffisante provoque 
une usure plus importante. Si le fluide frigorigene liquide parvenait dans le cylindre, celui-ci 
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entrainerait, lors de la compression, un transport direct de l’energie du piston a la culasse du 
cylindre. La plaque de soupape peut etre endommagee par les coups de liquide. 

Afin d’eviter toute aspiration de liquide, l’etat d’ aspiration du compresseur est decale de 
la ligne de rosee vers la droite. On "surchauffe" la vapeur du fluide frigorigene. 

Sa temperature t 0 2 h est situee au-dessus de la temperature T a . La surchauffe est calculee 
de la maniere suivante (Rapin. 1996) : 

^02h =r ^02h~ T 0 ( 2 . 1 ) 

• AT 0 2 h : Surchauffe a la sortie de l’evaporateur en °K; 

• T Q 2 h : Temperature de fluide frigorigene a la sortie de l’evaporateur en °C; 

• T a : Temperature d’ evaporation en °C, la lettre h est utilisee pour la surchauffe 
(Rapin. 1996). 



Figure 2.3: Surchauffe du fluide frigorigene dans 1’ evaporateur, (Rapin. 1996). 

Outre la protection du compresseur contre les liquides, la surchauffe presente d’autres 
avantages. Etant donne que le pourcentage de liquide dans la vapeur d’ aspiration entraine une 
reduction du debit de transport du compresseur, ce pourcentage pourra etre augmente par des 
phases de chauffage precises. La surchauffe ameliore notamment le recyclage de l’huile 
(Rapin. 1996). 

La valeur de surchauffe usuelle atteint environ 5-8°K, le rendement maximal de 
T installation est atteint lors de cette surchauffe, (Recknagel et al. 2001 ; Rapin. 1996). 

2.2.4. Le sous -refroidissement 

L’organe de detente a pour fonction de reduire, apres la condensation, le fluide 
frigorigene a un niveau de pression le plus has (pression d’ evaporation). La soupape 
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fonctionne de facon optimale lorsque le liquide presente a 1’ entree un niveau de purete 
egalement optimal, (Recknagel et al. 2001 ; Rapin. 1996). 

Si l’on compare une certaine masse de fluide frigorigene a l’etat liquide et a l’etat de 
vapeur (a pression constante), le fluide frigorifique sous forme gazeuse requiert un volume 
bien plus important. H s’ensuit que le fluide frigorigene sous forme de vapeur necessite plus 
de temps pour se flayer un passage a travers "l’etranglement", (Recknagel et al. 2001 ; Rapin. 
1996). 


Le fluide frigorigene sous forme gazeuse reduit, avant de parvenir a la soupape de 
detente, le debit et entraine une sous-alimentation de l’evaporateur en fluide frigorigene. 

Si l’installation frigorifique est exploitee de telle sorte que l’etat “Entree detendeur“ se 
trouve directement sur la courbe de separation gauche (ligne d’ebullition), les moindres 
variations des conditions d’ exploitation peuvent provoquer la formation de bulles a 1’ avant du 
detendeur (Recknagel et al. 2001 ; Rapin. 1996). 

Pour cette raison, l’etat “entree detendeur“ est deplace de la ligne d’ebullition vers la 
zone liquide et on parle alors de sous refroidissement. qui garantit une alimentation en liquide 
a 1’ avant de la soupape de detente. Le sous-refroidissement est ainsi calcule : 


AT , =T -T , 

c2u c c2u 


( 2 . 2 ) 


• AT c2u : sous-refroidissement a la sortie du condenseur en °K 

• T C 2 u '■ Temperature du FF a la sortie du condenseur en °C 

• T c : Temperature de condensation en °C 



Figure 2.4: Sous-refroidissement du fluide frigorigene. 
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Le Sous-refroidissement optimal serait la plupart du temps situe entre 4-7°K, 
(Recknagel et al. 2001). 

2.2.5. Dimensionnement du systeme de production frigorifique de la chambre positive 

Le but de ce travail est de dimensionner les differentes parties d’une chambre froide 
positive pour repondre aux besoins de conservation des produits pharmaceutiques, 
halieutiques et alimentaires en vue d’une consommation ulterieure pour une duree limitee. 
Une etude thermique detaillee des quantites de chaleur a extraire de l'interieur d'une chambre 
froide permettra de determiner la puissance frigorifique qui sera la base de dimensionnement. 
(voir figure 2.5). 

Ces quantites de chaleur sont calculees sur 24 h. Les dimensions de la chambre a 
concevoir sont resumees dans le tableau suivant: 


Tableau 2.3: Dimensions de la chambre froide positive 


Volume 
total (m 3 ) 

Dimensions exterieures (mm) 

Dimensions interieures (mm) 

1.56 

Longueur 

Largeur 

Hauteur 

Longueur 

Largeur 

Hauteur 

2000 

1295 

1600 

1210 

1055 

1220 


La fabrication et 1’ installation des chambres froides repondent a des normes de securite 
et d’hygiene. Les normes en vigueur sont NF C 15-100 pour l’installation electrique et NF E 
35-400 pour F installation frigorifique. 

Les calculs d’une chambre froide doivent satisfaire a trois conditions suivant le produit 
a traiter: 

• L’hygrometrie 

• La ventilation 

• La temperature en chambre de refrigeration entre 0°C et 8°C. 

Le groupe frigorifique sera adapte aux conditions du sud algerien (Temperature 
ambiante avoisinant de 50°C) et la puissance frigorifique du groupe doit etre adaptee au 
volume de la chambre pour avoir un froid positif compris entre 0 et +5°C. 

2.2.6. Caracteristiques techniques (materiaux) 

• Isolation des murs exterieurs et plafond en panneaux sandwich polyurethane 

• Isolation du sol en panneaux sandwich polyurethane renforcee 
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• Structure metallique exterieure en comieres 

• Porte iso thermique pivotante simple avec rideaux laniere 

La chambre contient un espace pour equipements frigorifiques et electriques tel que: le 
compresseur, condenseur a air, evaporateur a air, armoire electrique de commande et les 
organes annexes. Le choix de polyurethane comme isolant est du a ses caracteristiques 
interessantes (stable dans le temps, leger et resistance mecanique moyenne) et son cout 
abordable. L’ isolation doit limiter le coefficient global de transmission thermique a K=0.36 
W/(m 2 .°C). 



Figure 2.5: Vue exterieure de la chambre a realiser 

2.2.7. Isolation des parois verticales et du plafond 

La construction de 1’ isolation des parois se fait en parois sandwich constitues de : 

• Panneaux sandwich en mousse de polyurethanne 40 kg/m 3 

o Revetement interieur en tole inoxydable 
o Revetement exterieur en tole laque verte 

• L’ossature exterieure se compose de cornieres 100x100 mm en acier peint en vert, 
fixees sur les coins exterieurs des panneaux sandwich par des boulons a ecroues. 

• Porte iso thermique (700x1100) mm, equipee avec panneaux sandwich en mousse de 
polyurethane de 40 kg/m 3 de 120 mm d’epaisseur, les joints d’etancheite en 
caoutchouc. 
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La plancher (panneaux de sol) est constitue de : 

• Panneaux sandwich en mousse de polyurethanne 40 kg/m 3 et une conductivite 
thermique : A i SO iant = 0.027 W/m.°C. 

• Revetement interieure en resine antiderapante sur une plaque en contreplaque (10 
mm). 

• Revetement exterieur du plancher est en tole d’acier galvanise e ag =1.8 mm fixe sur la 
structure metallique de l’ossature de la chambre. 

L’epaisseur de la mousse a injecter doit etre comme suit, (Breidert. 1995) : 


K 


J_ 

1 e ag_int e isolan t | e ag_ext 1 
^int ^ ag_int ^isolant ^ ag_ext h ex £ 


= 0.36 


w 


m 2 .°C 


(2.3) 


Ou : 


Q s 1 BXt 

^isolant ~ ^isolant 36 — 


* 

a 9-int 


hi 


A, 


) 


(2.4) 


■a-dext hext ,l int ' La 9int 

Pour une paroi de chambre froide positive on adopte : l/hi nt =0.12 m 2 . °C/W et 
l/hext=0.03 m 2 .°C/W, li u „ et h ext sont des coefficients de convection interieur et exterieur. 

Toute machine frigorifique est constituee de differents composants a savoir dans le cas 
d’une machine frigorifique classique a compression d’une vapeur avec changement de phase : 
un compresseur, un condenseur, un detendeur, un evaporateur, un reseau de tuyauteries et 
enfin un certain nombre d’organes de commande et de regulation. 


entree de fair 



sortie de fair 


Condensateur a 
refroidissement 
a air 


Detendeur 


Compresseu 


M 


c6te haute pression 
de la machine fngonfique 




c6te basse pression 
de la machine fngonfique 


iitiiMimz 


appcrts dus 
aux personnes 







0 o''» <JU L 





entree de rair 


Evaporateur 


sortie de Pair 


apports dus 
£ I'dclairage 



Chambre froide 




- apports dus au renouvellement d air 

- apports dus & I'cuverture des portes 


Figure 2.6: Schema de principe d’une installation frigorifique utilisee en refroidissement 

d’une chambre froide (Breidert. 1995). 
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Figure 2.7: Methodologie de calcul des besoins energetiques des elements principaux d’une chambre froide positive 


Chapitre 2: Conversion de Venergie solaire - refrigeration solaire 

2.2.8. Bilan thermique d’une chambre froide 

Les charges thermiques autour d’une chambre froide se resument en deux composantes : 
les charges internes et externes, la figure 2.7 resume la methodologie de calcul des besoins 
energetiques des elements principaux d’une chambre froide positive. 

a. Charges externes: elles sont constituees essentiellement de : 

• Charge thermique due aux apports de chaleur par transmission a travers l’enveloppe de 
chambre froide (parois verticale, planchez bas, planchez haut). 

• Complement d’apport par transmission a travers les parois exposees directement au 
rayonnement solaire 

• Charge thermique par renouvellement d’air 

• Charge thermique par ouverture des portes. 

b. Charges internes: elles sont subdivisees en deux categories: 

• Charges thermiques independantes des denrees entreposees 

o Charge thermique due a 1’ eclairage 
o Charge thermique due aux personnes 
o Charge thermique due a des machines diverses. 
o Charge thermique due au materiel roulant 

• Charges thermiques internes dependantes des denrees entreposees (respiration des fruits 
et des legumes, etc. . .). 

o Charge thermique due aux denrees entrantes 
o Charge thermique due a la respiration des denrees 
o Charge thermique due aux moteurs des ventilateurs 
o Charge thermique due aux resistances de degivrage 

2.2.9. Determination des charges externes 

a. Charge thermique par transmission a travers les parois 

Les charges thermiques dues aux apports de chaleur par transmission sont: 

• Le plafond 

• Le planchez ou sol (mais par hypothese, on considere qu’il n’y pas d’apport) 

• Les deux parois verticales orientees Nord et Sud. 

• Les deux parois verticales orientees Quest et Est 
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Soient S mt et S ext les surfaces interieures et exterieures de ces differentes parois. Les calculs 
de surfaces d’echange sont donnes dans le tableau 2.4: 


Tableau 2.4: Calcul des surfaces interieures et exterieures 


Parois 

Surfaces interieures 

Surfaces exterieures 

Plafond 

1.210x1.055 = 1.28 

2.0x1.295 = 2.95 

Paroi Nord 

1.220x1.210= 1.48 

1.6x2.000 = 3.20 

Paroi Sud 

1.220x1.210= 1.48 

1.6x2.000 = 3.20 

Paroi Ouest 

1.055x1.220= 1.287 

1.6x1.295 = 2.07 

Paroi Est 

1.055x1.220= 1.287 

1.6x1.295 = 2.07 

Somme des surfaces 

S int = 6.814 m 2 

S ext = 7.090 m 2 


Les charges dues aux apports de chaleur par transmission sont donnees par la relation 
suivante (Breidert. 1995) : 


Qtp = 24 x K x 5 x AT, (Wh) 


(2.5) 


Avec K est le coefficient global d’echange de chaleur (W/m 2 .K) donnee par la formule 2.3, 
S Aire de la paroi a travers laquelle se fait le transfert (m 2 ), AT ecart de temperature en °K 
propre a chaque paroi entre la temperature ambiante (T a ) et la temperature de la chambre froide 
(Tf) et 24 est la duree d’une journee en heures. 

Les parois des chambres froides etant souvent bien isolees thermiquement, par souci de 
simplification, on ne tient compte dans les calculs, que la resistance thermique de l’isolant 
(Breidert. 1995). 


K - 


1 | e isolant | 1 

h-int ^ isolant ^ ext 


( 2 . 6 ) 


Concernant la surface de transmission est donnee par la formule suivante: 


S V ^int^ext 


(2.7) 


b. Complement d’apport par transmission a travers les parois exposees directement au 
rayonnement solaire 

Comme la chambre sera placee dans un endroit a l’abri de la lumiere, le complement 
d’apport par transmission a travers les parois exposees directement au rayonnement solaire est 
suppose nul (Breidert. 1995). 
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c. Charge thermique par ouverture des portes 

Dans une chambre froide comportant une seule porte on calcule simplement la charge 
thermique par renouvellement d’air. Mais si la chambre contient plusieurs portes il faut calculer 
la charge thermique due a l’ouverture des portes. Dans notre cas on a une porte, done Q op =0, 
(Ndoye. 2002 ; Breidert. 1995). 


d. Charge thermique par renouvelement d’air 

Dans de nombreuses chambres froides, il est prevu de renouveler plus ou moins Fair 
ambiant, e’est-a-dire de remplacer une partie de Fair de la chambre froide par de Fair exterieur. 
La quantite d’air neuf admise doit etre refroidie de la temperature exterieure a la temperature de 
la chambre froide constitue done une charge thermique. La Charge thermique par renouvelement 

d’air (Q re en kW (~)) est donnee par F equation suivante (Breidert. 1995): 

Qre TH-a.e- ( 2 . 8 ) 

Avec m a . e : Debit-masse d’air exterieur admis en kg/s, Ah : difference d’enthalpie entre 
Fair exterieur h a . e (T a i r -ext=35 °C avec un taux d’humidite egale a 70%, done h a . e = HO kJ/kg 
determine a partir du digramme de Fair humide presente dans la figure 2.8) et Fair ambiant h a . a 
(T a ir-int=4°C avec un taux d’humidite egale a 80%, done h aM =17.5 kJ/kg determine a partir du 
digramme de Fair humide presente dans la figure 2.8) de la chambre froide en kJ/kg. Pour des 
conditions de temperature et d’humidite relative de Fair donnee, l’enthalpie correspondante peut 
etre lue soit dans un diagramme de Fair humide (Figure 2.8) soit dans des tables specialises. 


On a par ailleurs (m a , e e n kg/s) (Breidert. 1995): 


m 


a.e 


Va.e-Pa.a 

86400 


(2.9) 


Avec 86400 nombre de secondes dans une journee, V ae debit- volume d’air exterieur en 
m 3 /s, p a .a masse volumique de Fair de la chambre froide en kg/m 3 (Breidert. 1995) : 

Pa.a ~ ~T~r ( 2 - 10 ) 

1 + 

273.15 

La valeur de po etant prise egale a 1.293 kg/m 3 , T etant la temperature de la chambre froide 
exprimee en °C. Enfin, pour le debit-volume d’air exterieur (renouvellement d’air), on a : 


Vz.e Vcf-ft 


( 2 . 11 ) 
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Avec V c f volume de la chambre en m , n en s' 1 est le taux de renouvellement d’air 
journalier donne par l’equation suivante (Breidert. 1995): 



Figure 2.8: Diagramme de l’air humide utilise pour le calcul des variations d’enthalpie de 

l’air, (Breidert. 1995). 


2.2.10. Determination des charges internes 

2.2.10.1. Charges thermiques independantes des denrees entreposees 
a. Charges thermique due a l’eclairage 

Le niveau d'eclairement moyen a atteindre dans les lieux de stockage est de 125 a 250 lux. 
Les fabricants prevoient, en general, une puissance de 10 W/m 2 avec des lampes 
d’ incandescence. La charge thermique apportee par les lampes est donnee par (Breidert. 1995 ; 
Ndoye. 2002): 

Q ec = 860. VK.t.n (2.13) 

Avec W est la puissance unitaire des lampes en kW est de 10 W/m 2 x2.95m 2 =29.5 W. et n 
est le nombre des lampes, t est le temps d’eclairage en heures par jour (24h). 
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b. Charge thermique due aux personnes 

La charge thermique due aux personnes se calcule par la formule (Ndoye. 2002): 

Q P = qtn-t-N (2.14) 

Avec N est le nombre de personne, q th est la puissance totale emise par individu en kW 
(sensible + latente), cette valeur est de l’ordre de 372 W/personne et t est le temps de sejour en 
heures (lh). 


2.2.10.2. Charges thermiques internes dependantes des denrees entreposees 
a. Charge thermique due aux denrees entrantes 

L’eau a ete utilisee comme charge avec un volume maximale de 140 Kg et avec une 
temperature de 30°C avant introduction dans la chambre froide (doit refroidir a 15 °C) (Ndoye. 
2002 ). 

Qde = Qdel + Qd2 (2.15) 

Avec Qdei refroidissement entre la temperature d’ entree T e et de la congelation T c , Q d2 
charge due aux emballages (8% Q de i )• 

Qdel = m C (Je — T c ) (2.16) 

Avec m masse des denrees, C chaleur denrees avant congelation donnee par la formule 
suivante : 


~ _ a+OAb 
~ 100 

Avec a teneur en eau de la denree (%), b la teneur en matieres solide (%). 


(2.17) 


2.2.11. Puissance frigorifique 

Cette etude a permis de calculer les charges susceptibles d’etre prise en compte dans les 
calculs thermiques de la chambre. Ces charges etant tres difficiles a determiner proprement sans 
faire des hypotheses, le tableau suivant resume les differentes charges calculees. 


La duree de fonctionnement de compresseur est 12 heures sur 24 heures (Ndoye. 2002). 



charge quotidienne 

duree de fonctionnement de compresseur 


0.46 kW 


(2.18) 
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Tableau 2.5: Differentes charges thermiques en (kWh/j) calculees sur la chambre froide 

positive 


Charge thermique par transmission a travers les parois 

2.28 

Charge thermique par ouverture des portes 

0 

Charge thermique par renouvelement d’air 

0.123 

Charge thermique dues a l’eclairage 

0.61 

Charge thermique dues aux personnes 

0.73 

Charge thermique dues aux denrees entrant 

1.75 

Apport majoree avec 5% 

0.265 

Somme des charges quotidienne 

5.76 


2.2.12. Dimensionnement des differents elements de la chambre froide 
a. Determination du compresseur 

Le cycle frigorifique a ete trace sur un diagramme enthalpie h=log(p ) en basant sur la 
valeur de la production frigorifique nette, la temperature d’ evaporation et de la condensation du 
R12, la valeur de la surchauffe dans l’evaporateur et dans la tuyauterie d’ aspiration et la valeur 
de sous refroidissement, en suite on a determine le taux de compression et l’exposant poly 
tropique pour calculer les temperatures a 1’ entree et a la sortie de chaque point de cycle. 

En utilisant le cycle trace on a determine premierement le debit massique a deplacer (eq. 
2.19), debit volumique qui doit aspirer le compresseur (eq. 2.20), debit theorique du compresseur 
(eq. 2.21-2.22), dimensions du compresseur (eq. 2.23), calcul de la puissance effective sur 
l’arbre du compresseur (eq. 2.24), calcul de la puissance du moteur d’entrainement (eq. 2.25), 
calcul des conduites (eq. 2.26), (Breidert. 1995). 



Qo 

Ahev 


Avec Qo est production frigorifique nette et A h ev 
evaporateur 


(2.19) 


difference d’ enthalpie entre 


Qvr QmPl ( 2 . 20 ) 

Avec vi est le volume massique d’air exterieur (m 3 /kg) 

Qvm = \ r (2.2i) 

'IV 


Avec t] v est le rendement volumetrique donnee par : 

rj v — 1 — 0.05.— (2.22) 

Pa 

Avec P a et P r sont les pressions d’ aspiration et de refoulement respectivement 
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q v m = ^.i.ve ff -z. ^.icr 9 (2.23) 

Avec d diametre de piston en mm, l la course du piston en mm, rj e ff nombre d’effets du 
compresseur, Z nombre de cylindre et N vitesse de rotation en tr/min. on prend comme 
parametres de conception rj e ff=l, Z=1 avec une vitesse de rotation de 3000 trs/min. 

p = Qm^zhi) (2.24) 

Vm 

Avec r] m rendement mecanique du compresseur (90%) et Ah a travers le compresseur 

P m — 13.P e ff (2.25) 

A t = (2.26) 

w a 

Avec Wi la vitesse du fluide dans la conduite. 


b. Determination de l’evaporateur 

La surface d’echange thermique S ev de l’evaporateur ventile de marque FRIGA BOHN est 
donnee par la formule suivante, (Ndoye. 2002): 



Qpf 


K gy.AT 


(2.27) 


Avec K ev est le coefficient global de transmission thermique de l’evaporateur en W/m 2 .°C, 
AT est la difference de temperature entre la temperature moyenne du milieu a refroidir et la 
temperature d’ evaporation du fluide frigorigene. On adapte pour notre cas un evaporateur a gaz a 

circulation d’air avec des tubes a ailettes (K PV = 47 ). 

v ev m 2 .°C 


c. Determination du condenseur 

Determination de la surface d’echange (Ndoye. 2002): 

c Qcond 

cond ~ Kfc.AT 


(2.28) 


Avec Qcond puissance calorifique cedee au condenseur, K k coefficient d’echange du 
condenseur variant de 20 a 30 W/m 2 .°C. AT est la difference de temperature entre la temperature 
de condensation et la temperature a 1’ entree du condenseur. 

La selection du detendeur s’effectue en fonction de la puissance frigorifique et le debit 
volumique du liquide. 
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Tableau 2.6: Dimensionnement des differents elements du systeme frigorifique 


Surface d’echange thermique de l’evaporateur 

1.33 m 2 

Surface d’echange thermique du condenseur 

0.5 m 2 

Puissance du compresseur hermetique 

3/8 cv (275 W) 

Detendeur 

11.84 bars 



Figure 2.9: Differents equipements constituant la chambre froide 
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Avec 1: mousse de polyurethane de 40 Kg/m3 epaisseur 120mm ; 2 revetement interieur en 
resine antiderapante sur une plaque en Contreplaque (10 mm) ; 3 porte iso thermique (700x1100) 
mm ; 4,5 et 6 mousse de polyurethane avec de 40 Kg/m3 epaisseur 120mm avec revetement 
interieur en tole inoxydable et revetement exterieur en tole laque vert ; 8 evaporateur ventile 
positif de marque FRIGA BOHN ; 9 cache en grille d’ aeration en tole d’acier galvanise peint en 
vert. 10 et 19 groupe electrogene et systeme solaire a prevoir ; 11 groupe de condensation avec 
compresseur hermetique TECUMSEH ; 12 ; 13 et 20 ossature exterieur ; 14 et 16 structure 

metallique ; 15 revetement interieur de plancher est en tole galvanise 18/ 10mm; 17 ; 18 tableau 
de commande de regulation DIXELL. 

2.2.13. Conclusion 

II ressort de cette etude theorique que la refrigeration a compression de vapeur est le 
procede le plus repandu pour la production du froid. Les systemes de refrigeration 
conventionnels necessitent de l'energie electrique pour leur fonctionnement, les sources de cette 
energie sont responsables a remission de C02 et d’autres gaz a effet de serre tels que les CFC et 
les HFC qui sont consideres la cause principale de l'appauvrissement de la couche d'ozone. Ces 
problemes poussaient les chercheurs a trouver des autres sources d’ electricite, tel que la 
conversion du rayonnement solaire par les panneaux photo voltaiques. 

Dans cette partie, une etape de dimensionnement des differents elements de la chambre 
froide positive de 1.56 m 3 a ete faite. La realisation de cette chambre a ete faite dans une usine de 
fabrication des chambres froides a Tizi-ouzou. Cette chambre est adaptee pour etre alimentee 
avec une source d’energie renouvelable dont l’impact positif sur 1’ environnement est 
incontestable (energie solaire). 
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2.3. Refrigeration solaire a adsorption solide-gaz 

2.3.1. Introduction 

Les machines frigorifiques a adsorption font partie des solutions alternatives pour 
remplacer les systemes de refrigeration et de climatisation conventionnels. Les 
refrigerateurs solaires a adsorption utilisent la faculte de certains materiaux microporeux a 
fixer des fluides a l'interieur de leurs pores. Beaucoup de materiaux poreux tels que le 
charbon actif, les silica gels et les zeolithes satisferont a cette condition. 

Les couples principalement utilises sont le charbon actif associe au methanol, le 
silica gel et l’eau ou encore la zeolithe et l’eau. 

2.3.2. Description du modele de base 

La figure (2.10) montre le schema de base de refrigerateur solaire a adsorption qui 
est compose essentiellement de (Jannot. 2003). 

• Un capteur solaire constituant le generateur; celui-ci est rempli d'adsorbant 
chauffe directement par les rayons du soleil ; 

• Un condenseur a eau ou a air ; 

• Un vase de recuperation de F adsorbat condense ; 

• Un evaporateur contenant le refrigerant (l'adsorbat) et qui est place dans une 
caisse isolante constituant la chambre froide ; 

• Des vannes permettant d'isoler les differents elements et d'assurer le deroulement 
du cycle de refrigeration. 



Figure 2.10: Schema de base d’un refrigerateur solaire a adsorption, (Jannot. 2003). 
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2.3.3. Principe de fonctionnement 

Le refrigerateur solaire a adsorption fonctionne selon un cycle intermittent jour/nuit 
alternant des phases de desorption/condensation (diurnes) et des phases d’ adsorption/ 
evaporation (nocturnes). Pendant le jour, le generateur est en communication avec le 
condenseur. Les rayons du soleil chauffent l’adsorbant et l’adsorbat qu’il contient, la 
temperature et la pression augmentent progressivement dans le capteur. Lorsque cette 
derniere atteint la pression de saturation correspondant a la temperature du condenseur, la 
desorption commence et les vapeurs desorbees se condensent et sont acheminees par 
gravite vers le vase de recuperation. L’adsorbant a joue le role de compresseur chimique, 
(Duffie et al. 1991). 

En fin d’ apres-midi, la temperature et la pression baissent dans le generateur, le 
refrigerant contenu dans le vase de recuperation est achemine vers 1’ evaporateur. Ce 
dernier est mis en communication avec le generateur. La temperature et la pression 
continuent de baisser. Lorsque la pression atteint la valeur de saturation correspondant a la 
temperature de 1’ evaporateur, l’adsorption commence avec l’evaporation qui produit le 
froid necessaire. L’adsorbant a joue le role de detendeur chimique. Cette phase se prolonge 
jusqu’au matin (Duffie et al. 1991). 



9 8 7 


Figure 2.11: Schema de la machine frigorifique solaire a adsorption 
l-Vitrage ; 2-Adsorbeur ; 3- Boitier du capteur ; 4- Collecteur ; 5- Support ; 6 - Chambre froide ; 7- 

Evaporateur; 8 - Condenseur ; 9- Volets de refroidissement. 


2.3.4. Coefficient de performance (COP) 

Le coefficient global joumalier de performance des systemes a adsorption se calcule 


au moyen de la formule suivante, (Hildbrand et al. 2004): 


COR 


Energie frigorifique (Q EV ) 


solaire — 


Energie solaire globale journaliere refue sur la surface des modules (G) 


COP , 


m ad x Am (kg de fluide/kg d'adsorbant) x {L v (T evap ) - C p (T cond 

^evap)} 


solaire ~ 


2 X S cav X f 

CUfJ Jy 


heure de coucher du soleil 


heure de lever du soleil 


G(t). dt 


(2.29) 

(2.30) 
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Avec Am la quantite de fluide desorbe par unite de masse de l’adsorbant (kg de 
ri ui de/kg d’adsorbantX • masse d adsorbant. 

II est tres difficile de donner un COP s defini pour chaque systeme. En effet, le COP s 
evolue grandement en fonction de la temperature exterieure journaliere moyenne et de 
l'irradiation. Le COP s moyen des systemes publies dans la litterature se situe entre 0.083 
et 0.21, (Hildbrand et al. 2004). 


2.3.5. Dimensionnement de la machine frigorifique 

Actuellement, la technologic nous offre une variete de moyens permettant d’obtenir 
le froid. Ainsi, de nombreux scientifiques et ecologistes se prononcent done pour une 
technologic capable d’ assurer le futur ecologique de notre planete. La disponibilite et la 
gratuite de l’energie solaire, d’une part, et la protection de 1’ environnement, d’ autre part 
pourraient rendre les machines a adsorption de plus en plus competitives dans le marche 
relatif a cette industrie. Comme l’Algerie dispose d’un gisement solaire important, il est 
interessant d’ exploiter cette ressource naturelle dans le domaine de la production de froid 
surtout dans les zones rurales isolees. C’est dans ce concept que se situe notre travail dont 
le but est de dimensionner afin de construire une machine frigorifique solaire. Les etapes 
de dimensionnement ont ete prises principalement a partir de (Chekirou. 2008 ; Fadar. 
2008) selon F organigramme presente dans la figure 2.12. 


2.3.6. Calcul des besoins en eau glacee (stockage) 

Pour determiner la quantite d’eau glacee, que la machine frigorifique doit produire, il 
faut calculer le bilan journalier du local considere. La masse totale d’eau est calculee a 
partir de la formule suivante : 


m e = 


Q 


Fj 


C Petise -he) 


(2.31) 


Ou Qfj : charge frigorifique journaliere (calculee a partir de la batterie) 

La machine doit refroidir l’eau d’une temperature initiale jusqu’a la temperature 
voulue a 1’ entree de la batterie (charge) humide : 

Q fj =me.c p AT e (2.32) 


Ou AT e : ecart entre la temperature initiale de l’eau et la temperature a F entree de la 


batterie humide 
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Figure 2.12: Methodologie de calcul des besoins energetiques des elements principaux 

d’une machine frigorifique a adsorption 

2.3.7. Estimation des besoins journaliers 

L’ estimation des besoins energetiques quotidiens d'une machine frigorifique est 
resumee dans le tableau 2.7. 


Tableau 2.7: Besoins energetiques journaliers 


Puissance de la batterie Qb 

3.6 kW 

Duree de fonctionnement 

3 h 

Charge frigorifique journaliere Qfj 

38880 kJ 

Masse d’eau glacee me ((T s -T ee )=5°C) 

1860 kg 

Quantite de froid foumie par la machine 
frigorifique (evaporateur) Qu (AT=9°C) 

69984 kJ (19.44 kWh) 
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2.3.8. Dimensionnement des differents composants de la machine frigorifique a 
simple effet 

2.3.9. Methodologie de calcul 

Les donnees qui nous permettent le calcul de bilan et les differentes surfaces 
d’echange. Les donnees qui nous permettent le calcul du bilan et les differentes surfaces 
d’echange sont les suivantes : 

1 . masse du fluide frigorigene evapore 

2. masse de l’adsorbant 

3. surface du capteur 

4. chaleur specifique massique de l’adsorbant 

5. chaleur specifique massique du fluide frigorigene liquide 

6. chaleur specifique massique du fluide frigorigene gazeux 

7. chaleur specifique massique des metaux constituants le generateur 

8. chaleur latente de desorption 

9. chaleur latente de vaporisation du fluide frigorigene 


1. Calcul de la masse du fluide frigorigene evapore 


m ff ~ 


Q, 


L 


(2.33) 


Avec Qe la charge frigorifique globale journaliere (kJ/jour) et L v la chaleur latente 
d’ evaporation a T ivap (kJ/kg). 


2. Calcul de la masse de l’adsorbant 


m 


m ads — 


ff 


Am 


(2.34) 


Avec mjf : masse du fluide evapore (kg), Am la quantite de fluide desorbe par unite de 
masse de l’adsorbant (kg de fluide/kg d’adsorbam)- 

Cette quantite est obtenue a partir du reseau d’isosteres sur lequel on trace un cycle 
caracterise par quatre temperatures (T m i s , T reg , T evap et T com [ ou a partir de F equation de 
DIB ININ caracterisant le couple (adsorbant/adsorbat). 


3. Calcul de la surface de captation du generateur 

L’ensoleillement global pour la region consideree Ig et le COP s et de l’ordre delO% 
par m 2 done l’ensoleillement regu est : (/qx 0.1) 

Ceci correspond a une surface de (m 2 ): 
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Q F 

c — (2 

captation j q | ' — ' / 

Le bilan thermique est ecrit de la maniere suivante : 

Qg + Qev = Qcd + Qa ( 2 . 36 ) 

Qg : quantite de chaleur re§ue par l’adsorbant 
Qev '■ energie frigorifique 

Qcd '■ quantite de chaleur cedee par le condenseur 
Qa : quantite de chaleur cedee par l’adsorbant pendant 1 ’ adsorption 
Ou : 


Tddes Treg 

1 ) Qg = \Yj m i c P i dT+ \H m i c P i dT + m Amm ads ( 2 - 37 ) 

Tads Tddes 

2) Q ev = Lv.(T evap ).A m.m ads - Am.c p (T cond - T evap ).m ads (2.38) 

3) Q cd = Lv.(T cond )Am.m ads + Am.c p (T reg - T cond ).m ads (2.39) 

Tads (initial) Tads (fin) Tads (fin) 

4 ) Qa = \Yi m i c Pi dT+ lH m i c Pi dT+ \ MiAmm a ds ( 2 - 4 °) 

Treg Tads (initial) Tads (initial) 

Z">. C ,. rfr = Z& = +®2 +®3 (2-41) 


Avec SQi : chaleur sensible du metal ; 6Q2 : chaleur sensible de l’adsorbant et SQs : 
chaleur sensible du fluide. 


4. Calcul de l’epaisseur de l’adsorbant 


Connaissant la quantite de l’adsorbant et la masse volumique de celui-ci on peut 
calculer l’epaisseur de l’adsorbant dans le capteur : 


S cap = 2 -(5 


capatation 


) 


( 2 . 42 ) 


(Face avant et face arriere) 


* ads 


m 


ads 


P 


ads 


( 2 . 43 ) 


= 


V 


ads 


ads 


cap 


( 2 . 44 ) 


5. Calcul de la surface d’echange de l’evaporateur 

L’evaporateur sera constitue de tubes (cylindres horizontaux). L’echange avec l’eau 
s’effectue par convection naturelle : 
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Q,=K.S, mp AT 


(2.45) 


Q E : Charge frigorifique globale par unite de temps 
Sevap : Surface d’echange de l’evaporateur 

AT : Ecart de temperature entre celle du fluide frigorigene et de l’eau glacee. 


K s’exprime par la relation suivante: 


hi. di 


Done : 


1 


1 de 1 de , 

+ — + .In 


K hi. di he 2.X 


f de') 
\di j 


T 1 de de . 

Les termes et .In 


2.X 


l/ de\ 1 

— sont negligeables devant — 

di J he 

1 1 


K he 


(2.46) 


(2.47) 


Le coefficient he est determine par un calcul des nombres de Prandtl, Graschoff et 


Nusselt 


Nu = 


hc.D 


(2.48) 


he : coefficient d’echange convectif ; D : diametre du tube ; A : conductivite thermique de 
l’air ou de l’eau. Macadams propose F equation suivante (Convection naturelle laminaire 
ou turbulente autour d’un cylindre) pour le calcul de Nusselt pour les tuyaux horizontaux: 


Nu d = c.(Gr. Pr) 


n 


(2.49) 


Ou : 


Pr = 


c p -h 

~T~ 


(2.50) 


c p : chaleur specifique de l’eau ou de Fair 
r) : viscosite dynamique de l’eau ou de Fair 


p 2 .g.B.(T-TJ.D 


Gr = 


M 


2 


(2.51) 


Avec L est le diametre du tube. Dans notre cas D=de 


6. Calcul de la surface d’echange et la puissance echangee par le condenseur 

Le calcul se fait par F equation suivante: 


Q f = he.S F .AT 


(2.52) 


5, = 


Q 


he AT 
-38- 


(2.53) 
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Avec : Q E est la charge globale rapportee a F unite du temps (a la duree) de 


fonctionnement. On precede a un calcul similaire a celui utilise pour determiner la surface 
d’echange de l’evaporateur. Le condenseur doit evacuer la chaleur absorbee par la surface 
de captation et au moins sur une duree egale a celle de F ensoleillement 



1 COP S 

± solaire captation 


T 


(2.54) 


2.3.10. Application 

1. Calcul de la surface d’echange de l’evaporateur (cylindre horizontal) 

La temperature de la paroi estimee a 0°C est T evav = -2°C et d e = 40x1 O' 3 m, si on 
considere que la duree de F adsorption est egale a celle de la desorption (6h) done la 
puissance de l’evaporateur est : 




T 


adsorption 


(2.55) 


Done: Q E 


69984 
et 

6.3600 


Q e = 3.24 kW 


Tableau 2.8: Calcul de la surface d’echange de l’evaporateur 


Parametres de l’eau calcules a une temperature 
entre celle de la paroi et de l’eau glacee(r ra =5°C) 

Nombres adimensionnels avec leurs 
valeurs correspondantes 

c p =4.203 kJ/kg.K 

|0,=1.520 xlO' 3 kg/m.s 

X=0.5723 W/m.K 

C „ ./I 

Pr = ' =11.6 

/l 

p= 1000.0 kg/m 3 

g=9.81 m/s 2 

5=0.052xl0" 3 1/K 

T-T«=6 K 

Gr= P -8-B.(T-TJ.d e _ 84?84 48 

10 3 <Gr. Pr =946195 < 10 9 
convection laminaire n=l/4 
C=0.53 

Nu D =C.{Gr.Pr) n = 16.53 

Coefficient convectif, h c 
W/m 2 .K 

Surface d’echange, S e m 2 

Longueur du tube L m 

hc= X .Nu =263 50 
d e 

• 

Qr 

S F = ^ E =2.049 
hcAT 

L= Se = 16.3 
7t.d e 
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2. Calcul de la surface de captation et estimation de la puissance echangee par le 
condenseur (cylindre horizontal) 

Si on considere que la duree de captation ret que l'energie recue par une surface 
horizontale est de Ig et le COP s est de 10% : le condenseur doit done echanger la puissance 
suivante: 


Tableau 2.9: Estimation de la puissance echangee par le condenseur 


Le rayonnement global regn sur une surface horizontale Ig en W/m 2 

23174 

COP so laire CI1% 

10 

Duree de l’ensoleillement x en h 

6 

Surface de captation en m 2 

30.2 

Puissance du condenseur en kW 

3.24 


Tableau 2.10: Calcul de la surface de captation a T am =34°C, T c =49°C (temperature de la 

paroi exterieure du condenseur estimee a 47°C) et d e = 40x1 O' 3 m. 


Parametres de Pair calcules a une temperature 
entre celle de la paroi et de Pair ambiant 
(T m =40°C) 

Nombres adimensionnels avec leurs valeurs 
correspondantes 

c p = 1.0096 kJ/kg.°K 

(1=1.91 144x1 O' 5 kg/m.s 

X=0.02709 W/m.°K 

p r = Cp '^ = 0.7123 
A 

p=l. 11259 m 3 /kg 

g=9.81 m/s 2 

B= 1/Tm, 5=0.00319 1/K 

T-Too= 1 3 K 

Gr = P -S- B -( r ~ T J- d e =27652939.14 

A 

10 3 < Gr.Pr =19697188.45 < 10 9 
Convection laminaire, n=l/4 ; c=0.53 

Nu D =c.{Gr.Prf = 35.308 

Coefficient convectif, h c 
(W/m 2 .K) 

Surface d’echange S e 
(m 2 ) 

Longueur du tube L 
(m) 

hc=^ Nu =23.91 
d e 

• 

S = Qc = 10.42 
e hcAT 

L= Se = 82.9 

7t.d e 
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Tableau 2.11: Calcul de la surface de captation a T am =34°C, T c =44°C (temperature de la 


paroi exterieure du condenseur estimee a 42°C) et de =40x1 O' 3 m 


Parametres de Pair calcules a une temperature entre 
celle de la paroi et de Pair ambiant (T m =38°C) 

Nombres adimensionnels avec 
leurs valeurs correspondantes 

c p = 1.007 kJ/kg.K 

|i= 1.9x1 O' 5 kg/m.s 

X=0.0272 W/m.K 

C n .JU 

Pr = p =0.703 
/l 

p=1.13 m 3 /kg 

g=9.81 m/s 2 

B= 1/Tm, 5=0.00317 1/K 

T-T~=8 K 

„ p 2 .g. BfT-TJ.de 3 

Gr = 

M 

56318.09 

1 0 3 < Gr. Pr =39592< 1 0 9 
convection laminaire n=l/4 c=0.53 

Nu d = c.(Gr. Pr)" = 7.47 

Coefficient convectif h c 
(W/m 2 .K) 

Surface d’echange 
Se (m 2 ) 

Longueur du tube L (m) 

hc = ^' NU =5.08 
de 

• 

S= Qc =79.7 
hcAT 

L= Se = 634.4 
7t.de 


Tableau 2.12: Calcul de la surface de captation a T a m=Teau-ref=34°C, T c =44°C (temperature 
de la paroi exterieure du condenseur estimee a 42°C) et d e =40xl0~ 3 m. 


Parametres de l’eau calcules a une temperature entre 
celle de la paroi et de l’eau de refroidissement (T m =38°C) 

Nombres adimensionnels avec 
leurs valeurs correspondantes 

c p =4.179 kJ/kg.K 

p=0.68xl 0' 3 kg/m.s 

X=0.624 W/m.K 

C n-M 

Pr = p =4.6 
/l 

p=993 m 3 /kg 

g=9.81 m/s 2 

5=0.36xl0" 3 1/K 

T-T~=8 K 

^ p 2 .g.B.(T -TJ.de 3 

Gr - 2 

M 

3855867.8 

1 0 3 < Gr. Pr = 1 7736992< 1 0 9 
convection laminaire n=l/4 
c=0.53 

Nu d =c.{Gr. Pr)" =34.39 

Coefficient convectif h c 
(W/m 2 .K) 

Surface d’echange Se 
(m 2 ) 

Longueur du tube L 
(m) 

hc = ^' NU = 536.5 
de 

• 

S= Qc =0.75 
hcAT 

L= Se =6.01 

7t.de 


- 41 - 

































Chapitre 2: Conversion de Venergie solaire - refrigeration solaire 


Tableau 2.13: Bilan thermique et coefficients de performance 


Donnees permettant le calcul du bilan 

La masse du fluide frigorigene evapore niff 

28.1 kg 

la masse de l’adsorbant m a d S 

289.15 kg 

la surface du capteur S cap teur 
(30.2m 2 , (6.04x5x0.016) m) 

60.75 m 2 

la chaleur specifique massique de l’adsorbant cp a d S 

0.836 kJ/kg.K 

la chaleur specifique massique du fluide frigorigene liquide cpfa 

4.18 kJ/kg.K 

la chaleur specifique massique du fluide frigorigene gazeux cpff g 

1.91 kJ/kg.K 

la chaleur specifique massique des metaux constituants le generateur 
CPcapteur (cuivre) 

0.385 kJ/kg.K 

la chaleur latente de desorption Ldc S 

2490x1.2 kJ/kg 

la chaleur latente de vaporisation du fluide frigorigene L v 

2490 kJ/kg 


Avec : 

Qg + Qev = Qcd + Qa (2.56) 


Tableau 2.14: Calcul de coefficient de performance 


Quantites de chaleur 

Valeurs 

Qg 

136327 kJ 

Qev 

64918 kJ 

Qg+Qev 

201245 kJ 

Qcd 

75658 kJ 

Qa 

138533 kJ 

Qcd+Qa 

214191 kJ 

Coefficient de performance 

Formule et valeur 

COP (coefficient de performance 

Q p 

COP = £ = 0.47 

thermodynamique) 

Qg 

COP S oi a ire (coefficient de performance 
solaire) 

Qf 

COP solaire ~ j „ 

cap 

= 9.27% 
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2.3.11. Dimensionnement des equipements solaires d’une machine frigorifique 
fonctionnant avec le couple charbon actif/methanol 
1. Calcul de la surface de captation du capteur generateur 

Le capteur solaire est aussi le generateur des gaz du fluide frigorigene. La surface de 
captation est calculee par la formule suivante : 


QfJ 


I G .COP solaire 


(2.57) 


Le capteur generateur comprend une surface avant et une surface arriere. La surface 
totale de capteur est determinee a partir de la connaissance de la surface de captation 
(determinee a partir de la charge frigorifique) ainsi que la profondeur necessaire pour 
contenir une epaisseur donnee de l’adsorbant calculee par les formules suivantes : 

Connaissant la quantite de l’adsorbant et la masse volumique de celui-ci on peut 
calculer L epaisseur de l’adsorbant dans le capteur : 





(2.58) 




(2.59) 




(2.60) 


Avec P M( „ . : Perimetre de la surface de captation. 

' ^ capatation ' 

La profondeur du capteur a ete augmentee de 40% pour laisser un espace libre pour 
que la vapeur du fluide frigorigene puisse se deplacer librement a l’interieur du capteur- 
generateur. La masse de l’adsorbant est calculee par 1’ equation suivante: 



Am 


(2.61) 


Ou niff : masse du fluide evapore (kg) ; Am : la quantite de fluide desorbee par unite 
de masse de l’adsorbant (kg d efluide/kg d’adsorbant). 

Cette quantite est obtenue a partir du reseau d’isosteres sur lequel on trace un cycle 
caracterise par quatre temperatures (7%/ v , T reg , T evap et T con f) ou a partir de 1’ equation de 
DIB ININ caracterisant le couple (adsorbant/adsorbat). 

Si on considere que la duree de captation x et que l'energie re 5 ue par une surface 
horizontale est de Ig et le COP s est de 10% : la surface de captation est la suivante: 
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Tableau 2.15: Calcul de la surface de captation du capteur generateur 


Le rayonnement global refu sur une surface 
horizontale Ig en kJ/m 2 

23174 

COP so laire % 

10 

Duree de l’ensoleillement x en h 

6 

Surface de captation en m 2 

30.2 


2. Calcul de la masse de l’adsorbat 


L’ equation qui caracterise le couple (AC 3 5/methanol) est la suivante : 


In 




f Ps^ 

- n 

= -D 

T. In 



Wj 



{ P J 



Avec q (g/kg); p. (kg/1) ; P: (bar) ; Wo=0.425 1/kg; D=5.02x 1 0‘ 7 ' n=2. 1 5. 


(2.62) 


La masse desorbee Am=q est calcule entre la temperature ambiante (pression initiale 
de vaporisation a -2°C) et la temperature de regeneration sous la pression de condensation 
(49°C). 


Tableau 2.16: Calcul de la masse de l’adsorbant 


Conditions au debut de la desorption 

P (-2°C)) 

0.025 bar 

Tam 

34 °C 

Ps (T am ) 

0.233 bar 

p(T am ,P(-2°C)) 

665 g/1 

mi (masse initiale du methanol) 

150.75 g/kg 

Conditions a la fin de la desorption 

P (49°C) 

0.492 °C 

T 

1 rege 

90 °C 

Ps (Trege) 

2.558 bars 

p(T re ge,P(49°C)) 

614 g/1 

m 2 (masse finale du methanol) 

163.17 g/kg 

Am (masse desorbee) 

-12.42 g/kg 

Nouvelles conditions 1 

P(49°C) 

0.492°C 

T 

1 rege 

100 °C 

Ps(T r ege) 

3.572 bars 

p(T re2 e,P(49°C)) 

603 g/1 

m 2 (masse finale du methanol) 

121.23 g/kg 

Am (masse desorbee) 

29.52 g/kg 

Nouvelles conditions 2 

P (44°C) 

0.387 °C 

T 

1 rege 

100 °c 

Ps (trege) 

3.572 bars 

p(W,P (44°C)) 

603 g/1 

m 2 (masse finale du methanol) 

98.4 g/kg 

Am (masse desorbee) 

52.35 g/kg 
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69984 

1160 


(2.63) 


Done m„ 

.ft 


69984 

1160 


= 60.33 kg/ jour 


Avec Qe : charge frigorifique globale joumaliere (kJ/jour) et L v : chaleur latente 
d’ evaporation a T evap (kJ/kg). 


3. Calcul de la masse de l’adsorbant 

La masse de charbon actif necessaire est calculee par 1’ equation suivante: 


m 


ads 


Mrr 

JJ_ 

Am 


Avec m ads = 


60.33 

52.35 


60.33 


m ads ~ 


(52.35)10 


— Done m . = 1152.4 kg 


(2.64) 


2.3.12. Dimensionnement des equipements solaires d’une machine frigorifique 
fonctionnant avec le couple zeolithe/eau 

La masse du methanol ainsi que du charbon actif est tres importante done il y’aura un 
probleme de volume et du cout tres eleve de 1’ installation ce qui nous oblige a choisir un 
autre couple. A cet effet nous avons choisi l’eau/zeolite. L’eau a une chaleur latente de 
vaporisation de l’ordre de 2490kJ/kg.K (1160 kJ/kg.K) ce que signifie que la masse qui va 
etre utilisee pour le refroidissement de l’eau sera deux fois moins du methanol. Les 
donnees de L equation qui caracterise le couple (zeolite/eau) sont: Wo=0.263 l/kg ; 
D= 1.80x1 O' 7 et n=2.0. 


Tableau 2.17: Calcul de la masse de l’adsorbant 


Conditions au debut de la desorption 

P(1°C) 

0.0066 bar 

Tam 

34 °C 

Ps (T am ) 

0.0532 bar 

p(Tam,P(l°C)) 

994 g/1 

mi (masse initiale du methanol) 

242.78 g/kg 

Conditions a la fin de la desorption 

P (44°C) 

o 

• 

o 

h— ^ 
h— * 

o 

n 

T 

1 rege 

150 °C 

Ps (Trege) 

4.758 bars 

p(T re ge, P(44°C)) 

916.6 g/1 

m 2 (masse finale du methanol) 

145.58 g/kg 

Am (masse desorbee) 

97.2 g/kg 
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La masse de l’adsorbat (eau) = 28.1 kg et la masse de l’adsorbant (zeolithe) = 289.15 kg. 
Pour une masse volumique apparente, de la zeolithe, egale a 600 kg/m 3 le V Z e=0.482 in et 
h Z e=0.0159 m ( h Z e=1.6 cm). 

2.3.13. Conclusion 

Ce travail a permis de dimensionner les differents elements d’une machine 
frigorifique solaire a adsorption. Le couple adsorbant/adsorbat est l’eau/zeolithe, a travers 
cette partie on conclut que. 

La quantity d’energie absorbee par le systeme est presque celle cedee au niveau du 
condenseur et pendant la phase d’ adsorption. La difference apparente est due a 
l’incertitude dans la determination de la temperature au debut de la desorption d’une part 
d’ autre part a la temperature du debut de 1’ adsorption qui est consideree dans notre cas 
egale a la temperature ambiante etc. 

Nous avons considere dans cette etude une duree de fonctionnement de trois heures 
avec un regime de consommation maximal. II est done indispensable de recourir a 
1’ experience (utiliser des petites quantites du couple (eau/zeolithe)) pour determiner avec 
precision la duree exacte de 1’ adsorption et de la desorption, ce qui nous permettra ensuite 
d’etudier la possibility de resoudre le probleme de 1’ intermittence du cycle par l’utilisation 
de plusieurs capteurs generateurs. Cette possibility peut aboutir a une diminution de la 
masse d’adsorbant et de l’adsorbat utilises. 

Le calcul a ete fait en programmant toutes les equations citees precedemment sur un 
programme en Fortran. 

Les conclusions quand on peut tirer de cette partie concernant la production de froid 
par adsorption solaire sont : 

• II ressort done de cette etude que le methanol est un excellent adsorbat pour les 
applications envisagees, et le charbon actif un tres bon adsorbant. Ces deux 
constituants forment actuellement le couple ideal, selon le critere etudie. 

• L’energie solaire constitue une reserve d'energie inepuisable a l'echelle de la vie 
de l'humanite. Cette energie n'est pas exploitee par les pays qui en disposent, en 
majorite des pays en voie de developpement. Car le froid solaire est surtout 
interessant dans les regions ou relectricite n'est pas disponible, et dans les regions 
chaudes et isolees. 

• Les unites de refrigeration a adsorption qui utilisent les energies renouvelables, 
pourraient contribuer au developpement durable et preserver l’environnement. 
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• Etant donne la concordance etroite de la demande de refroidissement avec la 
disponibilite du rayonnement solaire, le refroidissement solaire permet de reduire 
judicieusement la consommation electrique en ete. 

• La modelisation des cycles a adsorption solide necessite une bonne connaissance 
du phenomene d’ adsorption. Done, il est essentiel de maitriser des informations 
concernant l’equilibre, ce dernier necessite un choix convenable de modele 
d’ adsorption (isotherme d’ adsorption) qui doit decrire correctement les proprietes 
d’ adsorption du couple utilise. 
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3. Les reseaux de neurones artificiels (RNA) 

3.1. Introduction 

Les reseaux de neurones artificiels (RNA) sont des systemes paralleles et distribues de 
traitement de 1’ information inspires par le fonctionnement du cerveau humain. Le present chapitre 
decrit le principe de fonctionnement des RNA, leurs differentes architectures, les families 
d’algorithme d’apprentissage disponibles ainsi qu’un resume des principaux avantages des reseaux 
de neurones en tant que systeme de traitement de 1’ information. 

II ressemble au cerveau humain sur deux aspects: 

• La connaissance est acquise par le reseau de son environnement a travers un processus 
d’apprentissage. 

• Les forces de liaison entre neurones (poids synaptiques) sont utilisees pour stocker la 
connaissance acquise. 

Une presentation succincte des notions fondamentales des RNA et du vocabulaire necessaire a 
la comprehension de la partie modelisation par les RNA a ete faite. Ensuite, nous examinons 
1’ application des RNA de type feedforward, en particulier les perceptrons multicouches ( MLP , 
Multi Layer Perceptron) comme approximateurs universels de fonctions, (Joduin. 1994). 

3.2. Notions fondamentales sur les reseaux de neurones artificiels 
3.2.1. Analogie entre le neurone biologique et le neurone artificiel 

L'idee des reseaux neuraux artificiels a ete inspiree de la maniere par laquelle les neurones 
biologiques traitent les informations. Ce concept est utilise pour mettre en oeuvre des logiciels de 
simulations des processus massivement paralleles qui impliquent des elements de traitement 
connectes dans une architecture de reseau. L'apprentissage dans le cerveau humain se fait a travers 
un reseau des neurones qui sont connectes par l’axone, les synapses et des dendrites, (Bouzidi. 
2006 ; Rutkowski. 2004) (figure 3.1). 



Figure 3.1 : Neurone biologique (Bouzidi. 2006 ; Rutkowski. 2004). 
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Le neurone biologique est une cellule composee d’un corps cellulaire et d’un noyau. Le corps 
cellulaire se ramifie pour former les dendrites. C’est par les dendrites que 1’ information est 
acheminee de l’exterieur vers le soma, corps du neurone. L’information traitee par le neurone 
chemine ensuite le long de l’axone pour etre transmise aux autres neurones. La jonction entre deux 
neurones est appelee la synapse. Les dendrites, recueillent, via les synapses, les signaux 
electrochimiques ou les potentiels d’ action venant d’ autres cellules ou de l’exterieur et les 
acheminent vers le corps cellulaire ou le soma. Le soma, en outre son role concernant le 
metabolisme de la cellule, collecte et concentre les informations et en fait une sommation. Si le 
potentiel somatique depasse un certain seuil, il y a emission d’un potentiel d’ action, qui correspond 
a une oscillation electrique tres breve. L’ axone est un prolongement cellulaire charge de faire passer 
1’ information nerveuse a d’ autres cellules a travers les synapses grace a son arborisation terminale. 
Les synapses se trouvent au point de contact entre neurones, elles convertissent les impulsions 
nerveuses en signaux chimiques c’est-a-dire en molecules de neurotransmetteur. En effet, l’arrivee 
d’un potentiel d’action a l’extremite d’un axone entraine le passage d’une quantite de 
neurotransmetteurs dans l’intervalle synaptique. Le nombre de synapses par neurone varie 
considerablement de plusieurs centaines a une dizaine de milliers, (Si-moussa. 2010 ; Mahmoudi. 
2008). 

Le neurone artificiel (figure 3.2) refoit les entrees qui sont analogues aux impulsions 
electrochimiques que les dendrites de neurones biologiques re§oivent d'autres neurones. II est done 
un operateur mathematique tres simple qui traduit la representation mathematique d’un neurone 
biologique. II possede des entrees qui peuvent etre les sorties d’autres neurones, ou des entrees de 
signaux exterieurs, et une sortie. La valeur de la sortie resulte du calcul de la somme des entrees, 
ponderees par des coefficients (poids synaptiques) et du calcul d’une fonction non lineaire bomee 
de cette somme ponderee. Cette fonction (fonction d’ activation) est le plus souvent une fonction 
sigmoide ou un echelon, (Bouzidi. 2006 ; Abraham. 2000). 



Signaux ou 
variables 
d’entree 


Poids 

synaptiques 


Figure 3.2 : Neurone artificiel, (Bouzidi. 2006 ; Mouhamdi. 2008). 
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La figure 3.3 resume la correspondance entre le neurone artificiel et le neurone biologique. 
On remarque dans cette description schematique, que chaque composant ou constituant du neurone 
artificiel trouve son inspiration biologique. 



Synapses 


Foids 
synaptiques 


Dendrites 


Fonction 
d’ activation 


Noyau 


Corps cellulaire 
ou Soma 

^ Axone 


Signaux ou 
variables 
d' entree 


Figure 3.3 : Comparaison entre un neurone biologique et artificiel, (Si-moussa. 2010; 
Bouzidi. 2006). 


3.2.2. Concept du reseau de neurones de type feedforward 

Les neurones artificiels sont arranges en couches (figure 3.4), ou la couche d’entree refoit des 
signaux d’entree (n, ) du monde reel et chaque couche successive recoit les sorties ponderees ( w ij u i ) 

de la couche precedente comme son entree resultant ainsi en un RNA de type feedforward (cas de 
perceptrons multicouches), dans lequel chaque sortie d’une couche est alimentee en avant a la 
couche qui la succede ou elle est traitee. Les sorties de la demiere couche constituent les sorties au 
monde reel (Rafinia. 2014). 

n 

V/ 

◄ ► ◄ ► < ► 

Couche d ’entree Couche cachee Couche de sortie axec 

axec n neurones axec m neurones l neurones ei l sorties 

ei n entrees. 

Figure 3.4: Reseau de neurones artificiels de type feedforward. 
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Dans un tel RNA de type feedforward un neurone, dans une couche cachee ou une couche de 
sortie, a deux taches, (Si-moussa. 2010; Rafinia. 2014): 

II additionne le biais aux entrees ponderes de plusieurs connexions et applique ensuite une 
fonction d’ activation (dite aussi fonction de transfert ou encore fonction de seuillage) a cette somme 
comme donnee par l’equation 3.1 (pour le neurone j de la couche cachee de la figure 3.4), (Hornik. 
1990) : 



f n ^ 

f h +I, m 

V i = 1 J 


j = l,2,...,m 


(3.1) 


• II propage la valeur resultante par les connexions partantes aux neurones de la couche 
suivante ou elle subit le meme processus tel que donne par 1’ equation 3.2 (par exemple les 
sorties de la couche cachee alimentee au neurone k de la couche de sortie donnent la sortie 
v k ), (Hornik. 1990) : 



/< 


o 


f m 

z 

W=1 


3 


KjZj+Kk 


j 


k = 1,2,...,/ 


(3.2) 


En combinant des equations 3.1 et 3.2 on obtient la relation entre la sortie v, et les entrees 


u t du RN: 


1 

^ m 

f n ' ^ 

\ 

II 

YAjh 

H w ji u i +b hj 

+ b ok 

1 

\ j =i 

V i= 1 ) 

J 


k = 1,2,...,/ 


(3.3) 


Le nombre de neurones dans la couche d’ entree et la couche de sortie est fixe par le nombre 
de variables independantes et dependantes respectivement. Le modelisateur definit le nombre de 
couches cachees et le nombre de neurones dans chaque couche cachee. Le developpement du 
modele est realise par un processus d’ apprentissage dans lequel un ensemble de donnees 
experimentales des variables independantes est presente a la couche d’ entree du reseau. Les sorties 
de la couche de sortie comprennent une prediction des variables dependantes. Le reseau apprend les 
relations entre les variables independantes et dependantes par la comparaison iterative des sorties 
calculees et des sorties experimentales et l'ajustement posterieur de la matrice de poids et le vecteur 
des biais de chaque couche par un algorithme d’ apprentissage de retro propagation du gradient de 
l’erreur. Ainsi le reseau developpe un modele neuronal capable de predire avec precision acceptable 
les variables de sortie comprises dans le domaine des donnees d’ apprentissage. Par consequent, 
l'objectif de la modelisation neuronale est de reduire au minimum les erreurs de prediction de 
L ensemble des donnees de validation presentees au reseau a la fin de l'etape d’ apprentissage, 
(Hornik. 1990). 
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3.2.3. Fonctions d’ activation 

La fonction d’ activation calcule la valeur de l’etat du neurone qui est transmise aux neurones 
avals. La fonction d’ activation sert aussi a introduire une non-linearite dans le fonctionnement du 
neurone. A la difference des neurones biologiques dont l’etat est binaire, la plupart des fonctions de 
transfert sont continues, offrant une infinite de valeurs possibles comprises dans l’intervalle [0, +1] 
ou [-1, +1]. II existe de nombreuses formes possibles pour la fonction de transfert. Le tableau 3.1 
resume quelques fonctions de transfert utilisees dans le programme de calcul techniques MatLab® 
(Hagan et al. 1996 ; Rafinia. 20014). 

En outre, du fait que la continuite et la differenciabilite des fonctions d’ activation a tous les 
points sont les caracteristiques requises par les algorithmes d'optimisation actuels, les fonctions 
d'activation typiques qui accomplissent ces exigences sont la fonction lineaire, la fonction tangente 
hyperbolique et la fonction sigmoide logarithmique. Selon (Samarasinghe. 2006), la fonction 
lineaire est plus appropriee pour les neurones de sorties alors que les fonctions sigmoides sont plus 
appropriees pour les neurones des couches cachees. 


Tableau 3.1: Fonctions de transfert usuelles, a=f(n), (Samarasinghe. 2006) 


Nom de la fonction 

Relation entree / 
sortie 

Graphe de la 
fonction 

Fonction MATLAB® 
correspondante 

Lineaire 

a = n 

a 

i 

/ 

-i 

purelin 

Logistique 

(Sigmoide 

logarithmique) 

i 

a = 

l + e~ n 

i 

i 

a 

•N+l 

■ 

0 

-1 

logsig 

Tangente hyperbolique 

e + n _ e n 

a 

+n . -n 

e + e 

( 

J 

a 

M 

n 

0 

-I 

tansig 


3.2.4. Architecture des reseaux de neurones artificiels 

Nous avons introduit precedemment le concept de RNA de type feedforward, il existe 
cependant d’autres types de RNA. Une description succincte, basee sur le guide d’utilisateur de la 
boite a outils des RN de MatLab® (Demuth et al. 2000), des principaux types de RNA est donne 
dans ce qui suit, (Demuth et al. 2000): 

• Les reseaux de neurones a fonctions radiales de base 

• Les reseaux competitifs ou cartes auto organisatrices et les reseaux recurrents 
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3.2.5. Selection d’un reseau de neurones artificiel 

II est clair que l'application reussie des RNA pour la modelisation des problemes d'ingenierie 
est fortement affectee par quatre facteurs majeurs, (Pham. 2006a): 

• Type de reseau, 

• Structure de reseau (nombre de couches cachees, nombre de neurones par couche cachee), 

• Fonctions d’ activation, 

• L’algorithme d’apprentissage. 

3.2.5. 1. Selection du type de reseau 

Parmi les reseaux de neurones de type forward qui possedent la propriete d’ approximateur 
universel de fonction, on trouve les MLP. Ils presentent ces caracteristiques, (Pham et al. 1995) : 

• Un MLP peut avoir une ou plusieurs couches cachees. 

• Pour un MLP, les neurones des couches cachees ou de la couche de sortie ont le meme 
modele. 

• Les couches cachees et la couche de sortie d’un MLP, sont toutes non lineaires. Cependant, en 
cas de regression non lineaire la couche de sortie est de preference lineaire. 

3.2.5.2. Selection de la structure du reseau 

II est bien etabli que la variation du nombre des neurones de(s) couche(s) cachee(s) a un effet 
significatif sur la capacite de prediction du reseau. Le nombre adequat de neurones caches depend 
d'une facon complexe de, (Behrang. 2010) : 

• nombres de neurones de la couche d’ entrees et de sorties ; 

• la taille ou la dimension de l’echantillon d’apprentissage ; 

• la quantite de bruit dans les cibles (sorties du reseau) ; 

• la complexity de la fonction ou de la classification a faire apprendre au reseau ; 

• l'architecture du reseau ; 

• type de fonction d'activation des neurones caches ; 

• l’algorithme d’apprentissage ; 

• la regularisation 

Dans la plupart des situations, il n’y a aucune facon de determiner le nombre adequat de 
neurones caches sans faire l’apprentissage de plusieurs reseaux et d’evaluer ensuite l'erreur de 
generalisation de chaque reseau. Tres peu de neurones caches donnent des erreurs d’apprentissage 
et de generalisation trop elevees en raison du sous apprentissage. Par contre un nombre est trop 
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eleve de neurones caches, donne une erreur d’apprentissage tres faible mais une erreur de 
generalisation top elevee en raison du sur apprentissage, (Behrang. 2010). 

Plusieurs regies heuristiques ont ete proposees pour fixer le nombre de neurones caches. Selon 
(Curry et al. 2006), a present il n’y a aucune methode qui permet de determiner le nombre approprie 
de neurones caches afin de realiser une performance donnee. Done, il est necessaire de se servir 
d'une methodologie specifique. Pour (Plumb et al. 2005) la fa£on la plus commune d'optimiser la 
performance d’un RNA est en variant le nombre de neurones dans la couche cachee et choisissant 
l'architecture avec la meilleure capacite de generalisation. 

3.2.5.3. Selection de la fonction d’activation 

Les fonctions d'activation pour les unites cachees sont necessaires pour introduire la non 
linearite dans le reseau. Pour 1’ apprentissage par retro propagation, la fonction d'activation doit etre 
differentiable et il est preferable qu’elle soit bomee; les fonctions sigmoides comme logsig et tansig 
sont les choix les plus habituels. La fonction tansig qui produit des valeurs positives et negatives a 
tendance a donner un apprentissage rapide que des fonctions qui produisent seulement des valeurs 
positives comme la fonction logsig, (Behrang. 2010) 

3.2.5.4. Selection de type d’apprentissage : Apprentissage supervise 

Dans 1’ apprentissage supervise, un superviseur fournit une valeur ou un vecteur de sortie 
(appele cible ou sortie desiree) que le reseau de neurones doit associer au vecteur d’ entree, done 
L apprentissage supervise implique L existence d’un "professeur" qui a pour role d’ev aluer le succes 
(ou l’echec) du reseau quand il lui est presente un stimule connu (Abraham. 2000 ; Mahmoudi. 
2008). 



La regie de Widrow-Hoff "regie de delta" est une methode de minimisation de 1’ erreur entre 
la sortie reelle et la sortie desiree. Cette regie est generalement utilisee dans le cas d’apprentissage 
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supervise. Done les fonctions d’ activations sont des neurones lineaires. Elle converge vers la 
solution des moindres carres qui minimise la fonction d’erreur E. Son but est de faire evoluer le 
reseau vers le minimum de sa fonction d'erreur (erreur commise sur l'ensemble des exemples). Elle 
est utilisee dans le modele de l'adaline (ADAptive LINear Element). L'apprentissage est realise par 
iteration (les poids sont modifies apres chaque exemple presente), et on obtient le poids a l'instant 
t+1 par la formule (Abraham. 2000 ; Mahmoudi. 2008) : 

W(t + i) = W(t)+n.(T-0').E (3.4) 

Avec: E = (0. - L,.) 2 (3.4) 

^ i 

Avec W est le poids, T la sortie theorique et O la sortie reelle. E l'entree et n un coefficient 
d'apprentissage (entre 0 et 1 ) que l'on peut diminuer au cours de l’apprentissage. C'est en fait un cas 
particulier de l'algorithme de retro propagation du gradient, (Abraham. 2000 ; Mahmoudi. 2008). 

3.2.5.5. Selection de l’algorithme d’apprentissage 

L’apprentissage est la procedure par laquelle les valeurs des poids synaptiques du reseau sont 
determinees et la structure de reseau est etablie. Plusieurs algorithmes d’apprentissage utilisent la 
procedure de retro propagation et bien que chacun ait ses propres avantages, comme le taux de 
calcul et les exigences de calcul et de memoire de stockage, il n’y a aucun algorithme qui convient a 
tous les problemes. La performance de chaque algorithme depend du processus modelise, de 
l'echantillon d’apprentissage et le mode d’apprentissage. Cela explique que le choix des algorithmes, 
quoique d'une importance primordiale, depend etroitement de l'application concernee et demeure 
principalement experimental, (Hagan. 1994). 

L’algorithme de retro propagation du gradient comporte, done, deux phases (Hagan. 1994): 

• La premiere phase "forward" concerne le signal fonction qui se propage vers l’avant du 
reseau, d’un neurone a un autre, apres initialisation des poids synaptiques et des biais. 
Pendant cette phase, les poids synaptiques restent inchanges. Cette phase commence, done, a 
la premiere couche cachee par le calcul du signal fonction de tous les neurones qui la 
constitue et se termine a la couche de sortie par le calcul de la fonction d’erreur de tous les 
neurones qui la constitue aussi. 

• La deuxieme phase "backward" concerne le signal d’erreur qui se propage, dans le sens 
contraire, couche par couche, a partir de la couche de sortie. Pour chaque neurone, le 
gradient local de 1’ erreur est calcule pour permettre la correction des poids synaptiques, 
liaison par liaison, selon une regie connue sous le nom de la regie Delta. Dans le cas d’un 
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neurone de sortie, le gradient local est egal au signal d’erreur multiplie par la derivee de la 
fonction d’ activation appliquee a F entree nette de ce neurone. 

L'algorithme de retropropagation peut prendre deux formes ou deux modes. La premiere, dite 
retropropagation incrementale ou sequentielle, consiste a modifier les poids apres chaque 
presentation. Au contraire, la retropropagation par lots ( batch backpropagation ) n'ajuste les poids 
qu'au bout d'un certain nombre de presentations (une epoque), generalement egal au nombre 
d'exemples de l’echantillon d'apprentissage. Sous cette forme, on evite l'interference entre les 
modifications successives des poids : en calculant le gradient sur la somme des erreurs obtenues 
pour chaque exemple du lot, on effectue une sorte de moyenne. Cependant, si la base 
d'apprentissage est tres importante, le besoin en capacite de memoire de stockage et la redondance 
des informations rend la procedure par lot moins efficace que la procedure sequentielle, (Hagan. 
1994 ; Hassani. 2003). 

Puisque la retropropagation se resume a un probleme de recherche de minimum, son 
utilisation souleve au moins trois questions liees a F optimisation multidimensionnelle, (Henrique. 
2000 ; Plumb. 2002) : 

• Vers quel minimum la convergence s’effectue-t-elle ? 

V 

• A quelle vitesse l’algorithme converge-t-il ? 

• Quelles sont les capacites de memoire necessaires ? 

Les algorithmes de retro propagation implantes actuellement dans les logiciels de RNA 
varient en fonction de la methode d’ optimisation utilisee et ont pour but de porter la meilleure 
solution a ces trois questions, (Si-moussa. 2010). 

Un probleme important a tendance a apparaitre lors d’un apprentissage est celui du sur 
apprentissage. Si le reseau de neurone apprend par coeur, il donnera de mauvais resultats quand on 
lui presentera des donnees exclues du jeu d’apprentissage. Des methodes existent pour optimiser la 
phase d’apprentissage afin que le phenomene de sur ou sous apprentissage disparaisse, dont la 
technique de F arret precoce (earlystopping) et la technique de regularisation. La premiere technique 
consiste a diviser les donnees disponibles en trois lots distincts. Le premier lot sert a F apprentissage 
du RN, le second lot sert a tester le reseau et le troisieme lot sert a verifier si le reseau a une bonne 
capacite de generalisation, (Mellit. 2009 ; Cyril. 2014). 

Pour le deuxieme lot l’erreur (entre les predictions et les cibles) de test, calculee au cours de 
la phase d’apprentissage, doit normalement diminuer (la variance diminue). Mais quand le reseau 
commence a apprendre par cceur (le biais augmente, alors l’erreur de test recommence a croitre), 
F apprentissage est alors arrete avant la fin d’oii le nom de "apprentissage avec arret precoce" 
(stopped training) donne par (Bourquin et al. 1997) ou "apprentissage attenue" (attenuated training) 
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donne par (Plumb et al. 2002).L’apprentissage avec arret precoce est generalement applicable avec 
les algorithmes de retro propagation de Levenberg-Marquardt. La technique de la regularisation 
Bayesienne utilise une fonction de performance modifiee con§ue pour reduire au minimum le sur 
apprentissage par le lissage de la surface d'erreur de l’echantillon d’apprentissage. 

Le tableau 3.2 resume les algorithmes d’apprentissage par retro propagation implantes dans la 
boite a outils des RNA du programme MatLab®(Demuth et al. 2000). 


Tableau 3.2 : Resume des caracteristiques des algorithmes d’apprentissage de retro 
propagation implantes dans MatLab® (Demuth et al. 2000). 


Algorithme 

Caracteristiques 

Levenberg-Marquardt 

Levenberg-Marquardt ( trainlm ) 

C'est l'algorithme d’apprentissage le plus rapide pour les reseaux 
de taille moderee. 11 a la caracteristique de reduction de memoire 
pour l'utilisation quand l’echantillon d’apprentissage est grand. 

Regularisation 

Regularisation Bayesienne 

(trainbr) 

Mise a jour des poids et des biais selon l'optimisation de 
Levenberg-Marquardt. 11 reduit au minimum une combinaison 
d'erreurs quadratiques et des poids et determine ensuite la 
combinaison correcte afin de produire un reseau avec une 
capacite de generalisation meilleure. Ce processus est appele la 
regularisation Bayesienne. 11 reduit la difficulte de determiner 
l'architecture de reseau optimum. 


3.3. Conclusion 

Tout au long de ce chapitre, nous avons donne des informations sur les reseaux de neurones 
d’une fa£on generale tout en mettant 1’ accent sur les reseaux de neurones de type feedforward, les 
plus utilises dans les modelisations statiques du fait qu’ils sont doues de la propriete 
d’ approximation universelle. Nous avons presente un algorithme d’apprentissage tres utilise dans la 
communaute connexionniste: 1’ algorithme de retro propagation du gradient. Plusieurs travaux lui 
ont ete consacres et un certain nombre de modifications ont ete proposees pour accelerer le temps 
d’apprentissage, reduire les besoins en memoire de stockage et ameliorer la capacite de 
generalisation ainsi que les performances du reseau obtenu. 

Vers la fin, nous avons egalement presente la strategic implantee dans les programmes 
elabores pour la modelisation de la performance d’un refrigerateur solaire a adsorption. 
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Chapitre 4 

Partie experimentale 



Donnees de depart 


Dimensions exterieures et interieures de la 
chambre froide (longueur, largeur et hauteur) 
normalisees suivant ISO 1496-6 de 1988 



Bilan thermique d’une chambre froide 


Les charges thermiques autour d’une chambre 
froide se resument en deux composantes: les 
charges internes et externes 



Fabrication et installation des chambres 

froides 


Elies repondent a des normes de securite et 
d’ hygiene. Les normes en vigueur sont NF E 
35-400 pour 1’ installation frigorifique et NF 
C 15-100 pour l’installation electrique. Le 
respect de cette norme a une influence 
primordiale sur la qualite des produits. 



Calcul de l’epaisseur de la couche d’ isolation en mousse de 
polyurethane (stable en fonction de temps, moins chere et 
legere) et en limitant le coefficient de transfert thermique 
global U < 0.36 W/m 2 /K pour les chambres froides positives. 


1. Les charges externes sont constitutes essentiellement de : 

La charge thermique echangee par transmission a travers chaque 

paroi de la chambre froide, eq. 2.5-2.7 

La charge thermique par ouverture des portes (Q 0 =0). 

La charge thermique due au renouvellement d’air, eq. 2.8 a 2.12 

2. Les charges internes sont subdivisees en deux categories : 

^ Les charges thermiques independantes des denrees entreposees : 

- La charge thermique due aux personnes, eq. 2.14 

- La charge thermique due aux eclairages, eq. 2.13 

- La charge thermique due au materiel roulant (inexistante). 

- La charge thermique due aux machines di verses (inexistante). 
S Les charges thermiques dependantes des denrees entreposees : 

- La charge thermique due aux denrees ententes, eq. 2.15-2.17 

- La charge thermique due a la respiration des denrees 

(negligeable parce que on a utilise l’eau comme charge). 

- La charge thermique due aux moteurs des ventilateurs des 
evaporateurs (negligeable). 

- La charge thermique due a la resistance de degivrage 

(inexistante). 

^ Calcule la puissance frigorifique necessaire 

- Calcul de dimensionnement des differents elements de la 

chambre froide englobant : Compresseur, evaporateur, 

condenseur et detendeur, eq. 2.19-2.28 


Source d’energie electrique 


Utilisation d’un panneau photo voltaique pour 
la conversion de l’energie solaire en 
electricite 



Calcul de 1’ irradiation solaire regue sur un plan incline 
Calcul de la puissance crete necessaire et la meilleure 
disposition et inclinaison 
Choix des modules a installer (poly-cristallin) 


Protection de l’environnement 


Recherche des composes chimiques qui ne 
polluent pas pour les refrigerateurs a 
compression mecanique ou a adsorption 



Energie de fonctionnement 


Besoins energetiques necessaire au bon 
fonctionnement des differents organes du 
refrigerateur 



Le refrigerant R134a est un nouveau refrigerant 
ecologique. Son ODP (potentiel d'appauvrissement de la 
couche d'ozone) est de 0, aussi il ne cause aucun dommage 
a la couche d'ozone 

Le methanol peut facilement s'evaporer a des temperatures 
inferieures a 0 °C ce qui favorise la production de glace. 
L'electricite est conservee sous forme electrochimique afin 
de resservir a la demande quand le soleil ne brille plus. 


Maximiser la conversion d’energie solaire en electricite 
Minimiser la capacite de stockage et son autonomie 
Maximiser le coefficient de performance 


Conforts et cout 


Bruit, vibration et encombrant 
Le cout total du systeme en DA 
Schema de l’installation solaire 



Le niveau sonore correspond a la directive 2000/14/EC 
Le cout total en DA de systeme en tenant compte le cout de 
transport au cours de la phase de test, 
schema du systeme experimental pour un refrigerateur 
solaire PVAC, Figure 4.1.10 

Schema de la machine frigorifique solaire a adsorption, 

Figure 2.11. 


Figure 4.1.1a: Organigramme detailles de la Problematique du sujet 






















Contribution 1 


Realisation d'un refrigerateur 
solaire photovoltaique 
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Contribution 1: Realisation d’un refrigerateur solaire photovoltaique 


4.1.1. Introduction 

Les pics de consommation electrique nationale en Algerie se succedent a un rythme 
rapide. En ete 2012, la Sonelgaz (entreprise Algerienne de distribution d’electricite) a 
enregistre de nouveaux records de consommation electrique durant les periodes de canicule. 
Les consequences sont connues : delestages, coupures inopinees, . . . etc. Face a cette situation 
causee par le reseau conventionnel, une solution alternative pourrait etre utilisee par la 
Sonelgaz, c’est la rationalisation de l’utilisation de l’electricite et l’efficacite energetique qui 
reste une notion inconnue en Algerie et la diversification des ressources energetiques. Des 
estimations de l’ANPPRUE (Agence Nationale Pour la Promotion et la Rationalisation de 
l’Utilisation de l’Energie) ont montre que pres de 46% de la consommation d’electricite des 
menages est due a la refrigeration et la climatisation. Comme F Algerie est un pays riche en 
potentiel energetique renouvelable, elle planifie de s’orienter d’avantage vers les energies 
renouvelables pour varier ses ressources energetiques. L’utilisation de l’energie solaire pour 
la refrigeration et la climatisation est une application a promouvoir, car les besoins en froid 
coincident la plupart du temps avec la disponibilite du rayonnement solaire. Pour contribuer 
au developpement des systemes frigorifiques photo voltaiques en Algerie, un projet commun a 
ete conduit entre les deux laboratoires : le Laboratoire de Biomateriaux et Phenomenes de 


Transport (LBMPT- Universite de Medea) et le Laboratoire de Froid et Climatisation par 
Energies Renouvelables (LFCER- Unite de Developpement des Equipements Solaire de Bou- 
Ismail). Deux refrigerateurs a compression de vapeur de meme volume (160 litres) et une 
chambre de 1.5 m 3 seront dimensionnes pour etre alimente soit directement ou indirectement 
par energie solaire photovoltaique. En suite ces systemes de refrigeration seront testes sous 
deux differents climats: la cote mediterraneenne (Bou-Ismail-Tipaza) et le Sahara 
(Tamanrasset). Plusieurs modifications et ameliorations seront envisagees afin de fabriquer un 
prototype de refrigeration solaire photovoltaique avec une consommation minime et a 
moindre cout. Le premier refrigerateur (1) a courant alternatif sans modification du cycle 
frigorifique sera dimensionne pour etre alimente directement en energie solaire 
photovoltaique. Le deuxieme (2) refrigerateur sera alimente indirectement par energie solaire 
ou les equipements du cycle conventionnel ont ete remplaces par un kit englobant un 
compresseur/condenseur BD50F a courant continu et un evaporateur eutectique AC 15. En fin, 
la chambre (3) sera alimentee indirectement par energie solaire. 
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Realisation d’un refrigerateur solaire photovoltaique 


1. Dimensionnement d’une installation solaire. 

2. Estimation des besoins energetiques 

3. Estimation de l'ensoleillement 

4 . Estimation de la capacite de stockage 

5. Estimation du champ photovoltaique requis 

6 . Dimensionnement des cables electriques solaires 

7. Dimensionnement d'un regulateur de tension et de l’intensite 

8. Dimensionnement du convertisseur (Onduleur) 

9. Logiciel PVSyst 

10 . Caracteristiques de systeme de refrigeration 

11. Dispositif experimental 


1. Optimisation du dimensionnement de 1’ installation 
frigorifique solaire DC-AC sans modification dans 
le cycle frigorifique 


2 . Optimisation du dimensionnement d’un systeme 
frigorifique solaire direct DC-DC avec accumulation 
de froid 


Optimisation du dimensionnement d’une 
chambre frigorifique de 1.5 m solaire direct 
DC-AC. 


Optimisation du dimensionnement par PVSyst 

Methodologie et techniques de mesure 

Calcul de la surchauffe et de sous-refroidissement 

Bilan energetique du cycle frigorifique et coefficient de performance (COP) 
Diagramme enthalpique par le logiciel CoolPack 
Etude du comportement du compresseur 


1. Fabrication de la chambre frigorifique 

2. Dimensionnement et test de systeme frigorifique solaire photovoltaique a 
l’UDES 

3. Test sous les conditions climatiques d’ Ain-Guezzam (Tamanrasset) 


Dispositif experimental 

Kit BD50F de compression/condensation et evaporateur eutectique 
Optimisation du dimensionnement par PVSyst 
Resultats et discussion 

Comparaison entre le systeme de refrigeration DC-AC et le systeme DC-DC 
avec stockage du froid. 
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L'optimisation du dimensionnement de ces trois systemes photovoltai'ques a ete 
effectuee avec le logiciel de dimensionnement des systemes et des installations solaires 
photovoltai'ques PVSyst V5.59. Les differentes taches de la premiere contribution sont 
regroupees dans la figure 4.1.1. 

4. 1.1.1. Dimensionnement d’une installation solaire 

Pour une installation autonome avec stockage de l’energie produite par le generateur 
photovoltaique et sans appoint d'energie, le dimensionnement sert a calculer la puissance des 
modules solaires a mettre en place, la capacite des batteries et des caracteristiques des 
convertisseurs a utiliser pour satisfaire les besoins en energie electrique. 

Le dimensionnement d'une telle installation solaire se fait en sept grandes etapes : 
(Jimmy et al. 1998, Gilbert et al. 2007). 

1. Estimation des besoins energetiques 

L'energie a produire a 1’ installation selon les besoins de chaque composant, tel que (le 
compresseur, le ventilateur de condenseur, la lampe interne, l’onduleur et le regulateur). En 
multipliant la puissance electrique par le temps de fonctionnement journalier, on obtient les 
consommations journalieres. 


2. Estimation de l'ensoleillement 

Pour etre certain de disposer d'energie en toute saison, les calculs sont faits dans les 
conditions d'ensoleillement les moins favorables de la periode d'utilisation. Actuellement, en 
Algerie, la mesure des donnees du rayonnement solaire est effectuee d’une maniere 
instantanee et sur certains sites seulement de l’Algerie. Le calcul du rayonnement moyen 
global journalier refu sur un plan incline ( RGJI) est base sur l'integration de la fonction de 
RGJI. La duree d’ensoleillement a ete definie par l'Organisation meteorologique mondiale 
comme le temps pendant lequel le rayonnement solaire depasse le niveau de 120 W/m 2 , 
(Kalogirou. 2009). 


RIGJ = f 


heure de coucher du soleil 


G(t). dt 


heure de lever du soleil 


(4.1.1) 
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3. Estimation de la capacite de stockage 

Les batteries a usage solaire autorisent une decharge de 80% contre environ 50% pour 
les batteries classiques au plomb. Ce qui fait qu’une batterie de 100 Ah ne represente qu’une 
capacite de 80Ah. Les batteries a usage solaire autorisent une decharge maximale de 80% : 
une batterie de 100 Ah foumit une capacite maximale de 80Ah. A une tension elevee lors de la 
charge d'une batterie, les reactions electrochimiques deviennent tres rapides et donnent lieu a 
un degagement de gaz a l'interieur de la batterie : c'est la tension de gazeification. H faut 
maintenir un degre de ventilation car une concentration d'hydrogene superieure a 4% dans un 
local de batteries presente des dangers d'explosion. A 25°C, quand la batterie est dechargee, la 
densite specifique de l'electrolyte est dans la gamme de 1.150. Si la batterie est entierement 
dechargee, l'electrolyte est essentiellement de l'eau avec une densite specifique de 1. On 
mesure la densite specifique de l'electrolyte d'une batterie ouverte grace a un densimetre ou 
pese-acide, (Emmanuel. 2008). 

La formule simplifiee du pare batteries est, (Emmanuel. 2008): 


C = 


_ ( NjEj ) 


(dpV) 


(4.1.2) 


Ou C : Capacite du pare batterie (Ah) ; Nj: Nombre de jours d’autonomie desire; d p : 
Coefficient de decharge profonde des batteries solaires ( d p =0.8 ); V: Tension de 
fonctionnement nominal (V); E } : Consommation journaliere (Wh/jour). 


4. Estimation du champ photovoltaique requis 

Les pertes sont inherentes a tout processus de conversion d'energie. Les systemes 
photovoltaiques doivent fournir toute l'energie y compris celle qui est perdue. Ces pertes ont 
plusieurs origines (Dialgui. 2009): 

• Pertes par accumulation de poussiere sur les modules ; 

• Pertes dues a l'echauffement et a la variation de l'eclairement des modules ; 

• Chutes de tension dans le cablage et les equipements de conditionnement de l'energie 
(boite de derivation, regulateur) et les pertes dues aux batteries. 

La formule simplifiee du champ solaire est (Dialgui. 2009): 


Pr = 


(0.6 xEj) 


(4.1.3) 


2 ,- 


Ou P c : Puissance crete du champ solaire (W c ) ; E , : Gisement solaire (kWh/m /jour). 
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5. Dimensionnement des cables electriques solaires 

a) Resistance d'un cable electrique 

La resistance d'un cable electrique ne depend ni de la tension ni de l'intensite du courant 
qui le traverse, mais depend de la resistivite ip) du materiau utilise (cuivre, argent, fer, ...), de 
la longueur du cable, de sa section, et de la temperature. Le cuivre est de loin le conducteur le 
plus utilise, et sa resistivite oscille entre 16xl0" 9 a 0°C et 17xl0" 9 Q.ma 25°C. L'equation 


permettant de connaitre la resistance est la suivante (Dialgui. 2009): 

R = p x L/S 


(4.1.4) 


Ou R: Resistance (Q); p : Resistivite (O.m); L : Longueur du cable (m); S : Section du 
cable (m 2 ). 


b) Pertes dans un cable electrique "Pertes Joules" 

Le cable a une resistance, une partie de l'electricite qu'il transporte se transforme en 
chaleur, comme dans un radiateur. II se produit done une baisse de tension qui peut poser 
probleme. C'est ce que l'on appelle "l'effet joule". La resistance du cable vient egalement 
s'aj outer aux resistances des autres appareils de l'installation (resistance interne de la batterie 
par exemple). Plus ces resistances sont elevees, moins le courant pourra facilement circuler. 
En plus d'une chute de tension, il y aura done une reduction de l'intensite. En courant continu, 
comme en courant alternatif, l'equation permettant de connaitre les pertes est la 
suivante (Dialgui. 2009): 

P = Rxl 2 (4.1.5) 

Ou P : Perte (W); R : Resistance (Q) ; I : Intensite en (A). 


c) Reduire les pertes dans les cables 

Dans une installation photovoltaique, plusieurs solutions permettent d'ameliorer 
l'efficacite des cables (Dialgui. 2009) : 

• La premiere solution est de prendre une section de cable plus grosse, mais cela peut 
avoir un cout. 

• La deuxieme solution, c'est de faire en sorte que l'intensite du courant ne soit pas trop 
elevee dans le cable. Pour une meme puissance, reduire l'intensite consiste a 
augmenter la tension. Autrement dit, il faut mettre les modules en serie au lieu de les 
mettre en parallele. Dans ce cas, ils ne produisent plus 12V mais 24 voir 48V. Le 
probleme, c'est que le 24 et le 48V ne sont pas directement exploitables contrairement 
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au 12Y avec lequel certains appareils fonctionnent. De plus, les autres elements de 
l'installation (regulateur, batterie, . ..) utilisant ces tensions sont plus rares, et done 
legerement plus chers. 

• La troisieme est de reduire les distances le plus possible, que l'electricite soit produite, 
stockee, transformee et consommee au meme endroit. La distance entre le champ 
solaire et le local technique abritant la regulation et les batteries doit etre inferieure a 
10 m (jusqu’a 10 kW c ). Pour les generateurs a faible puissance (< 600 W c ) cette 
distance doit etre inferieure a 5 m. 

• Enfin, si les longueurs de cable sont trop importantes, le mieux est alors de regrouper 
les modules, le regulateur, les batteries et l’onduleur dans un meme endroit afin de 
transporter l'electricite en courant alternatif. 


6. Dimensionnement d'un regulateur solaire 

Le regulateur a pour fonction de stopper la charge des batteries lorsque celles-ci sont 
deja completement chargees et aussi de couper L alimentation des appareils lorsque la charge 
des batteries devient critique 20% et evitent ainsi la surcharge et la decharge profonde. Ces 
deux phenomenes sont prejudiciables aux batteries. 

Le regulateur comporte trois circuits entre (panneaux solaire), sortie (appareils 
d’utilisation), circuit de batteries. La caracterisation du regulateur se fait par les deux circuits 
d’ entree et de sortie. L’intensite maximale admissible par le circuit d’ entree doit etre 
superieure a l’intensite fournie par les modules solaires, meme chose pour le circuit de sortie. 

Entree: puissance des modules 

Sortie: nous estimerons que tous les appareils fonctionnent en meme temps (dans notre 
cas on rappelle qu’il y a que le refrigerateur). 

Les modules foumissent de l’electricite sous une certaine tension (12, 24 ou meme 
parfois 48VDC). C’est ce que l’on appelle la tension d’entree. On peut egalement trouver, 
dans leur fiche technique, l’intensite maximale qu’ils sont capables de delivrer. II faut done, 
en premier lieu, choisir un regulateur fonctionnant sous cette meme tension, et qui puisse 
supporter l’intensite maximale d’entree. Si vous optez pour un regulateur controlant 
egalement la decharge, il faut se renseigner sur l’intensite dont vos appareils ont besoin, e'est- 
a-dire l’intensite maximale de sortie. Parmi les onduleurs les plus utilises, on trouve Shunt ou 
MPPT (Yighetti. 2010), Le premier le plus courant car il offre un bon rapport qualite/prix. H 
assure une recharge complete de la batterie, par contre le MPPT qui signifie Maximum Power 
Point Tracker, autrement dit le "detecteur de point de puissance maximum". La tension des 
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modules varie en fonction de la temperature : plus il fait chaud et plus la tension diminue. 
Sachant qu’un module peut atteindre facilement 50 a 60°C meme en etant ventile, les 
constructeurs les ont done surdimensionne afin que la tension soit toujours suffisamment 
elevee. Cependant, si le module delivre une tension de 17V, par exemple, pour alimenter une 
batterie 12V, les 5V de difference seront perdus. un regulateur MPPT transforme le surplus de 
tension en intensite, done il reduit la tension lorsque c’est necessaire, tout en augmentant 
l’intensite pour maintenir la puissance du systeme au maximum. Le cout de cet appareil et sa 
consommation d’energie ne le rendent pas avantageux pour des installations de petites tailles. 

7. Dimensionnement du convertisseur (Onduleur) 

Le convertisseur de courant se dimensionne en fonction de plusieurs criteres (Vighetti. 

2010 ): 

a) La tension d'entree : c'est la meme que la tension des batteries ou du regulateur (12, 
24 ou 48V DC). 

b) La tension de sortie : en Algerie nous utilisons du 220/230 VAC, 50Hz. 

c) La puissance nominale : C'est la puissance dont les appareils qui consomment votre 
electricite ont besoin pour fonctionner de fa§on "normale". Pour connaitre cette 
puissance nominale, il suffit de faire la somme des puissances des appareils 
electriques. Il faut toujours choisir un convertisseur dont la puissance est legerement 
superieure a celle des appareils. 

d) La puissance maximale : L'onduleur est capable de fournir une grande puissance sur 
un court laps de temps (generalement 2 ou 3 fois la puissance nominale). Cette 
particularite est utile pour les appareils qui possedent un moteur (refrigerateur, micro- 
onde, lave-linge, ...), car leur consommation augmente tres fortement lors du 
demurrage. En general, si vous avez correctement dimensionne la puissance nominale, 
la puissance maximale est suffisante. 

e) Le rendement : Une partie de l'electricite transformee est consommee par le 
convertisseur de courant (entre 80 et 95% de l'energie est restituee). Il est important de 
controler ce rendement, sachant qu'un bon produit se situe autour de 90%. De plus, la 
plupart des convertisseurs consomment de l'energie meme lorsqu'ils ne fonctionnent 
pas. Heureusement, certains sont equipes d'un systeme de marche/arret qui permet de 
grandes economies dans les petites installations photovoltaiques. 

f) La forme du signal : La sinusoide (l'onde) produite par le convertisseur peut avoir 
plusieurs formes (de la plus couteuse a la moins chere : pur sinus, pseudo sinus ou 
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carre). Certains appareils supportent tres bien les ondes carrees, mais d'autres non. En 
revanche, n'importe quel appareil peut fonctionner avec une onde pur sinus. 

Vu de la complexity de prevision des conditions meteorologiques et des performances 
des composants de l'installation solaire photovoltaique, des methodes informatiques ont ete 
mises au point, ces methodes de dimensionnement fournissent directement une installation 
optimisee, c'est-a-dire la plus petite installation en modules et en batteries qui repond aux 
exigences du taux de couverture solaire voulu et a la profondeur de decharge maximale 
admissible par les batteries et respectant le fonctionnement des eventuels onduleurs du 
generateur photovoltaique. On peut citer des logiciels comme PVSyst, PV Design, etc. 

L'optimisation du dimensionnement de nos installations de refrigeration solaire 
photovoltaique a ete effectuee avec le logiciel de dimensionnement de systemes et 
d'installations solaires photovoltaiques PVSyst V5.59. 

4.I.I.2. Logiciel PVSyst 

Le logiciel PVSyst est un logiciel de simulation et de dimensionnement d'installations 
solaires photovoltaiques autonomes et relies au reseau. Ce logiciel a ete elabore par 
l'Universite de Geneve (en Suisse). Le logiciel PVSyst dispose de plusieurs entrees : entree 
flux solaires moyens mensuels, temperatures moyennes mensuelles, besoins energetiques, 
choix de modules PV et de leur inclinaison, choix des batteries, regulateurs de charge, 
onduleurs, entree du nombre de jours d'autonomie, du taux de couverture solaire et du cout 
d'investissement (achat materiel, cout d'installation du systeme). Les principaux resultats de la 
simulation sont : la puissance du champ requis, la capacite de stockage, les caracteristiques 
des composants sous des conditions precises et le cout de revient du kilowattheure (kWh). 

Ce logiciel se decompose en 2 parties. La premiere sert a faire un dimensionnement 
rapide de l'installation a partir de la meteo, de la charge, de l'orientation des panneaux, de la 
profondeur de decharge maximale des batteries et du taux de couverture voulu. On obtient 
alors la puissance totale des panneaux et la capacite des batteries a installer. Dans la seconde 
partie une simulation du fonctionnement de l'installation solaire peut etre faite a partir de ces 
premiers resultats de simulation qui prend en compte des facteurs comme les pertes de 
puissances dues au regulateur, aux cables. Le logiciel fonctionne aussi bien en simulation 
qu'en pre dimensionnement. II dispose de sorties detaillees et de base de donnees des 
caracteristiques de materiel (capteurs, batteries, onduleurs, ..), (Manuel PVSyst de cours 
photovoltaique, 1996). Les etapes de simulation par PVSyst sont illustrees dans 
1’ organigramme suivant, (figure 4.1.2). 
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Entrees 

1. Conditions meteorologiques 

■ Irradiation moyenne mensuelle globale re§ue sur une surface horizontale (kWh/j) 

■ Irradiation diffuse moyenne mensuelle re§ue sur une surface horizontale (kWh/j) 

■ Temperature ambiante moyenne mensuelle (°C) 

■ Vitesse de vent moyenne mensuelle (m/s) et l’azimut 

2. Coordonnees geographiques 

■ Latitude 

■ Longitude 

■ Altitude 

3. Consommation moyenne mensuelle du refrigerateur (KWh/j) 

4. Caracteristiques de batterie solaire 

■ Marque de batterie solaire 

■ Nombre de jours d’autonomie (03 jours) 

■ Tension des batteries/systeme (24 V) 

5. Caracteristiques de champs photovoltaic! ues 

■ Marque de module photovolta'ique et leur rendement 

■ Orientation et angle d’inclinaison 

■ Tension de champs photovolta'ique requis (24 V) 

6. Type (marque) du regulateur solaire de charge/decharge 

7. Type (marque) de l’onduleur solaire DC-AC 

8. Coordonnees geographiques 


$ 


■ Puissance crete et 1’ assemblage possible des modules pour 
avoir la tension et le courant de systeme. 

■ Capacite de stockage (Ah) avec T assemblage possible. 

■ Caracteristique de l’onduleur solaire utilise (Courant, tension et 
puissance). 

■ Caracteristique du regulateur de charge/decharge utilise 
(Courant, Tension et puissance). 

■ Cout d’investissement (achat et installation) 

■ Cout energetique (EUR/kWh) 

■ La section des cables d’ interconnexion et la distance entre le 
champ solaire et le local technique abritant la regulation et les 
batteries. 


Dimensionnement de systeme solaire par PVsyst V5 avec 
integration d’une base de donnees meteorologique soit generee de 
Meteonorm V6 soit enregistree dans la station radiometrique. 

5 

Sorties 
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Figure 4.1.2: Methodologie de dimensionnement d’un systeme solaire photovolta'ique autonome avec stockage electrochimique (Manuel PVSyst de cours 

photovolta'ique, 1996). 
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4.I.I.3. Caracteristiques de systeme de refrigeration DC-AC 


a) Dispositif experimental 

Un refrigerateur classique de marque ENIEM de 160 litres a ete choisi pour etre 
alimente en energie solaire. II est compose de deux parties distinctes, une partie de 120 litres 
assure la refrigeration et une autre de 40 litres assure la congelation, il fonctionne avec un 
compresseur de puissance moyenne 90W, une alimentation de 220-230/50/Hz, il consomme 
0.65kWh/j a vide (237kWh/an) (Figure 4.1.3 et Figure 4.1.4). 



Figure 4.1.3 : Refrigerateur a compression de vapeur de 160 litres. 



Figure 4.1.4 : Groupe de compression du refrigerateur de 1601itres 


4.1.2. Optimisation du dimensionnement de l’installation frigorifique solaire DC-AC 
a) Optimisation du dimensionnement par PVSyst 

Fa consommation energetique a vide de ce refrigerateur est 0.65kWh/j (puissance du 
compresseur est 90W et la puissance de la lampe est 15W). Apres un test de fonctionnement 
sous une charge frigorifique de 12 litres d’eau, le refrigerateur a fonctionne pendant 9h/j en 
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simulant le temps d’ouverture a lh/jour, la consommation du regulateur et l’onduleur est de 
40W. Le bilan energetique de ce refrigerateur est : 

CEJ = (90 x 9) + (15 x 1) + 40 = (865 Wh/j) (4.1.6) 

Le tableau 4.1.1 illustre les resultats de dimensionnement du systeme photovoltaique 
DC-AC sous les conditions climatiques mediterraneenne de Bou-Ismail et saharien de 
Tamanrasset. 


Tableau 4.1.1: Dimensionnement du systeme de refrigeration solaire photovoltaique DC-AC 


Region 

Coordonnees 

geographiques 

Puissance 

du 

Champ PV 
(W C /24V) 

Capacitede 
stockage (Ah / 24V) 
pour 03 jours 
d’ autononlie. 

Cout global 
du 

systeme (DA) 

UDES - Bou- 
Ismial (Tipaza) 

Latitude : 36.34 
Longitude : 02.57 
Altitude : 5.00 m 

300 

119 

463872.61 

Bechar 
(Centre ville) 

Latitude : 22.47 
Longitude : 05.31 
Altitude: 1377 m 

200 

119 

323919.45 


Le generateur photovoltaique de 300Wc est compose de deux panneaux de 150Wc-12V 
ou (04 panneaux de 75 We- 12V) relies en serie. Quatre batteries anciennes de 100 Ah 
chaqu’une avec un montage hybride (serie - parallele), ce qui offre une capacite 
surdimensionnee de 200 Ah. Pour dimensionner le regulateur on calcule son courant d’ entree 
et de sortie. Le champ photovoltaique genere un courant de 300/24= 12. 5 A. La puissance 
totale des appareils est: 90+15+40=125W done le courant de sortie est 145/24=6A. Nous 
prendrons done un regulateur de 20 A/20 A (la sortie sera largement sous dimensionnee). En 
fin, un onduleur de marque Victron Energy de 750W et de 20A a ete selectionne, il foumit a 
partir d'une tension continue 24VDC, une tension alternative de 230V/50Hz. Le choix de tous 
ces equipements est porte sur le rapport de la disponibilite sur le marche Algerien et le prix. 
Le systeme de regulation est schematise sur la figure 4.1.5. 



Figure 4.1.5: Regulation et conversion du systeme de refrigeration DC-AC. 
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Le principe de montage en parallele consiste a additionner les courants (I) et de garder 
le meme voltage (VDC) du systeme, par contre le montage en serie consiste a additionner le 
voltage du systeme (VDC) et de garder le meme courant (I). Pour un montage en parallele il 
faut toujours connecter ensemble les polarites de la meme phase, par contre il faut toujours 
connecter ensemble les polarites plus (+) et moins (-). 


Tableau 4.1.2 : Caracteristiques des equipements solaires utilises pour le systeme DC-AC 


References 


Regulateur 
solaire STECA 
12/24 DC 20 A 


Onduleur 


Quelques caracteristiques generates et 

techniques 

Pour panneaux solaires 200W a 500W. 
Regulation de tension et de courant 
Adaptation automatique au niveau de 
tension 12 V / 24 V 
Efficacite maximale : >98 % 

Consommation propre : 10 mA 

Tension sinusoidale pure 
Puissance: 750 W/ 24 VDC / 220VAC 
Puissance consommee a vide : 14 W 
Efficacite max : 89 % 

Tension / Frequence de sortie: 230 V 


VICTRON 


±3 % / 50 Hz ±0.1 %. 


Phoenix 


- Plage de tension d'entree: 10.5 - 15.5 


Photos 




VDC. 


Panneau solaire 
photovoltaique 
polycristallin, 


IS4000P- 12- 150- 
POLY 


Temperature de fonctionnement: de - 
40 a +50°C (refroidissement par 

ventilateur). 

Puissance max : 150W C . 

Gamme de temperature : de - 40 a + 
85°C 

Dimensions et poids : 

1490 x675x42mm et 13kg. 

Courant a puissance max, (Imp) : 

7.98A 

Courant de court-circuit, (Isc) : 8.54A 
Tension en circuit ouvert, (Voc)=22.5 




- Tension a puissance maximale, (Vmp) 
=18. 8V 

- Rendement : 15% 


b) Methodologie et techniques de mesure de la temperature 

Pour pouvoir negliger les fluctuations de temperatures a travers l’epaisseur de la 
conduite, il est necessaire que le nombre de Biot soit petit (Battaglia et al. 2010 ; John. 2005) : 



hxe 

~ 


Resistance interne ( Conduction ) 
Resistance externe ( Convection ) 


« i 


(4.1.7) 
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Avec h-10 W/m 2 .°K le coefficient de convection entre l’air et le cuivre sous une 
temperature ambiante comprise entre 15-30°C pour la convection libre, en general, il est 
suppose que h reste constant avec le temps, e est l’epaisseur de la conduite et A c’est la 
conductivity thermique du cuivre (380 W/m.°K). Le nombre Biot est calcule pour une 
conduite de l A inch (6.35 mm) et une epaisseur de 1mm, la valeur de ce dernier est done de 
l’ordre de 3xlO" 5 . 

D’ apres la valeur de Biot trouvee, on peut conclure que la resistance interne est 
negligeable et la conduite est supposee thermiquement mince, raison pour laquelle, la mesure 
de la temperature par les thermocouples a ete faite en collant la sonde des thermocouples soit 
a la surface exterieure soit en les injectant dans les tubes capillaires soudes apres pelage de la 
conduite frigorifique, mais le soudage du capillaire a l’interieur de la conduite est difficile, car 
le diametre est trop petit, de plus ce capillaire va causer une obstruction a l'interieur du circuit 
de gaz refrigerant (particule de cuivre, de soudures et le capillaire). 

Afin d’evaluer la performance du refrigerateur solaire photovoltaique DC- AC, une 
acquisition de donnees a ete effectuee par un Data logger (Hydra 2700), relie a un ordinateur 
via une interface graphique de communication. Le Data logger enregistre les temperatures 
mesurees par les thermocouples de type K (de composition Chrome/ Alumel possedant une 
gamme de temperatures qui s’etale de -250 a 1372 °C) apres un regime etabli et cela toutes les 
30 secondes. Les donnees ainsi recoltees sont sauvegardees sous forme de fichier Excel. II 
permet aussi un suivi instantane de 1’ evolution des temperatures a l’aide d’un trace de celles- 
ci en fonction du temps de marche. Ces thermocouples ont ete etalonnes avant d’etre utilises, 
et la sonde a ete bien isolee et collee sur la surface exterieure de la tuyauterie (Voir la figure). 

Afin d’evaluer la performance du refrigerateur solaire photovoltaique DC- AC, une 
acquisition de donnees a ete effectuee par un Data logger (Hydra 2700), relie a un ordinateur 
via une interface graphique de communication. Le Data logger enregistre les temperatures 
mesurees par les thermocouples et cela toutes les 30 secondes. Les donnees ainsi recoltees 
sont sauvegardees sous forme de fichier Excel. 
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Conduile de 'A inch et 1mm d’tpaisseur 


Figure 4.1.6: Mise en place des thermocouples 
La mesure des temperatures s’est operee avec des thermocouples de type K de 
composition Chrome/ Alumel possedant une gamme de temperatures qui s’etale de -250 a 
1372 °C, (Figure 4.1.7). Ces thermocouples ont ete etalonnes avant d’etre utilises. 



E3 

Nom Hydra: Hydra 1 

Description: (Aucun) 


Communication: COM1 

Type de declenchement: Lecture par intervalle 

Intervalle de lecture: 5 sec. 

Lectures en memoire: Toutes les lectures 

Destination des donnees: Memoire 

Mode de stockage: Bouclage en fin 

Voie moniteur: Aucun 

Verrouillage panneau avant: Desactiver les verrouillages 
Vitesse de mesure: Lent 


Fichier de donnees: umbbl .csv 

Format du fichier: ASCII 

Lectures enregistrees: Toutes les lectures 

Mode du fichier: Ajouter 

Rollover automatiquement: Aucun 

Espace libre minimum: ND 

Unite temperature: Celsius 

Detection TCO: Oui 

Totalisateur antirebond: Oui 


Voie 

Fonction 

Gamme 

Alarme 1 

Alarme 2 

Mx+B 

Unite 

Etiquette 

□ 

ARRET 







1 

VDC 

1 50/300 V 

ARRET 

ARRET 

ARRET 

vcc 

V champs 

2 

VDC 

1 50/300 V 

ARRET 

ARRET 

ARRET 

vcc 

V batterie 

3 

VDC 

1 50/300 V 

ARRET 

ARRET 

MARCHE 

vcc 

capteur champs 

4 

VDC 

1 50/300 V 

ARRET 

ARRET 

MARCHE 

VDC 

Capt Batt Reg 

5 

VDC 

1 50/300 V 

ARRET 

ARRET 

MARCHE 

VDC 

Cap Batt Ond 

6 

TC 

K 

ARRET 

ARRET 

ARRET 

°C 

Compresseur 

7 

TC 

J 

ARRET 

ARRET 

ARRET 

°C 

Sort condenseur 

8 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Aspiration 

9 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Refoulement 

10 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Milieu Cond 

1 1 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Ambiante 

12 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Thermostat 

13 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Entree evapor 

14 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Sortie evaporat 

15 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Evaporateur 

16 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Evap en dessous 

17 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Haut chambre 

18 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Mil Chambre Eau 

19 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

Bas Chambre 

20 

TC 

K 

ARRET 

ARRET 

ARRET 

°c 

T PV 


Fluke Hydra Logger - . 




-a- -b- 

Figure 4.1.7:a : Acquisition de donnees en utilisant Hydra 2700 Datalogger avec un 
enregistrement toutes les 30 secondes, b : Data Logger et avec carte de 20 thermocouples 


Pour la mesure du courant on utilise un capteur de courant a effet Hall, il exploite l'effet 
Hall (un courant electrique traversant un materiau baignant dans un champ magnetique 
engendre une tension perpendiculaire a ceux-ci) pour produire une tension qui est l'image 
exacte du courant a mesurer ou a visualiser (Figure 4.1.8). 



Figure 4.1.8: a : Thermocouple de type K, b : Capteur de courant et c : et alimentation DC 

stabilise. 
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La tension est mesuree directement par le Data Logger grace a un simple branchement 
sur la boite d’ acquisition et une configuration sur 1’ interface graphique. Enfin la 
consommation du refrigerateur est controlee par un compteur d’energie place entre la sortie de 
l’onduleur et le refrigerateur qui nous indique l’energie consommee en (kWh), (Figure 4.1.9). 



Figure 4.1.9: a: Emplacement des capteurs de courant et de tension, b: Compteur d’energie. 

L’eclairement solaire sur le plan d’ouverture du champ photovolta'ique est mesure par 
un Pyranometre de type CIMEL 5.31mV/W/m 2 , sur le site de Bou-Ismail (inclinaison : 36 ° et 
L orientation vers le sud du pays), par contre celui de Tamanrasset a ete estime par le Logiciel 
meteonorm'9 qui est integre dans la base de donnees de PVSyst. La duree d’ensoleillement a 
ete definie par l'Organisation meteorologique mondiale comme le temps pendant lequel le 
rayonnement solaire depasse le niveau de 120W/m 2 , (Solar Energy Engineering. 2009). 

c) Calcul de la surchauffe et de sous-refroidissement 

Le but de ce test est la verification de la surchauffe et le sous-refroidissement, ces 
derniers sont des renseignements importants sur les conditions de fonctionnement du circuit 
frigorifique en particulier le compresseur. La surchauffe indique 1’ ecart de temperature du gaz 
refrigerant a la sortie de l’evaporateur (gaz) et la temperature d’ evaporation lue au manometre 
BP, (relations : pression, temperature). Cet ecart represente la qualite de remplissage de 
l’evaporateur. Une surchauffe convenable est comprise entre 5°C a 8°C. Le sous- 
refroidissement est 1’ ecart de temperature du refrigerant liquide a la sortie du condenseur et la 
temperature de condensation. Cet ecart represente la qualite de remplissage du circuit 
frigorifique. Un sous-refroidissement convenable est compris entre 4°C a 7°C. 

D'apres les resultats de ce test, les valeurs de la surchauffe et le sous-refroidissement 
sont bien confondues aux valeurs trouvees dans la litterature. 
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LUMIfRE 

(PHOTONS) 


Pyranometre 



Figure 4.1.10: Illustration schematique du systeme experimental pour un refrigerateur solaire 

PVAC 

d) Bilan energetique du cycle frigorifique et coefficient de performance (COP) 

Les tests ont ete conduits sur mois sous differentes position de thermostat avec et sans 
charge frigorifique. Les temperatures des differents organes de notre refrigerateur relatives a 
la premiere position du thermostat et pour un fonctionnement sans ouverture de la porte et 
avec charge ont ete enregistrees. Ensuite, on a evalue les valeurs d’ enthalpies relatives aux 
temperatures mesurees du fluide frigorigene utilise (R134a), en tra§ant le cycle dans le 
diagramme de refrigerant, afin de caracteriser notre machine frigorifique. 
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Dans une premiere etape, on a procede au calcul des valeurs de la surchauffe et du sous 
refroidissement afin de determiner les performances energetiques associees au systeme a 
savoir, le travail fournie par le compresseur, la production frigorifique au sein de 
l’evaporateur et enfin la quantite de chaleur evacuee au niveau du condenseur. Ces valeurs 
energetiques vont nous permettre dans une seconde etape, de determiner le coefficient de 
performance de notre refrigerateur reel (COP^ei) ainsi que le theorique (COP t h). 


e) Diagramme enthalpique par le logiciel CoolPack 

Le logiciel de simulation CoolPack pour les systemes frigorifiques permet une analyse 
du cycle et le dimensionnement des principaux elements du circuit. II permet egalement de 
tracer le cycle frigorifique relatif a P installation etudiee afin de calculer les performances 
energetiques associees. On a introduit les valeurs des temperatures experimentales de 
condensation et de vaporisation respectivement 41.5°C et -19. 8°C, ainsi que la surchauffe et le 
sous refroidissement moyen AT SC =5 °K, AT sr = 4.5 °K ;D’apres P analyse des performances de 
la machine frigorifique durant les deux premiers jours de fonctionnement par le logiciel 
CoolPack donne un COP r e e i de 3.07, alors que le COP t h est donne par la relation suivante : 

COP th — ~~~ (4.1.8) 



Figure 4.1.11: Diagramme enthalpique du refrigerateur DC-AC 


Tableau 4.1.3: Proprietes du cycle du refrigerateur DC-AC 


Production frigorifique : Qe (Kj/Kg) 

144.64 

Chaleur echangee avec le condenseur : Q c (Kj/Kg) 

191.76 

Travail fourni au compresseur, W (Kj/Kg) 

47.13 

Taux d’ aspiration (-) 

8.93 

Surchauffe (°K) 

7 

Sous-refroidissement (°K) 

10 

COPreei de Pinstallation (-) 

3.07 
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f) Etude du comportement du compresseur 

Afin d’etudier le fonctionnement de ce refrigerateur, on a presente a titre d’exemple la 
variation de la puissance consommee par le compresseur au cours de fonctionnement. La 
puissance donnee par le fournisseur est de 90W mais elle est variable selon 1’ arret et le 
fonctionnement du compresseur. Au demarrage il a besoin d’une puissance moyenne de 125W 
et un courant de 5 A. La puissance moyenne lors du fonctionnement du refrigerateur en charge 
maximale est estimee a 115W.ha Figure 4.1.10 nous montre la variation de la puissance de 
fonctionnement du compresseur dans le temps. 


Puissance (W) 


Puissance die compresseur 



emps (mu) 


Figure 4.1.12: Variation de la puissance du compresseur du systeme DC- AC en fonction de 

temps. 


La temperature est un parametre important dans le comportement des modules, une 
cellule exposee a une irradiance de lkW/m 2 ne transforme en electricite que 12% au plus, le 
reste soit 88 % etant dissipe en chaleur et en reflectance. Ainsi, si les modules ne sont pas 
correctement ventiles, les cellules peuvent monter tres haut en temperature et voir leurs 
performances electriques se degrader. II est done important de verifier que vos panneaux 
photovoltaiques soient bien ventiles et qu’ils soient en mesure de supporter des montees en 
temperature sans des baisses de rendement importantes. 

D’apres la figure 4.1.13 la puissance generee par notre champ n’est pas stable. Afin de 
tirer le maximum, il est preferable d’utiliser des regulateurs integrants la technologie MPPT 
(Maximum Power Point Tracking) qui va extraire toute la puissance du champ solaire et 
augmente son performance jusqu'a 30% par rapport a un regulateur conventionnel. 
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Temps 


Figure 4.1.13: Puissance generee par le champ photovoltaique du systeme DC-AC 


4.1.3. Optimisation du dimensionnement d’un systeme frigorifique solaire 

direct DC-DC avec accumulation de froid 

4.I.3.I. Dispositif experimental 

Partant d’un refrigerateur 1601 de marque ENIEM, on precede au remplacement de tous 
les organes, compresseur, condenseur et evaporateur, en gardent la chambre frigorifique. Le 
choix du compartiment a refrigerer de volume de 1601 est base sur les caracteristiques 
techniques du groupe frigorigene utilise et 1’ evaporateur adequat a cet equipement. 
L’ adaptation de ces equipements au volume etant donne par le fabriquant (tableau 4.1. 4). 




Figure 4.1.14: Schema du refrigerateur modifie a) Evaporateur eutectique, b) kit BD50F 

compresseur/condenseur. 
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Type de l’evaporateur 

Volume 

maximale 

k 

refrigerer 

(Litres) 

Type de compresseur k air 

BD35F 

BD50F 

GD30 

2000RPM 

2500 RPM 

3000RPM 

2000RPM 

2500 RPM 

3000RPM 

3500RPM 

2000RPM 

2500RPM 

3000RPM 

3500RPM 

Evaporateur du type 
rat 

PT1 

40 

X 







X 




PT2 

60 


X 


X 





X 



PT3 

r 80 


X 


X 





X 



PT4 

130 


X 


X 





X 



PT5 

150 


X 


X 





X 



PT13 

170 



X 


X 




X 



PT8 

210 






X 




X 


PT7 

300 







X 




X 

c 

U 

3 

2 i 

o 

CL 

3 
> 

O 

SI 

80 


X 


X 





X 



S3 

150 


X 


X 





X 



S6 

210 






X 




X 


evaporateurs a 
plaques 

AC1 

40 


X 


X 





X 



AC5 

60 


X 


X 





X 



AC6 

80 



X 


X 




X 



AC 10 

100 






X 




X 


AC 15 

140 






X 




X 


AC20 

170 







X 




X 

AC30 

230 







X 




X 

3 

£ *01 
3 S 

© c 
a. % 

3 ^ 

> 

AIR10 

200 



X 


X 




X 




Tableau 4.1.4: Fiche technique utilisee pour 1’ adaptation des equipements frigorifiques au 

comportement a refrigerer. 


4.I.3.2. Kit BD50F de compression/condensation et evaporateur eutectique 
a) Compresseur 

Le compresseur DANFOSS BD50FF utilise a une puissance de 59W et une cylindree de 
2.50CC avec un rotor a aimants permanents alimente directement en courant continu 
12V/24V et un courant de 4.9A/2.45A. II fonctionne avec le R134a (0.050 kg), mais peut 
egalement etre utile avec le fluide frigorifique R12. Le passage du 12Y jusqu’au 24V se fait 
par le changement d’un fusible sur le boitier ou F unite electronique entierement integree avec 
une commande de variation de vitesse couvrant une plage de 2000-3500tr/min avec un pas de 
variation de 500 tr/min (Figure 4.1.15). 
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Figure 4.1.15: Unite electronique du compresseur BD50F 
Ce compresseur etant dote d’un bas couple de demarrage, ("Low Starting Torque" = 
LST) doit uniquement etre utilises pour les systemes de refrigeration equipes d’une detente 
par tube capillaire, ou l’equilibrage des pressions d’ aspiration et de refoulement se fait durant 
chaque periode d’ arret du compresseur. 


b) Condenseur 

Le systeme a compression de vapeur comporte un condenseur compact enroule avec 
trois rangees. Equipe d’un dispositif de refroidissement par ailettes, avec ventilation 
horizontale d'air assuree par un ventilateur de 12V relie a l’unite electronique. 



Figure 4.1.16: Condenseur compact a ventilation d’air horizontal utilise dans le systeme 

DC-DC 

Le refroidissement se declenche instantanement au demarrage du compresseur, 
refroidissant a la fois le condenseur et la cloche du compresseur, une gris de protection sur le 
cote permet de proteger le condenseur et facilite 1’ aeration au niveau du ventilateur (Figure 
4.1.17). 
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Figure 4.1.17: Groupe de Condensation utilise dans le systeme DC-DC 

Le raccordement entre le systeme a compression et T evaporateur se fait de maniere 
simple et pratique par des raccords speciaux (Figure 4.1.17). Le fonctionnement du 
compresseur est pilote par le thermostat a Tinted eur de la chambre situee sur la face laterale 
gauche de T evaporateur. Ce dernier, declenche la mise en marche du compresseur en fonction 
de la temperature de la chambre pour les differentes positions de 1 jusque a 7. On a adapte 
une lampe de 24V et 5W pour le systeme d’eclairage du refrigerateur. 

c) Evaporateurs a accumulation du froid (eutectiques) 

Un evaporateur a accumulation de froid est un evaporateur enferme dans un caisson 
metallique rempli d’une solution eutectique. Le but est d’accumuler le froid et de le restituer a 
T arret du compresseur. La particularite de la solution eutectique est d’ avoir un ecart de 
temperature de 6 ou 7 degres entre le moment ou elle est congelee et le moment ou elle va 
decongeler et d’etre non corrosive. 6 ou 7 degres Celsius de variation de temperature sont 
acceptables dans un conservateur, mais discutables dans un refrigerateur. La solution utilise et 
de l’Ethylene glycol dilue a 90% (10%Ethylene glycol et 90% d’eau). On utilisee un caisson 
en plaques d’inox polis (325x250x60) mm d’un volume effectif de 2250cm 3 le volume de la 
solution est de 2100cm 3 avec un evaporateur tubulaire en cuivre qui occupe un volume de 
150cm 3 . On a choisi initialement une position verticale de T evaporateur sur la face arriere de 
la chambre du refrigerateur. 

4. 1.3.3. Optimisation du dimensionnement par PVSyst 

La consommation d’energie est de 0.473 kWh/jour sous une charge frigorifique de 12 
litres d’eau et apres un temps de fonctionnement de 7.5h/jour. La lampe DC associe au 
refrigerateur a une puissance de 5W, le temps d’ouverture des portes est estime a lh/jour et la 
consommation du regulateur est de 26W/jour. On aura le bilan energetique suivant : (59 x 
7.5) + (5 x 1) + 26 = (473 Wh/jour). 
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Le tableau 4.1.5 illustre les resultats de dimensionnement du systeme photovoltaique 
DC-DC sous les conditions climatiques mediterraneennes (Bou-Ismail) et saharienne 
(Tamanrasset). 


Tableau 4.1.5: Dimensionnement du systeme de refrigeration solaire photovoltaique DC-DC 


Region 

Coordonnees 

geographiques 

Puissance 

du 

Champ PV 

(W/24V) 

Capacite de 
stockage (Ah / 24V) 
pour 03 jours 
d’autonomie. 

Cout global 
du 

systeme (DA) 

UDES - Bou- 
Ismial (Tipaza) 

Latitude : 36.34 
Longitude : 02.57 
Altitude : 5.00 m 

200 

65 

281098.78 

Bechar 
(Centre ville) 

Latitude : 22.47 
Longitude : 05.31 
Altitude: 1377 m 

100 

65 

195240.61 


Le generateur solaire ou le champ photovoltaique, comportant quatre panneaux solaires 
poly-cristallin de 50W, delivrant un voltage de 12Y par panneau, monte sur un cadre en 
aluminium resistant a la corrosion. Le cadre est place sur un support incline de 36° Sud. Le 
panneau solaire est dote d'une diode anti-retour integree dans la boite de jonction. En se 
basant sur le dimensionnement du systeme de refrigeration, Le montage et le cablage des 
panneaux sont realises de sorte a avoir 24V a la sortie. L’ utilisation d’un montage hybride 
serie parallele nous permet de passer du 12V au 24V et de garantir une puissance de 
4x50=200W c (Figure 4.1.18). 



Figure 4.1.18: Champ photovoltaique necessaire pour alimenter le systeme DC-DC. 


La capacite de stockage pour le systeme DC-DC est de 65Ah/24V. Dans notre cas on a 
opte pour deux batteries solaires de type Surrette Solar de 113Ah et 12V chacune, relie en 
serie pour avoir 24V a la sortie. La partie regulation est assuree par un regulateur ISOFOTON 
de 10 A, offrant ainsi, une protection pour les batteries de la decharge profonde. Afin de 
faciliter le deplacement et de proteger la partie alimentation on a installe les batteries dans un 
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coffre sur lequel on fixe le regulateur, on ajoute une prise monophasique de 24V pour 
alimenter les batteries par le champ photovolta'ique via le regulateur (Figure 4.1.19). 



ISOFOTON 





♦ u • it ♦ wi " n ♦ I 


OO op Qi 

1 1 i *V 


ISOFOTON 




Figure 4.1.19: Stockage et regulation pour le systeme DC-DC. 


4. 1.3.4. Resultats et discussion 

Afin d’etudier le comportement de notre systeme, on a place 20 thermocouples de type 
K d’une precision de 1%. L’emplacement se fait a l’interieur et l’exterieur du refrigerateur et 
sur les organes principaux tels que F evaporateur, condenseur et le compresseur afin de 
calculer le bilan de masse et de l’energie de chaque element du cycle pour calculer le COP, 
(Figure 4.1.20). 



Figure 4.1.20: Differents emplacements des thermocouples dans le refrigerateur (charge 


d’eau de 12 litres). 

Le tableau 4.1.7 montre F evolution de la temperature a l’interieur du refrigerateur en 
fonction de la position du thermostat. Apres fonctionnement du systeme sur un mois, les deux 
valeurs de la surchauffe et de sous refroidissement etaient dans la fourchette proposee par la 
litterature. Apres variation de la position de thermostat, on a constate que la position 3 est 
maximale car apres cette position, la solution eutectique sera completement gelee et va jouer 
le role d’un isolant entre le fluide frigorigene et Fair de la chambre, ce qui permet au 
refrigerateur de consommer plus d’energie. Done, les tests ont ete menes en fixant le 
thermostat sur la troisieme position. 
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Tableau 4.1.6: Representation de la temperature de la chambre et du thermostat. 


Position 

de 

thermostat 

Temperature a 
l’interieur du 
refrigerateur en 
(°C) 

Temps de marche 
du compresseur 

(hears) 

Consommation du 
compresseur et le 
ventilateur (0. 1 8 
A) 

(kWh) 

basse 

haute 

1 

03.50 

09.00 

07h 

0.443 

2 

02.50 

06.00 

06h50min 

0.411 

3 

02.65 

05.30 

07h30min 

0.473 

4 

-03.25 

-00.25 

13h30min 

0.854 

5 

-03.00 

- 01.00 

llh 

0.696 

6 

-05.00 

-02.80 

09h 

0.569 

7 

-08.13 

-05.40 

13h 

0.823 


4.I.3.5. Comparaison entre le systeme de refrigeration DC-AC et le systeme DC-DC avec 
stockage du froid 

A la lumiere de ce travail, deux refrigerateurs de meme volume ont ete choisis pour etre 
alimentes soit directement ou indirectement par energie solaire photovoltaique. Apres une 
phase de dimensionnement de deux systemes, le premier refrigerateur a courant altematif a 
ete indirectement alimente en energie solaire sans modification dans le cycle frigorifique en 
utilisant un convertisseur de tension 24V/230V. En revanche, les equipements du cycle 
frigorifique du second refrigerateur ont ete modifies afin d’etre alimentes directement (sans 
convertisseur de tension) par 1’ energie solaire photovoltaique. A travers les resultats obtenus, 
nous pouvons tirer les conclusions suivantes : 

• Les valeurs de la surchauffe et le sous-refroidissement trouvees ont confirmes le bon 
fonctionnement des deux refrigerateurs solaires photo voltaiques. 

• La simulation de la charge frigorifique a ete faite par l’emploi de 12 litres d’eau car ce 
dernier a la plus grande chaleur specifique. 

• L’ optimisation de dimensionnement de deux refrigerateurs solaires photovoltaiques 
sous deux differents climats (Tipaza et Tamanrasset) a ete faite manuellement et par le 
logiciel PVSyst. De plus le calcul thermodynamique a ete faite par le logiciel 
CoolPack. 

• La consommation journaliere sous le site de Bou-Ismail du systeme frigorifique DC- 
AC est de 0.865kWh/jour (les deux refrigerateurs ont ete places dans le laboratoire 
LFCER-UDES) avec un cout global du systeme environ 463872. 61DA. La 
consommation journaliere du systeme frigorifique DC-DC sous les memes conditions 
et le meme site est de 0.473kWh/jour avec un cout global du systeme environ 
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281098.78 DA, cela induit une reduction de 40%, c.a.d. une reduction de 182773 DA 
par rapport au systeme DC-AC. Cette reduction est due principalement au prix de 
convertisseur de tension. 

• La consommation joumaliere sous le site de Tamanrasset du systeme frigorifique DC- 
AC est de 0.865kWh/jour (Les deux refrigerateurs ont ete places dans un local 
climatise) avec un cout global du systeme environ 323919.45 DA. Le systeme 
frigorifique DC-DC sous les conditions et le meme site est de 0.473kWh/jour avec un 
cout global du systeme environ 195240.61 DA, cela induit une reduction de 40%, 
c.a.d. une reduction de 128678.85 DA par rapport au systeme DC-AC. Cette reduction 
est due principalement au prix de convertisseur de tension et vu que le site jouit d’un 
meilleur ensoleillement. 

• Le COP est le critere le mieux adapte pour la comparaison entre deux refrigerateurs, 
ce dernier est le rapport entre la quantite du froid produite sur 1’ energie donnee au 
compresseur. COP r e e i_Dc-Dc=l-54 COP r eei_Dc-Ac, done en termes de COP le systeme 
DC-DC est mieux adapte pour etre alimente en energie solaire photovoltaique que 
celui en DC-AC. 

• Avec l’utilisation des evaporateurs a accumulation de froid, on peut eviter le stockage 
electrochimique de 1’ energie dans les batteries solaire qui sont considerees les plus 
cheres dans une installation solaire avec une duree de vie moyenne de 5 ans. Ce qui 
traduit par la reduction de cout global du systeme frigorifique solaire DC-DC. 

• Le test du systeme par un pare de batterie afin de determiner la consommation de 
systeme en regime permanant et transitoire, vide et avec charge d’eau, ensuite 
caracteriser chaque equipement dans le but d’ assurer le bon fonctionnement de notre 
prototype et proteger ses equipements contre L endommagement. 

• Les tests sont faits a Tamanrasset dans une caserne de gendarmerie nationale, ce qui 
nous a obliges de ne pas collecter les donnees pour faire une evaluation 
thermodynamique. les conditions metrologiques par "meteonorm 4.49" integre dans le 
logiciel PVSyst. Au cours de la periode de test (Ete 2012), les deux refrigerateurs sont 
bien fonctionnes avec un test d’autonomie de 03 jours sans ensoleillement. 
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4.1.4. Optimisation du dimensionnement d’une chambre frigorifique de 1.5 
m 3 solaire direct DC-AC. 

4.1.4.1. Introduction 

L’objectif principal de ce travail est de fabriquer une chambre frigorifique positive 
(pour la preservation des denrees entre 4 et 8°C) de 1.5m 3 . L’ alimentation de ce dernier est 
basee sur la conversion photovoltaique de l'energie solaire en electricite (ressource gratuite et 
renouvelable). Pour la fabrication de la chambre, des panneaux sandwichs de polyurethane 
ont ete employes pour avoir une bonne isolation thermique afin d’ assurer son bon 
fonctionnement. II sera dimensionne et teste a Bou-Ismail, ensuite il sera transport^ a Ain- 
Guezzam - Tamanrasset pour etre exploite reellement par la Gendarmerie Garde-Frontieres 
(GGF) sous les conditions sahariennes. 

4.1.4.2. Fabrication de la chambre frigorifique 

Apres le dimensionnement de cette chambre frigorifique positif (Voir annexe C) en 
fixant les dimensions interieures de 1210x1055x1255mm avec un volume interieur de 1.56 
m 3 . Le compresseur employe est hermetique de 3/8 cv (275 W) avec condenseur et 
evaporateur a air. L’ isolation des murs exterieurs et plafond a ete realisee par des panneaux 
sandwich polyurethane de 40 kg/m 3 avec une epaisseur de e p =120 mm, en revanche F isolation 
du sol a ete realisee par les panneaux sandwich polyurethane renforces par une tole laquee+ 
bois + resine. La figure 4.1.19 represente une image de la chambre au cours de realisation 
dans un atelier de fabrication mecanique specialise dans le domaine frigorifique a Tizi-Ouzou. 
Le dimensionnement du condenseur doit tenir compte du sous-refroidissement. Le condenseur 
est alors legerement surdimensionne afin d'offrir au fluide frigorigene une surface d'echange 
suffisante pour condenser completement (meilleure performance du cycle frigorifique, 
reduction du risque de "flash gaz" au niveau du detendeur). Le sous-refroidissement 
represente la difference entre la temperature de condensation lue au manometre HP (relations 
: pression, temperature) et la temperature mesuree a la sortie du condenseur. Generalement la 
valeur du sous-refroidissement est comprise entre 5 et 10°C. Un ecart de temperature inferieur 
a 5°C signifie en general un manque de charge, un ecart de temperature superieur a 10°C 
signifie un exces de charge. La valeur de la surchauffe est 7°C, le sous-refroidissement est de 
4°C et le COP egale a 8.16. La figure 4.1.21 resume les differents equipements integres dans 
la chambre. 
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Entrees : La charge quotidienne echangee 

1. Calcul des charges thermiques extemes 

a. Charge thermique par transmission a travers les parois 

b. Charge thermique par renouvellement d’air 

c. Charge thermique par ouverture des portes 

d. Complement d’apport par transmission a travers les parois exposees 
directement au rayonnement solaire 

2. Calcul des charges thermiques internes 

A. Charges thermiques internes independantes des denrees entreposees 

a. Charge thermiques dues a l’eclairage 

b. Charge thermiques dues aux personnes travaillant dans la chambre froide 

c. Charge thermique dues au materiel roulant 

d. Charges thermique due a des machines diverses 

B. Charges thermiques internes dependantes des denrees entreposees 

e. Charges thermiques dues aux denrees entrant a temperature ambiante 

f. Charges thermiques dues aux moteurs des ventilateurs 

g. Charge thermique due a la respiration des denrees 

h. Charge thermique due aux resistances de degivrage 

3. Calcul de la puissance frigorifique en fonction de la charge quotidienne et le temps de 
fonctionnement journalier 



Dimensionnement des equipements frigorifiques soit 
manuellement par les correlations mathematiques existantes soit 
en utilisant des logiciels ou des applications Excel libres. 



Sorties 


■ La puissance du groupe a installer (Compresseur) 

■ La surface totale d’echange de l’evaporateur avec le milieu 
exterieur 

■ La surface totale d’echange thermique du condenseur 

■ Perte de charge du detendeur (A P) 

■ Le separateur de fluide : ses dimensions dependent de la 
puissance frigorifique et la nature de fluide frigorigene 
employe. 

■ Reservoir de liquide 

■ Bouteille d’ aspiration 

■ Les filtres, le voyant liquide et le deshydrater 


Figure 4.1.21: Methodologie de dimensionnement d’une chambre froide positive 
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Figure 4.1.22: Differentes etapes de fabrication de la chambre. 



Figure 4.1.23: Schema synoptique de la chambre realise. 


4.I.4.3. Dimensionnement et test de systeme frigorifique solaire photovoltai'que a 1’UDES 
a) Donnees de rayonnement global incline 

Le rayonnement global regu sur un plan incline est mesure par un Pyranometre de type 
CIMEL 5.31 fiV/W/in 2 incline 36° sud sur le site de Bou-Ismail. La figure 4.1.24 montre un 
exemple de la variation du rayonnement. Le rayonnement global joumalier calculer est de 

l’ordre 3.5 ^ W “ sur la region de test. 

m z .jour 



Figure 4.1.24: Variation journaliere du rayonnement solaire global incline (station de TUDES). 
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Figure 4.1.25: Station meteorologique situee a l’UDES 


Cette station meteorologique est composee de deux Pyranometres pour mesurer le 
rayonnement global sur un plan horizontal et incline, et une console Vantage Pro 2 avec des 
capteurs integres. La station meteo Vantage Pro 2 est composee d'un bloc de capteurs 
comprenant un anemometre, un pluviometre, dans l'abri blanc un capteur d'humidite et de 
temperature, un capteur de rayonnement UV et un capteur de rayonnement solaire, etc... 

b) Calcul de la consommation electrique de la chambre 

La puissance de chaque equipement dans le cycle frigorifique a ete tiree directement de 
la fiche technique correspondante. Apres la phase d’ evaluation de temps de fonctionnement et 
la consommation journaliere du systeme complet en utilisant une source conventionnelle a 
courant altematif et un compteur d’energie sous les conditions climatiques de Bou-Ismail et 
sous une charge frigorifique de 200 kg d’eau. Une etape de dimensionnement de differents 
constituants de 1’ installation solaire photovoltaique DC-AC a ete faite manuellement. Apres 
deux mois d’ analyse et de test, la chambre a ete transportee pour etre teste dans un climat 
aride a Ain-Guezzam (Tamanrasset), (tableau 4.1.8). 


Tableau 4.1.7: Puissance et temps de fonctionnement de chaque equipement dans le 


cycle frigorifique sous test a l’UDES. 


Equipements de la chambre 

Puissance 
unitaire (Watt) 

Duree de 
fonctionnement 
Quotidienne (heure) 

Energie 

consommee (E c ) 
par jour (Wh/j) 

Un compresseur 

3/8 cv = 276 

12 

3312 

Un condenseur 

53 

12 

636 

Un evaporateur 

30 

12 

360 

Un panneau electrique impermeable 

6 

24 

150 

Total 

380 

12 

4458 
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Pour tenir compte des variations meteo et des rendements des differents composants de 
[’installation on majore ce chiffre grace un coefficient correcteur k, en generate la majoration 


est de 50% done k=0.65 . on commence le calcul de l’energie a produire (E p ): 

Eq 4458 


T-f ^ Li A 

hp — — 

r K 0.65 


= 6858.5 W 


(4.1.9) 


Ensuite on calcule la puissance crete P c du generateur photovoltaique necessaire : (On 
prendra une irradiation moyenne de 3.5 kWh/m 2 /jour ) : 

E P 6858.5 


Pr = 


RGJI 


3.5 


= 1960 Wr 


(4.1.10) 


A partir de la puissance crete des panneaux Isofoton IS- 150/24, on determinera le 


nombre de panneaux solaires necessaires a 1’ installation : 

P c I960 


N - 


Puissance crete unitaire dmn panneau 150 


14 Panneaux PV (4.1.11) 


Le calcul de la capacite des accumulateurs necessaires a ce systeme ainsi que leur 


nombre grace a cette formule : 

Ec xN 


c = 


DxU 


4458x3 

0.8x24 


= 690 Ah 


(4.1.12) 


Ou C : capacite de la batterie en ampere-heure (Ah), Ec : energie consommee par jour 
(Wh/j), N : nombre de jour d’autonomie, D : decharge maximale admissible (0.8 pour les 
batteries au plomb), U : tension de la batterie (V). Le pare batterie sera constitue 
d’accumulateur de 2Y chacun, le systeme fonctionnant sous 24 Y les accumulateurs seront 
couples en serie. Le nombre = 24/2 = 12 batteries de 550Ahx2=24 batteries (deux systemes). 
Pour les systemes de petite puissance, une tension de 12V est preferee, les recepteurs 
alimentes en 12V etant plus frequents. Mais au-dela de 2.5 kWh stocke, on passera en 24V ou 
plus pour limiter les pertes ohmiques en gardant des sections de cable raisonnables. Les 


batteries sont commercialisees en cellules de 2V, de 200 Ah a 3000 Ah ou en ensemble 12V 


jusqu’a 400Ah. On evitera si possible de mettre en parallele des batteries. La figure 
4.1.26presente un schema synoptique pour P alimentation electrique autonome de la chambre 
solaire photovoltaique. 
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GENERATEUR l/1050Wc 


> Disjoncteur 60A 

TS 45 


GENERATEUR 2/l050Wc 




BAC1. 24V/550AH 


BAC2 .24V/550AH 


BUS 

24DC 



SORTIE 
220 AC 


CONVERTISSEUR DC/AC 
24V/750VA 


Figure 4.1.26: Schema synoptique du systeme photovoltaique necessaire pour alimenter la 

chambre del. 5 m 3 


Tableau. 4.1.8: Caracteristiques des equipements selectionnes apres dimensionnement 

A * • A • 1 / • A TU A 


Caracteristiques des equipements 


L’evaporateur employe de type GC4040/16, il maintient 
la temperature entre+2 et +10. Dans notre cas on a fixe 
cette gamme entre +4 et +8. Le compresseur employe est 
hermetique de 3/8 cv. 


L’onduleur Top Class 22/24, il transforme le courant 
continu de 24VDC en courant altematif de qualite reseau 
230V AC. Il fournit un courant parfaitement sinusoidal 
d’une puissance de 2000VA. Il a une faible 
consommation a vide. 


Photos 



La batterie solaire a plomb ouvert (OPzS Solar). Les 
batteries employees ont les caracteristiques suivantes : 
Capacite de stockage : C120: 550 Ah/2V 
Nombre de cycle:4000 cycles a 20 % de Cl 20 
Poids acide inclus: 26.7kg ; Resistance interne: 0.63 mO 
Courant de court-circuit: 3250 A 





Le regulateur TriStar TS-45 de Morning star est congu 
pour reguler la charge, il integre trois fonctions 
permettant de charger les batteries avec une grande 
fiabilite, de reguler la charge ou encore la derivation. H 
fonctionne au choix selon l'un de ces trois modes de 
fonctionnement. Le courant nominal solaire est de60A. 



Les panneaux selectionnes de type Isofoton IS 150/24 
possedent quelques caracteristiques : 

Puissance maximale : 150W/24V 

Courant de court-circuit : .45A ; Tension a circuit ouvert 
(V oc ) : 43.2V ; Courant a puissance maximale (I mp ) : 
4.35A et Tension a puissance maximale (I mp ): 43.6V 
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c) Test du courant d’appel et de la puissance 

L’objectif de cette partie est de mesurer le courant instantane de fonctionnement de la 


chambre a l’aide d’un capteur de courant a effet Hall. On a: 


le = — OUlc = - tL - 
b m * 1000 


(4.1.13) 


Et aussi on a : 


Tl .xRxi,-*) 

V = - 

1000 


(4.1.14) 


Ou i s : Courant secondaire et i p : courant primaire, n= 5 spires et la resistance R= 100 Q. Alors 
la tension sera: 


Ce qui donne : 


y = 

1000 


i p = 2xV 


(4.1.15) 

(4.1.16) 



Figure 4.1.27: Implementation des capteurs de courant sur le circuit d’ alimentation. 


La tension efficace de marche de la chambre froide a chaque seconde a ete calculee par 
1’ equation suivante : 

i P eff=2V eff (4.1.17) 

Cette valeur egale environ 4A en marche et environ de 0.02 A a 1’ arret. Cette derniere 
due a la consommation du panneau electrique qui commande la temperature a l’interieur de la 
chambre, (figure 4.1.28). 
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Figure 4.1.28: Evolution du courant en fonction de temps de marche du compresseur 

Une solution serait de mesurer le courant d'appel avec une pince ampere metrique sur 
l'une des phases. La valeur mesuree depend fortement a la difference de la pression entre 
F aspiration et le refoulement au niveau du compresseur. Ce courant d’appel est generalement 
important mais avec un temps de vie tres court, d’ environ 70 jus. 



Figure 4.1.29: Courant d’appel (demarrage) du compresseur 


Vu l’impossibilite d’enregistrer les valeurs du courant par le Data Logger de type Hydra 
car ce dernier il enregistre des valeurs inferieures a 3A seulement, de plus, voir la forme du 
signal a la sortie de l’onduleur en milliseconde a l’aide d’un oscilloscope, d’apres la figure 
4.1.30, il est claire que la forme est sinusoidale avec une periode de 20ms (50 Hertz ) et une 
amplitude de V max = 3 Volts. 

V ‘rr = v -W <4 ' U8) 
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Ou V e ff — 2.12 Volts et i pe ff = 4.24 A 



Figure 4.1.30: Verification graphique de la sortie du l’onduleur 


La valeur mesuree (I pe ff =4.24A ) represente le courant de marche du systeme 
(compresseur, condenseur et evaporateur) mais cette valeur varie a cause de courant de 
demarrage du compresseur. 

Quand il n’y a pas d’ensoleillement, le systeme est alimente a partir des batteries. La 
chambre fonctionne avec un courant moyen egal a 4A, mais le courant aspire par le pare de 
batterie egale a 18. 5A. Ce dernier sert a alimente les deux regulateurs, l’onduleur et la 
chambre. Dans le cas inverse ou la chambre est a F arret, les batteries delivrent un courant 
faible de 0.7 A pour alimenter les LEDs et le tableau de reglage de temperature dans la 
chambre. 
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Temps de fonctionnement, 26/04/2011 


Figure 4.1.31: Variation nocturne du courant en fonction de temps de fonctionnement de la 


chambre solaire. 
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Figure 4.1.32: Variation diurnal du courant en fonction de temps de fonctionnement de la 

chambre solaire. 

Quand il y de 1’ ensoleillement le courant delivre par les deux champs est suffisamment 
pour alimenter la chambre et charger le reste dans les batteries. 
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d) Caracterisation du temps de fonctionnement de la chambre 

Les thermocouples ont ete places comme suit : entree-sortie compresseur, entree-sortie 
condenseur, entree-sortie evaporateur, entree detendeur, chambre inferieure et l'ambiant. 
L'acquisition de donnees a ete faite en utilisant Hydra 2700Datalogger avec un enregistrement 
toutes les 30 secondes. 

Le systeme frigorifique met 27minutes pour baisser la temperature de la chambre a 4°C, 
en revanche, il met 15 heures pour baisser celle de l’eau jusqu’a 7°C. Durant cette phase le 
systeme fonctionne reellement pendant cinq heures sans arret et consomme une puissance 
active de 2kWh. Le but est de maintenir la temperature de la chambre dans une gamme de 
4°C a 8°C pour une periode la plus grande possible durant laquelle le systeme est a l'arret. 

Dans un premier temps, on a remarque que le temps d’ arret du systeme est de 5 minutes 
par cycle quand la temperature de l’eau est superieure a 7°C, par contre il est environ de 1 
heure quand cette demiere est environ de 7°C. Cela peut s'expliquer par la capacite calorifique 
de l’eau (140kg d’eau) qui cede sa chaleur dans la chambre entrainant une augmentation 
rapide de temperature 1 °C/min. 
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Figure 4.1.33: Test de systeme avec une charge de 200kg d’eau 
Afin de determiner le temps necessaire au systeme et de baisser la temperature de la 
chambre et de l’eau, on a met le systeme a l’arret et recharger le refrigerateur avec 200kg 
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d’eau a 25°C, et la meme temperature pour l’ambiant et de la chambre. Le role de systeme est 
ramene ces temperature a cette gamme (4°C - 8°C), pour cela il a met environ 6 heures pour 
baisser celle de la chambre, par contre, il a met 34 heures pour baisser celle de l’eau. On peut 
dire que le temps de refroidissement de l’eau est six fois plus grand que celui de Fair de la 
chambre, ce qui traduit en termes de consommation. 


e) Cout total de l’installation solaire 

Le cout global de la chambre frigorifique solaire photovoltaique realise a ete calcule en 
basant sur le prix unitaire d’ achat de chaque equipement. Le prix total de telle installation est 
de 2744859.63 DA, cette technologie reste chere en Algerie, raison pour laquelle elle est plus 
exploitable par les militaires que les services civils. 


Tableau 4.1.9: Devis estimatif de la chambre solaire photovoltaique realise. 


Equipements 

Nombre 

Prix unitaire 

(DA) 

Total en (€) 

Chambre 

01 

242787.82 

242787.82 

Panneaux photovoltaiques 

14 

71147.37 

996063.16 

Batteries solaires utilisees 

24 

34446.25 

826738.28 

Onduleur sinusoidal 

01 

100818.30 

100818.30 

Regulateur 

02 

23844.07 

47688.14 

Disjoncteur 

04 

17529.38 

70117.50 

Presses etoupes panneaux et 
boite de raccordement 

04 

325.22 

4553.08 

Cosses 

08 

42.28 

433.63 

prise europeenne 

03 

1951.32 

5853.97 

Prix total TTC 

{TV A 19.6%) 2580000 DA 


4.I.4.4. Test sous les conditions climatiques d’Ain-Guezzam (Tamanrasset) 

Apres le bon fonctionnement du systeme frigorifique solaire photovoltaique sous les 
conditions mediterraneennes, le meme systeme sans redimensionnement a ete met en test sous 
des conditions climatiques arides dans une caserne de GGF a Ain-Guezam. Apres deux mois 
de test, quelques problemes ont ete rencontres, tel que l’endommagement de l’onduleur du 
principalement a la hausse de la temperature enregistree qui a influence sur la consommation 
totale du systeme. Mais en generale, le systeme a bien fonctionne pour le refroidissement de 
l’eau pour plus de 100 gendarmes en periode caniculaire ou la temperature ambiante a 
depasse les 50°C a l’ombre. La figure 4.1.32 represente le systeme frigorifique solaire au 
cours de test sans collecte de donnees. 
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Figure 4.1.34: Unite de test experimental a Ain-Guezzam 


4.1.5. Conclusions 

Une chambre frigorifique de 1.5 m a ete realisee et dimensionne pour etre alimente en 
energie solaire photovoltaique, ce systeme a ete teste dans deux regions differentes en 
Algerie. Les resultats suivants ont ete obtenus a partir de l'etude experimentale: 

• Au cours de la journee, le pic le plus eleve de l'energie produite par les panneaux 
photovoltaiques correspond au pic de consommation de la chambre. 

• L'energie necessaire pour le refrigerateur est fournie par le solaire, done il n'y a pas 
d'impact negatif de ce systeme frigorifique sur l'environnement. 

• Les tests effectues sous les conditions climatiques de deux regions de test, nous 
montrent la performance de la chambre frigorifique realisee. 
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Contribution 2. : Prevision neuronale directe et inverse du coefficient de 
performance solaire optimal d’un systeme frigorifique a adsorption. 

4.2.1. Dimensionnement d’une machine frigorifique a adsorption 

Actuellement, la technologie nous offre une variete de moyens permettant d’obtenir le 
froid. Ainsi, de nombreux scientifiques et ecologistes se prononcent done pour une 
technologie capable d’ assurer le futur ecologique de notre planete. La disponibilite et la 
gratuite de l’energie solaire, d’une part, et la protection de l’environnement, d’ autre part 
pourraient rendre les machines a adsorption de plus en plus competitives dans le marche 
relatif a cette industrie. C’est dans ce concept que se situe notre travail dont le but est de 
construire une machine frigorifique solaire (la realisation a ete arretee dues aux problemes de 
soudage rencontres) et les etapes de dimensionnement sont presentees dans 1’ annexe B. 
Raison pour laquelle la partie de construction a ete remplace par la partie modelisation. 

Cette demiere partie consiste a utiliser les reseaux de neurones artificiels (RNA) directs 
pour predire les valeurs du coefficient de performance solaire (COP s ) d'un systeme de 
refrigeration a adsorption fonctionnant par intermittence et utilisant le couple charbon actif- 
methanol comme couple destine a la conservation des medicaments. Cette modelisation 
neuronale a pour but d’etablir un modele mathematique capable d’estimer la valeur de COP s 
de n’importe quelle machine frigorifique sous n’importe quelles conditions operatoires ou 
geographiques. 

La strategic d'obtention du RNA s'articule sur l'elaboration d'un programme sur 
MatLab®, comportant plusieurs boucles ou on fait varier les algorithmes d’apprentissages, les 
fonctions d’ activations, le nombre de couches cachees et le nombre de neurones dans chaque 
couche, afin de minimiser la fonction du cout sous contrainte d'une erreur quadratique 
moyenne (REQM) et un coefficient de correlation (R 2 ) fixes prealablement. Le modele 
neuronal confu a permis d'une part de reproduire avec une tres bonne precision les donnees 
experimentales tirees directement de litterature et d'autre part une estimation meilleure et plus 
precise des valeurs calculees par rapport aux modeles classiques de la thermodynamique. 

Apres la construction du modele neuronale, e’est-a-dire determiner leurs poids 
synaptiques et leurs biais, nous envisageons a developper une interface graphique conviviale 
en langage MatLab® pour rendre cette methodologie facile aux utilisateurs. 

En fin, la methodologie du reseau de neurones inverse (RNAI) sera adoptee pour 
trouver la tracabilite du COP s . 
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4.2.2. Les parametres devaluation de la machine frigorifique a adsorption 

Trois parametres principaux ont ete utilises par les chercheurs afin d'evaluer la 
performance d’une machine frigorifique solaire a adsorption (Wang et al. 2012). Ceux-ci 
etaient les suivants: l'energie solaire re£ue par le capteur de l’adsorbeur, la capacite de 
refroidissement ( SCP ) et le coefficient de performance solaire ( COP s ). 

L'energie re£ue a partir du rayonnement solaire est calculee comme la somme des 
contributions du produit de l'irradiation, le temps et la surface de captage, equation (4.2.1) 
(Baiju. 2013). 

Q k =^G,x(x4 (4.2.1) 

i = 1 

La capacite de refroidissement de l’evaporateur est donnee par l’equation (4.2.2) (Baiju. 
2013) : 

Qev ~ m methnol (^/ C p At) (4.2.2) 

Le coefficient de la performance solaire est defini comme la capacite de refroidissement 
de l'evaporateur a l'energie re£ue a partir du rayonnement solaire, equation (4.2.3) (Baiju. 
2013) : 

COP =9ev_ (4.2.3) 

Qr 

Alors que la puissance de refroidissement specifique par kilogramme d’adsorbant (SCP) 
est donnee par l’equation (4.2.4) (Baiju. 2013): 

SCP(W/kg) = 1000— (4.2.4) 

m AC 

Ou Mac est la masse du charbon actif. Vu le manque des valeurs de SCP dans les 
articles utilises, on a pris que le COP s comme parametre de control du systeme. 

4.2.3. Modelisation par les reseaux de neurones artificiels 
4.2.3. 1. Collecte des donnees 

L’objectif de cette etape est de recueillir des donnees pour developper et tester le RNA. 
La base de donnees experimentale (de mille valeurs (1000)) concemant les huit parametres 
experimentaux qui caracterisent la machine frigorifique a adsorption a ete obtenue a partir des 
articles scientifiques suivants, (Luo et al. 2005; Lemmini et al. 2007; Lemmini et al. 2005; 
Hassan et al. 2011; Dai et al. 2003; Antonio et al. 2000; Leite et al. 2004; Sumathy et al. 1999; 
Yong et al. 2004; Dai et al. 2003b; Sumathy et al. 2003;Vasta et al. 2008; Boubakri. 2003; 
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Boubakri et al. 1992a, b; Meunier. 1994; Boubakri. 2006; Anyanwu et al. 2003; Li et al. 
2004a, b; Li et al. 2005; Li et al. 2002; Li et al. 2003; Wang et al. 2002; Li et al. 2006). 

La base de donnee est realisee par une matrice rectangulaire de dimension (1000 ; 9), 
repartie de la fa§on suivante: 

• (1000 ; 8 ) pour la matrice d’ entree; 

• (1000 ; 1 ) pour la matrice de sortie ; 

Le choix des variables d'entree est motive par leurs valeurs mesurables, existantes et 
communes dans les articles selectionnes, tandis que la sortie est l'un des parametres 
difficilement mesurable, dont la valeur commune est donnee dans la litterature et qu’on peut 
utiliser pour evaluer la performance des machines a adsorption solaire. 

4.2.3.2. Analyse des donnees 

II est souvent preferable d’effectuer une analyse des donnees de maniere a determiner 
les caracteristiques discriminantes pour detecter ou differencier ces donnees. Ces 
caracteristiques constituent 1’ entree du reseau de neurones. Cette determination des 
caracteristiques a des consequences a la fois sur la taille du reseau (et done le temps de 
simulation), sur les performances du systeme (pouvoir de separation, ...), et sur le temps de 
developpement (temps d’apprentissage). 

Une etude statistique sur les donnees peut permettre d’ecarter celles qui sont aberrantes 
et redondantes. 

4.2.3.3. Preparation et partage de la base de donnees 

Afin de developper une application a base de reseaux de neurones, il est necessaire de 
disposer de deux bases de donnees : une base pour effectuer l’apprentissage et une autre pour 
tester le reseau obtenu et determiner ses performances. 

Afin de controler la phase d’apprentissage, il est souvent preferable de posseder une 
troisieme base de donnees appelee "base de validation croisee" (Pados et al. 1994). 

Afin d’arreter l’apprentissage juste avant que ne se produise ce phenomene de 
surapprentissage, plusieurs methodes ont ete proposees. La plus simple consiste a disposer de 
trois bases de donnees distinctes : une base d’apprentissage, une base de test et une base dite 
de "validation croisee". Cette derniere base est utilisee pendant l’apprentissage afin 
d’ examiner le comportement du reseau pour des donnees qui lui sont inconnues. Ainsi, 
l’apprentissage est arrete lorsque l’erreur atteint un minimum. 
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Notons que cette technique necessite d’ avoir suffisamment de donnees pour constituer 
trois bases a la fois representatives et distinctes. Dans le cas ou on ne disposerait pas 
suffisamment de donnees, une technique a ete developpee pour n’utiliser que la base 
d’apprentissage; elle consiste a effectuer deux apprentissages successifs. Pour le premier 
apprentissage, la base d’apprentissage est divisee en deux de maniere a constituer une plus 
petite base d’apprentissage et une base de validation croisee. Le nombre de cycles 
d’ apprentissage necessaire pour parvenir a la zone de generalisation optimale est conserve. Ce 
nombre est ensuite utilise pour stopper 1’ apprentissage avec la base d’apprentissage globale 
(Lang et al. 1990). 

II n’ y a pas de regie pour determiner ce partage de maniere quantitatif. II resulte souvent 
d’un compromis tenant compte du nombre de donnees dont on dispose et du temps imparti 
pour effectuer 1’ apprentissage. 

La figure 4.2.1. Illustre le pourcentage de partage de la phase d'apprentissage, test et 
validation. 
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Figure 4.2.1: Pourcentage de partage de la base de donnees 

4.2.3.4. Apprentissage et optimisation du reseau de neurones 

L' apprentissage supervise est retenu dans cette etude, il consiste a foumir au reseau de 
neurones des jeux de donnees d'entrees pour lesquels on connait les reponses souhaitees. 

On construit alors un critere d'erreur en comparant la reponse du RNA a celle de 
1’ experience. Un algorithme de correction des poids synaptiques va permettre de diminuer 
cette erreur et d'ameliorer ainsi la reponse du RNA. La mise en oeuvre de cette technique 
impose la constitution d'une base d'apprentissage formee de l'ensemble de jeux de donnees et 
des sorties associees. 

L’ architecture du RNA est optimisee par une methode de validation croisee appliquee a 
la base d'apprentissage/validation. Le nombre optimal de neurones caches est defini comme le 
nombre pour lequel l'erreur quadratique moyenne est minimale (essai-erreur (Lee et al. 1999)). 
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Sur la base de la theorie presentee au chapitre des RNA, un organigramme a ete developpe 
sous MatLab® afin d’optimiser les parametres de reseau de neurones artificiel (figure 4.2.2). 

Pendant la premiere etape de 1’ optimisation, nous avons teste les parametres 
fondamentaux afin de savoir lesquels donnent la meilleure performance. Parmi ces parametres, 
on cite la fonction d’ activation pour chaque couche et les algorithmes d’apprentissage. Les 
fonctions d’ activations qui ont ete employees pour ce type de reseau sont la fonction TANSIG 
pour la couche cachee et la fonction PURELIN pour la couche de sortie. Tous les deux sont 
des choix tres communs pour ce type de RNA-PMC, (Vogl et al. 1988; Hagan et al. 1996; 
MathWorks. 2004). 

Le choix d’algorithme d’apprentissage et une autre tache plus compliquee que le choix 
des fonctions d’ activation. En outre, il y a un grand nombre des algorithmes d’apprentissages 
qui peuvent produire un reseau precis, rapide et fiable. Afin de decouvrir quel algorithme 
execute mieux, dans notre cas on a teste six algorithmes. On a donne des noms aux reseaux 
crees tels que, netl, net2, net3, net4, net5, net6, les algorithmes d’apprentissages avec leurs 
correspondants codes sont presents dans le tableau 4.2.1. Plusieurs chercheurs mentionnent 
qu’il n’y pas de regie pour le choix de nombre de neurones dans la couche cachee mais on le 
determine empiriquement apres une phase d’apprentissage. Dans notre cas, il a ete estime 
avec douze neurones en appliquant la formule de (Kalogirou et al. 1996a, b ; Kalogirou et al. 
1997 ; Kalogirou et al. 1998 ; Kalogirou et al. 2000), Il s'agit d'une formule empirique, qui a 


ete utilisee dans un certain nombre de problemes d'ingenierie pour la modelisation et la 


prediction et a donne de tres bons resultats. Ce nombre sera question d'essais et d'erreurs dans 


les etapes qui suit. 


Tableau 4.2.1: Noms des codes et des algorithmes d’apprentissages correspondants 


Nom 

Nom complet de l’algorithme d’apprentissage 

Code MatLab® 

Acronyme 

Netl 

Levenberg-Marguardt 

Trainlm 

LM 

Net2 

Scaled conjugate Gradient 

Trainscg 

SCG 

Net3 

regularisation bayesienne 

Trainbr 

BR 

Net4 

Resilient Backpropagation 

Trainrp 

RP 

Net5 

BFGS Quassi-Newton 

Trainbfg 

BFG 

Net6 

Variable Learning Rate Backpropagation 

Traingdx 

GDX 


Bien que les RNA font partie des modeles de boites noires qui represented une relation 
empirique entre des ensembles d’ entrees et de sorties, le choix de variables d’ entrees et de 
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sorties est decisif pour que le modele neuronal se rapproche des modeles deterministes. Par 
consequent la selection des donnees experimentales des variables d’ entrees et de sorties et le 
pretraitement des donnees constituent la premiere etape de la modelisation neuronale. L’ etape 
suivante concerne le choix des parametres du reseau de neurones qui vont permettre au reseau 
d’apprendre a partir des donnees selectionnees la relation entre les entrees et les sorties 
choisies. La troisieme etape est consacree a 1’ operation d’ apprentissage et le test de 
generalisation du reseau choisi. Enfin la performance du reseau peut etre, eventuellement, 
comparee avec celles des autres types de modelisation. La figure 4.2.3 resume les etapes de la 
methodologie utilisee dans ce travail pour la modelisation neuronale de la performance solaire 
des machines frigorifiques solaires a adsorption. 



Figure 4.2.3: Procedure totale d’elaboration d’un RNAS. 
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Figure 4.2.2: Phase d’ optimisation des parametres du RNA. 
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4.2.3.5. Normalisation 

De maniere generate, les bases de donnees doivent subir un pretraitement afin d’etre 
adaptees aux entrees et sorties du reseau de neurones. Un pretraitement consiste a effectuer 
une normalisation appropriee, qui tienne compte de 1’ amplitude des valeurs acceptees par le 
reseau. Les donnees des variables selectionnees ont ete normalisees entre Y min =-1 et Y max =l 
par F equation de normalisation suivante : 



norm 


X — X 

old min f y _ y A.y 

v y v 1 max min / min 

max min 


(4.2.5) 


La fonction Premnmx deMatLab® realise la normalisation selon F equation 4.2.5, la 
fonction inverse de denormalisation est realisee par la fonction Postmnmx. Comme les 
reseaux de neurones ne peuvent pas etre extrapoles, nous avons represente dans le tableau 
4.2.2 les valeurs minimales et maximales des entrees utilisees. 


Tableau 4.2.2: Minimums et maximums des entrees utilisees 


Intervalle de validite des entrees 

Intervalle de validite de la sortie 

Parametres 

[Min-Max] 

Parametre 

[Min-Max] 

A c (m 2 ) 

0.50 a 1.50 

COPs 

0.04 a 0.24 

T e (°C) 

70 a 120 


m A c (kg) 

14.50 a 38 

I (MJ/m 2 ) 

11.60 a 29 

T c (°C) 

20 a 5 

Latitude 

-34.95 a +121.42 

Longitude 

7.70 a 29 

T a (°C) 

7.14 a 36 


4.2.3.6. Choix des parametres du reseau et d’apprentissage 

La base des donnees qui est constituee de 1000 points experimentaux, est divisee sur 
trois sous-ensembles bien distribues sur l’espace de donnees {670 points pour l’apprentissage, 
200 pour le test et 130 pour la validation}. F architecture du RNSO est donnee dans le tableau 
4.2.3. 


Tableau 4.2.3: Topologie du RNA optimise pour la prediction du COPs. 


Type du reseau 

RNFF-BP (fonction newjfde MatLab®) 

Couche 

Nb. de neurones 

Fonction d’ activation 

Couche Entree 

08 

— 

Couche Cachee 

08 

Sigmoide (fonction Tansig de MatLab®) 

Couche Sortie 

01 

Lineaire (fonction Purelin de MatLab®) 

Algorithme 

d’apprentissage 

BRBP utilisant F optimisation de Levenberg-Marquardt. 

(fonction trainlm de MatLab®). 
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Figure 4.2.4: Reseau de neurones multicouches pour le calcul du COP solaire. 



COP . 


de S neurones 


S neurones. 


Figure 4.2.5: MLP pour la modelisation de COP. 


La figure 4.2.2 represente les etapes d’ optimisation de differents parametres du RNSO, 
ces etapes sont presentees sous forme d’ organigramme ou elles sont traduites par un 
programme sous MatLab®. La figure 4.2.3 resume les differents parametres optimises d’un 
perception multicouche avec huit neurones d’ entrees, huit neurones caches et un neurone de 
sortie. 

Apres la phase de l'apprentissage, le reseau de neurones artificiel optimise (RNSO) est 
teste sur de nouveau jeux d'entree/sortie. Cette etape permet de verifier le comportement du 
RNA sur des cas non appris. 

Pour analyser plus precisement les performances du reseau de neurones statique 
optimise du RNSO, l'erreur relative absolue en pourcentage {REQM (%)} a ete calculee. 

f M 


Out , =COP s =^ ] 
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V 


k= 1 


J 


- 105 - 


Chapitre 4. Partie experimental - Contribution 2 : Prevision neuronale directe et inverse du coefficient de 
performance solaire optimal d’un systeme frigorifique a adsorption 


Ou s est le nombre de neurones dans la couche cachee (5 = 5), k est le nombre de 
neurones dans la couche d'entree (K = 8), l est le nombre de neurones dans la couche de sortie 
(1 = 1), Wi, W 0 et bl( S ), b2(i) sont des poids et des biais de la couche cachee et la couche de 
sortie respectivement. Dans ce travail, un RNA multicouche avec une couche cachee a ete 
utilise pour tous les ensembles de donnees. La performance du modele neuronale a ete 
statistiquement mesuree par l'erreur quadratique moyenne ( REQM) et le coefficient de 
determination (. R 2 ), ces deux sont calculees avec les valeurs experimentales et les predictions 
du RNSO. Ces calculs sont utilises comme critere d’ evaluation du modele. 



Ou Q est le nombre de points de donnees. 

Le RNSO est obtenu quand YREQM est minimale et le R 2 est maximale. La matrice des 
poids synaptiques et les vecteurs des biais, des trois couches cachees et de la couche de sortie, 
sont regroupes dans le tableau 4.2.4. 


Tableau 4.2.4: Matrices des poids synaptiques et vecteurs des biais du RNA optimise 






K = 8, S = 8 and / 

= 1 




Wi(s, k) 

5.4332, U) 

-2.3128,1,2) 

6.6 174,1,3) 

11,3074,1,4, 

9.9890(1,5) 

-7.0328(1,6) 

7.881 l(i,7) 

-13.3758(1,8) 


7.0175( 2 ,d 

2.6250(2,2) 

-9.2014,2,3) 

7.8745(2,4, 

1.5188(2,5) 

10.0217(2,6) 

3.7381(2,7) 

7.7021(2,8) 


10.2978,3,1) 

7.8674(3,2) 

0.1591(3,3) 

3.7913(3,4, 

-3.6712(3,5) 

9.4267(3,6) 

0.9406(3,7) 

-1 1.3688(3,8) 


3.9525,4,1) 

6.9891(4,2) 

9.9965(4,3) 

17.8020,4.4) 

3.8398(4,5) 

-6.4535(4,6) 

3.2026(4,7) 

2.8565(4,8) 


-6.9145(5,1) 

-0.5965(5,2) 

-3.6103,5,3) 

-3.6187(5,4) 

-8.6257(5,5) 

-1.4914(5,6) 

-3.7563(5,7) 

-19.5260(5,8) 


-0.1408(6,1) 

3.7399(6,2) 

10.0862(6,3) 

-13.4490(6,4) 

-3.9238(6,5) 

-4.2963(6,6) 

-4.1655(6,7) 

13.8741(6,8) 


-6.0809(7,d 

-8.6808,7,2) 

-6.7334,7,3) 

-8.61 13(7,4) 

2.5434(7,5, 

-7,1016(7,6) 

7,7391(7,7) 

-2.6564(7,8) 


2.6438,8,1) 

6.6025(8,2) 

-5.4891,8,3) 

0.9468(8,4) 

-6.7670(8,5) 

-7.3190,8,6) 

-9.1026(8,7) 

-9.2443(8,8) 

Wo(l, s) 

7.4348,1,1) 

-2.5090,1,2) 

-6.9018,1,3) 

-6.0313(1,4) 

-11,6041,1,5) 

-3.5476(1,6) 

7.7361(1,7) 

0.9714(1,8) 


bl(s) -10.1220(i) 


-16.5116 ( 2) 
-10.0202 ( 3) 
-21.0441(4) 
23.8909 C 5) 
0.6878(6) 
11.7904(7) 
16.6413( 8 ) 
b 2 ( i ) 2.4831,!) 
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4.2.4. Validation du modele neuronal 

La figure 4.2.6 montre une comparaison entre le COP solaire experimental tire 
directement de la litterature et celui predit par le modele de RNSO. Cette representation a ete 

faite par la fonction Postreg de MatLab® qui montre un accord entre les donnees estimees et 

2 

celles experimentales avec un vecteur proche de l'ideal, [a (pente), p (ordonnee a l'origine), R 
(coefficient de determination)] = [1.00000 ; 0.00014 ; 0.9985]. 

La courbe de regression lineaire de la sortie du RNSO avec la cible est visualisee dans 
la figure 4.2.6. La juxtaposition de la premiere bissectrice et de la droite du meilleur 
ajustement lineaire de la sortie avec la cible ainsi que la distribution de points experimentaux 
montrent P excellent accord entre le COP s calcule par le modele neuronal et le COP s 
experimental avec la REQM (%) egale a 0.0453. 



Figure 4.2.6: Comparaison entre les valeurs simulees par RNSO et celles experimentales 

tirees de la litterature. 


4.2.4. 1. Comparaison des resultats 

La figure 4.2.7 represente le COP solaire en fonction de 1’ irradiation, elle montre bien 
que les points calcules par le reseau de neurone sont juxtaposes sur celles de l'experimentales, 
ce qui confirme une autre fois la performance du reseau de neurone optimise. Une autre 
representation du COP solaire en fonction de la masse du charbon actif, de la temperature de 
condensation et de la longitude, est donnee sur les figures 4.2.7b, figure 4.2.7c et figure 4.2.7d 
respectivement. 
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Et afin de visualiser encore plus les capacites predictives de modele neuronal en phase 
de validation, nous avons etabli une deuxieme comparaison en termes d’ analyse statistiques a 
travers le calcul de la racine de l’erreur quadratique moyenne ( REQM) et l’erreur standard de 
prediction (ESP) en parallele avec le facteur de biais (Bf) et le facteur de precision (Af) 
(Kumar et al. 2001; Lu et al. 2011; Zafar et al. 2012a, 2012b). 

Avec ESP est l’Erreur Standard de Prediction: 


n 


Yy. 

REOM Z-f I.exp 

ESP (%)= — — xlOO Avec Y = 1=1 


Y 


n 


(4.2.9) 


Ou Y iyexp est la valeur experimentale; Y hpre d est la valeur predite; Y e est la moyenne des 
donnees experimentales ; n est le nombre de points experimentaux. 


B f = 10 


losiYi pred /Yj ex p )/ H 


(4.2.10) 


A, =10 


HJ iXi,pred / ,exp )| 


/ n 


(4.2.11) 


Avec ERAM est l’Erreur Relative Absolue Moyenne : 


n 


ERAM(%) = —Y 

n ,=i 


v 


exp 


3 ^ pred 




exp 


(4.2.12) 
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Figure 4.2.7a: Validation du RNSO tors de la simulation du COP solaire en fonction de 


l’irradiation journaliere. 
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Figure 4.2.8: Validation du RNSO en termes de differentes erreurs locales lors de la 

simulation du COP solaire. 
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Figure 4.2.9: Validation du RNSO en termes de differentes erreurs globales lors de la 

simulation du COP solaire. 


Les resultats de la comparaison entre le reseau de neurone statique optimise (RNSO) et 
les resultats experimentaux, en termes de la distribution des differentes erreurs locales, lors de 
la simulation du COP solaire sont presentes sur l’histogramme de la figure 4.2.8. Alors que la 
figure 4.2.9 represente la distribution des differentes erreurs globales. 
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Ces resultats indiquent que les deviations des predictions obtenues par le modele 
neuronal optimise sont nettement plus faibles, ce qui confirme la robustesse du modele etabli. 
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Figure 4.2.7b: Validation du RNSO lors de la simulation du COP solaire en fonction 

de la masse de charbon actif. 
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Figure 4.2.7c: Validation du RNSO lors de la simulation du COP solaire en fonction 

de la temperature de condensation. 
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Figure 4.2.7d: Validation du RNSO lors de la simulation du COP solaire en fonction 

de la longitude. 


4.2.4.2. Programme de calcul de COPs 

Apres avoir construit le modele neuronal, c’est-a-dire determiner ses poids synaptiques 
et ses biais, nous avons developpe un programme informatique en langage MatLab®. Cette 
demarche permet a l’utilisateur d’ avoir tous les inputs necessaires a 1’ execution de ce modele 
pour calculer la sortie. 

Cette demarche a ete faite dans le but d’apporter plus de souplesse a l’utilisation du 
modele neuronal elabore. Le programme va done permettre de faire un calcul simple et rapide 
du COP solaire. La methodologie suivie dans cette partie peut se resumer sur 1’ organigramme 
de la figure 4.2.10. 



Figure 4.2.10 : Etapes a suivre pour 1’ elaboration de 1’ interface graphique de RNSO par 

MatLab®. 
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Figure 4.2.11: Interface graphique developpee pour le calcul du COP solaire. 
4.2.4.3. Analyse de la sensibilite des entrees par rapport a la sortie 

Le RNSO a fourni une matrice de poids de connexion, les poids sont des coefficients 
relatifs aux neurones artificiels, sont analogues aux forces synapses entre les axones et les 
dendrites dans de vrais neurones biologiques. La matrice des poids issue du RNSO peut etre 
utilisee pour analyser la sensibilite des entrees par rapport a la sortie desiree. Garson (Garson. 
1991) a propose une equation basee sur le partitionnement de poids de connexion: 


(r \ \ 




(4.2.12) 


Ou Ij est sensibilite relative de la j-ieme variable d'entree sur la variable de sortie; N t et 
N/ r represented le nombre de neurones d'entree et cachees, respectivement; W s sont des poids 
de connexion, l'exposant i, h et o referent a la couche d'entree, cachee et de sortie 
respectivement, et les indices k, m, n et referent a la couche d'entree, cachee et de sortie 
respectivement. 

La sensibilite relative des differentes variables calculees par l'equation 12 est 
representee sur la figure 4.2.12. Comme on peut le constater, toutes les variables ont des effets 
importants sur le coefficient de performance. Par consequent, il peut etre observe que la 
contribution la plus elevee a ete obtenue avec I (t) (24.81%). La temperature de condensation 
et la masse de charbon actif ont la meme influence sur les valeurs COP s (22.11%), la 
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longitude (14.56%) et autres parametres qui ont moins d’influence sur la valeur de COP 
solaire, tel que la latitude. Par consequent, des efforts significatifs peuvent viser a 
l'amelioration de ce type de systemes de refrigeration au moyen de reseaux de neurones 
artificiels inverses. 


Longitude 

15 % 


TA (°C) Ac{m 2 ) Tg,(°C) 


Latitude 
4 % 


Tc{°C) 

22 % 



(AC)(kg| 
2 % 


/day) 


Figure 4.2.12: Importance relative des variables d’entrees sur la prediction du COP s de la 

machine frigorifique a adsorption charbon actif /methanol. 

4.2.4.4. Elaboration du RNA apres analyse de sensibilite 

Apres 1’ etude de la sensibilite des entrees par rapport a la sortie, il est done necessaire 
de prendre la meme demarche de F architecture du RNA, le meme algorithme d’apprentissage, 
les memes fonctions d’ activations et en eliminant les variables d’entrees qui presentent les 
faibles importances lors de la simulation neuronale du COP solaire, par consequent 
P architecture du reseau de neurones artificiels est resumee dans le tableau 4.2.5. 


Tableau 4.2.5: Topologie du RNA optimise apres la phase de sensibilisation 


pour la prediction du COPs. 


Type du reseau 

RNFF-BP (fonction newjfde MatLab®) 

Couche 

Nb. de neurones 

Fonction d’ activation 

Couche Entree 

04 

— 

Couche Cachee 

07 

Sigmoide (fonction Tansig de MatLab®) 

Couche Sortie 

01 

Lineaire (fonction Purelin de MatLab®) 

Algorithme 

d’apprentissage 

BRBP utilisant F optimisation de Levenberg-Marquardt. 

(fonction trainlm de MatLab®). 


Apres la phase d’apprentissage ou on a varie le nombre de neurones dans la couche 
cachee et le nombre de couche cachee. Un reseau de neurone a ete obtenu apres 1000 
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iterations avec la racine sur l’erreur quadratique moyenne est mois de 1% et un coefficient de 
correlation egale ( R=1 ). 



Figure 4.2.13: Reseau de neurones multicouches elabore apres 1’ analyse de sensibilite. 


4.2.5. Performance optimale au moyen d’RNAI 

Selon le modele RNA (equation 6), il est possible de simuler les performances du 
systeme lorsque les parametres d'entree sont bien connus. Dans la figure 4.2.12, il a ete 
constate que le rtiAc est le parametre d'entree le plus pertinent pour etre optimise par une 
methode inverse. Par consequent, il est important de savoir dans ce processus, quelle valeur 
mAc est necessaire pour un COP s souhaite. Toutefois thac n'est pas le seul parametre qui peut 
etre controle et optimise, depuis T c , longitude ou le prototype a ete teste sont egalement des 
parametres pertinents. I(t) et TA ont ete elimines afin d'optimiser car il n'est pas possible de 
les controler. 

Par consequent, nous avons developpe une strategic optimale pour estimer le mac et T c a 
partir du reseau de neurones artificiels inverses (RNAI). 

Les reseaux neuronaux artificiels inverses ont ete utilises avec succes dans des 
applications differentes dans lesquelles un parametre de sortie desiree est selectionne, puis les 
parametres d'entree sont determines (Hernandez et al. 2009a, b, c; Cortes et al. 2009; 
Hamzaoui et al. 2011). Un reseau general est constitue par la fonction de transfert TANSIG et 
PURELIN. Ensuite et afin d'eviter toute ambigui'te, la sortie de RNAI est representes ci- 
dessous etape par etape. 


COPs = 



1 + e 


-2(M, s ,+T (Wi^^x /« (i) )) 


-1 


+ b 2 


a) 


(4.2.13) 


Etant donne que les reseaux neuronaux artificiels sont pauvres en extrapolant, cette 
equation est valable dans la plage des donnees d'entree presentees lors de l’apprentissage 
(Tableau 4.2.2). Cette equation peut etre exprimee en Eq. (14). Ensuite, nous avons: 


COPs=b2 m -YW KU)+ Y 


2 XWr 


0 (l,s) 


2(bl is)+ Y(W iM x 


1 + e 


(4.2.14) 
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A cette etape, nous avons obtenu la fonction qui doit etre minimisee a zero pour obtenir 
les parametres optimaux d’ entrees : 


Fun(In(x )) = COPs + b2(l )~L s W o(m 


+L 


2xW, 


0 (l,s) 


1 + e 


~2(Wi ( , t) X/n(x)+J^ k W«s.t) X In fk) +b \s)'> 


(4.2.15) 


Ou x est la valeur rtiAc ou Tc qui peut etre calcule. Le modele propose par le RNA a 
plus d'un neurone dans la couche cachee, done, il est necessaire d'utiliser une methode 
d'optimisation (Hernandez et al. 2009b, c). Dans cette recherche, nous avons utilise une 
methode d'optimisation appliquee pour inverser le reseau de neurones artificiel avec huit 
neurones dans la couche cachee. L'optimisation du parametre d'entree est effectuee en 
utilisant l'algorithme du simplexe de Nelder-Mead pour l'optimisation sans contrainte de 
fonctions non lineaires (Hamzaoui et al. 2011). La methode de Nelder-Mead tente de 
minimiser une fonction objective non lineaire a variables multiples utilisant les valeurs de 
fonction /uniquement, sans aucun derive. La methode de Nelder-Mead est une adhesion a des 
methodes de recherche directe qui n'utilise pas de gradient numerique ou analytique. La 
fonction "fminsearch ()" est integree dans le MatLab® et travaille avec l'algorithme de Nelder- 
Mead (Lagarias et al. 1998). 

Un test a ete realise avec des donnees differentes pour optimiser le in,\c dans des 
conditions differentes de demontrer la faisabilite de l'utilisation RNAI. Les parametres 
simules ont ete compares avec les donnees experimentales afin de verifier l'exactitude des 
RNAI avec une l’erreur relative absolue moyenne donnee par 1’ equation suivante: 


ERA = 


QeXP(LU) Q 


SIM - par-ANNi 


Q 


xlOO 


EXP 


(4.2.16) 


4.2.5.I. Estimation de la masse du charbon actif par le RNAI 

Nous pouvons calculer la masse optimale mac pour un COP solaire requis, en utilisant 
le logiciel MatLab® R2008b et l’equation 17. Toutes les valeurs simulees par le RNAI sont 
presentees dans le tableau 4.2.6. 


f(m, c ) = -A + 
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(4.2.17) 


2 W 2 W 

^ vy 0(1,7) , ^^0(1,8) 
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A COPs b2 (l) +W 0(11) +Wo( 1>2 ) + 1^0(1, 3 ) "*"^o(i, 4 ) 
X t = -2 x (W j(u) x A c + W i(U2) xTg + Wj (M) X / + W i(US) 
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x 2 

X 3 

x 4 

X 5 

X 6 

X 7 

X 8 


— 2 x(Wj (21) x A c +W / ( 2 2) xTg + W^ 2j4 ) x/ + Wj ( 2 5 5> xT c + W/(2,6) x Lat + W i( ^ 2>7 ) xLong + W^ 2>8 ) xT^j+M 
— 2 x(w / ( 3j i ) xA c +W)( 3j2 ) xTg “*"^(3,4) x/ + W^ 3j5 ) "*"^/(3,6) x La/ - + W^ 3>7 ) x Long + W i( ^ 3 ^ xT A )+bl 

~ 2 x(w^( 4> i) xA c +^-(4 5 2) xZg + W]-( 4j4 ) x/ + W^ 4>5 ) xT c + W^ 46) x Lat + W i( ^ 4 7 ^ xLong +^-( 458) xT A )+M 
— 2 x(w /(5 i) xA c +W]. (5j2 ) xTg + W^ 5>4) x/+W^. (5>5) xT c +W]. (5>6) xLat + W i( ^ 51) xLong +W^. (5j8) x 7 ^ )+ bl 
— 2 x(w / ( 6 i) xA c +W]. (6j2 ) xTg +^-(6,4) "*"^(6,5) "*"^(6,6) x Lat + W]. (6>7) xLong +W]. (6>8) xT a )+M 

— 2 x(w]-( 7> i) xA c +^ (7i 2) xTg + W^ 74) x/ + W^- (7 5) xT c + W^ 7;6) xLat + W i(71 ^ xLong + W^ 7;8 ) xT A )+M 
— 2 x(W j -( 8 j) xA c +W]-( 8j2 ) xTg + W*( 8>4 ) x/+W^. (8 5) xT c +W J - (8 6) xLat +W i( ^ 7) xLong +W J - (8 8) xT A )+M 


( 2 ) 

(3) 

(4) 

(5) 

(6) 
(7) 

( 8 ) 


Tableau 4.2.6: Validation du modele neuronal inverse. 


A c (normalise par 1’ equation 5) 

00.525 

00.4000 

00.4000 

00.5250 

T g (normalisee par 1’ equation 5) 

00.600 

00.4400 

00.4400 

00.6000 

Mac simulee par l’equation 17 

23.9800 

37.9696 

18.9848 

14.4884 

I (normalisee par 1’ equation 5) 

00.442 

00.4010 

00.5448 

00.5240 

T c (normalisee par l’equation5) 

00.600 

00.5067 

00.5067 

00.6000 

Latitude (normalisee par l’equation5) 

00.550 

00.4014 

00.4014 

00.5501 

Longitude (normalisee par L equation 5) 

00.576 

00.4882 

00.4882 

00.5760 

TA (normalisee par L equation 5) 

00.400 

00.4743 

00.5293 

00.5239 

La valeur du COP s desiree 

00.1140 

00.1400 

00.1200 

00.0600 


Les masses de charbon actif estimees par le reseau de neurones artificiel inverse (RNAI) 
et celles de 1’ experimentale sont presentees sur la figure 4.2.14 avec une erreur relative de 
0.08%. Les resultats de la comparaison entre le RNAI et les donnees experimentales, en 
termes de la REQM et 1 'ESP, sont presentes sur l’histogramme de la figure 4.2.15 et la figure 
4.2.16. Ces resultats indiquent que les deviations des predictions obtenues par le modele 
neuronal optimise sont nettement plus faibles, ce qui confirme la robustesse du modele 
inverse etabli. Ce qui donne une erreur trop faible lors du calcul de COP solaire pour les 
valeurs de la masse de charbon actif estimee par le RNAI. 


40 


T 


T 


O masse experimentale 
• masse simulee par le RNSI, (ER=0,08%) 


<§> 


35- 


CD 


o 


30- 


o 

_Q 

cB 25 

sz 

o 


<§> 


"O 

<u 

% 20 
03 


<D 


15- 




0,06 


0,08 


T 

0,10 

COP 


0,12 


0,14 


- 116 - 


Chapitre 4. Partie experimental - Contribution 2 : Prevision neuronale directe et inverse du coefficient de 
performance solaire optimal d’un systeme frigorifique a adsorption 


Figure 4.2.14 : Simulation de la tragabilite de quelques valeurs du COP solaire en 

determinant la mase de charbon actif correspondante. 


3 r 5 

.1 

2 

1,5 

1 

0.5 

0 



■ REQM(°/c) 

■ FSP{S'n) 

■ Af [-)/ Bf(-) 

■ ERAM[%) 


Figure 4.2.15: Validation du RNAI en termes de differentes erreurs 



Figure 4.2.16: Validation du RNAI en termes de differentes erreurs 


4.2.5.2. Estimation du temps de calcul du RNA direct et inverse 

D’apres le tableau 4.2.7, on constate que le temps d’ optimisation des parametres du 
RNAI ne represente qu’un 1/10 de temps d’ optimisation du RNA, ceci est du aux differents 
parametres de 1’ architecture du RNA a optimiser (figure 4.2.2). 

Le temps d’ optimisation est plus langue en comparaison avec celui de 1’ execution, et 
cela reste pour un programmeur trop familier avec les RNA. 
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Tableau 4.2.7: Recapitulation de temps de programmation et d’execution pour le RNA direct 

et inverse 


En raison de huit heures par jour 

RNA 

RNAI 

Rapport 

(tRNAl/tRNA ) 

Temps de programmation, ( t p ) 

1 jour 

1 jour 

- 

Temps d’ optimisation, (t op ) 

20 jours 

2 jours 

10% 

Temps de calcul de simulation, (t e ) 

120 secondes 

10 secondes 

8.3 % 

Rapport, (top/t e ) 

14400 

17280 

- 


Une comparaison entre le RNSO et celle de la litterature a ete effectuee pour confirmer 
la performance du modele propose. 

Le temps de calcul est plus important par rapport a celle de la litterature, ce qui est du 
principalement aux proprietes de l’ordinateur employe par Hernandez. 


Tableau 4.2.8: Comparaison entre les resultats de RNSO et ceux de la litterature 


References 

Hernandez et al. 2012 

Laidi et al. 2013 

Topologie 

6-6-1 

8-8-1 

REQM(%) 

- 

0.0453 

Test statistique de la pente et 
l'ordonnee (%) 

99 

100 

R 2 (RNSO) 

0.9875 

0.9985 

Temps de calcul de RNAI (seconde) 

0.2900 

10 


La modelisation neuronale a ete faite sur un micro-ordinateur de type HP Pavilion 
dv2000, les caracteristiques essentielles de ce micro-ordinateur sont resumees dans le tableau 
4.2.9. 

Tableau 4.2.9: Proprietes du micro-ordinateur utilise pour la modelisation neuronale. 


Carte graphique 

Mobile Intel (R) 945 Express Chipset Family 

Processeur 

Intel (R) Core (TM) 2 CPU T5500@ 1.66 GHz 

Disque Dur 

Espace total : 37.6 Go / Espace utilise : 34.3 Go 

Memoire vive 

Memoire totale : 2068956 ko 


4.2.6. Conclusion 

Dans cette contribution nous avons decrit en detail une procedure pour la modelisation 
neuronale du COP solaire de la machine frigorifique a adsorption en utilisant un 
programme MatLab® elabore dans le cadre de ce travail, les donnees utilisees pour 
l’apprentissage/ test et validation ont ete tirees directement de nombreuses publications 


- 118 - 









































Chapitre 4. Partie experimental - Contribution 2 : Prevision neuronale directe et inverse du coefficient de 
performance solaire optimal d’un systeme frigorifique a adsorption 


scientifiques. La methodologie adoptee pour la determination du modele neuronal optimal est 
basee sur une recherche par essai et erreur de la topologie du RNA ayant une capacite de 
generalisation acceptable pour une erreur de prediction fixee a priori. Nous avons egalement 
decrit la methode devaluation de la performance du RNA optimise sous forme de calcul de 
vecteur du meilleur ajustement lineaire entre les donnees experimentales de la solubilite et 
celles calculees par le RNA optimise, et sous forme de calcul d'erreur des valeur calculees par 
rapport aux valeurs experimentales par le reseau de neurones statique optimise ( RNSO). 

De ce fait la modelisation neuronale pour le COP solaire s’avere une alternative, fiable, 
facile. Les conclusions relatives a la strategic a adopter pour arriver au modele neuronal avec 
une capacite de generalisation acceptable et une performance (en termes d’erreur) fixee a 
priori, sont comme suit : 

• Cette etude confirme les resultats de travaux anterieurs indiquant que les regies 
heuristiques, pour le choix du nombre de couches cachees et du nombre de neurones par 
couche cachees, n’ont aucune credibility, du moins pour les reseaux de type feedforward. 
La regie la plus raisonnable est la regie d’essai et erreur en allant du reseau a structure 
simple (une seule couche cachee et un nombre faible de neurones) vers des structures a 
plusieurs couches cachees, en fixant comme objectif une bonne capacite de generalisation 
et une performance fixee. 

• Pour attenuer le probleme de convergence vers les minimums locaux, il est conseille de 
realiser l’apprentissage du RN ayant une structure donnee avec plusieurs initialisations 
differentes des matrices des poids et biais synaptiques. 

• Le RNSO dont la conception de la topologie est optimisee en terme: 

- 08 neurones dans la couche d'entree; 

- 08 neurones dans la couche cachee avec la meme fonction d'activation (Tansig); 

- 01 neurone de sortie avec la meme fonction d'activation (Purline); 

- 01 algorithme d'apprentissage ( Levenberg-Marquardt). 

- 01 nombre d'iteration de 1000; 

• Le RNSO ainsi obtenu presente une capacite d’ estimation interessante pour le COP 
solaire avec la racine sur l’erreur quadratique moyenne (REQM) egale a 0.045% et un 
coefficient de correlation R 2 > 0.9985. 

• Une etude de sensibilite des entrees par rapport a la sortie a ete faite, parmi les huit 
parametres utilises, quatre ont une faible importance sur la sortie, ce qui permet de 
construire un autre RNA avec quatre entrees pertinentes. Le RNA a ete obtenu avec 07 
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neurones caches apres 1000 iterations, la racine sur l’erreur quadratique moyenne moins 
de 1 % et un coefficient de correlation unitaire. 

• Le RNSO a ete inverse et optimise pour trouver la tra§abilite du COP solaire par la 
methode d’ optimisation de simplexe Nealder-Mead, cette phase a ete realisee par la 
fonction de MatLab® "fminseach" integree dans les soubroutine de MatLab®. 

• Un programme informatique en langage MatLab® a ete developpe a partir des poids 
synaptiques et des biais du RNSO elabore, permettant a l’utilisateur de calculer le COP 
solaire a partir de huit entrees a savoir : "surface du capteur/adsorbeur, temperature 
maximale au sein de l’adsorbeur, masse du charbon actif, irradiation journaliere re£ue sur 
un plan incline et sur la surface de captation du capteur/adsorbeur, temperature de 
condensation du methanol, latitude, longitude et la temperature ambiante 
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Conclusion generale 

Notre travail de these a pour but d’etudier et de modeliser un refrigerateur solaire, 
l’objectif de la premiere partie est de realiser un refrigerateur domestique qui fonctionne en 
energie solaire photovoltaique avec une faible consommation et un moindre cout. Pour cela, 
trois prototypes de refrigerations a compression de vapeur ont ete realises, le premier (DC- 
AC) c’est un refrigerateur de volume de 1601, alimente indirectement par F energie solaire en 
utilisant un onduleur solaire ; le deuxieme (DC-DC) est de meme volume et de meme 
alimentation avec modification dans le cycle frigorifique, sans onduleur et avec une capacite 
de stockage electrochimique minimale ; le dernier prototype (DC- AC), c’est une chambre 
froide de 1.5m , alimentee aussi indirectement par F energie solaire. Ces prototypes ont ete 
dimensionnes avec le logiciel PVSyst. V5.59 pour une autonomie de 3 jours sans 
ensoleillement et testes sous les conditions climatiques sahariennes et mediterraneennes. La 
consommation joumaliere du systeme frigorifique DC- AC est de 0.865kWh/jour et un cout 
global du systeme estime a 462779.46 DA. La consommation joumaliere du systeme 
frigorifique DC-DC sous les memes conditions est de 0.47 3kWh/j our et un cout global du 
systeme de 280436.35 DA, cela induit une reduction de 40%, c.a.d. une reduction de 
182343.1 IDA par rapport au premier systeme DC- AC. Cette reduction est due principalement 
au prix de F onduleur et a la reduction de la capacite de stockage electrochimique. D’apres les 
tests menes sur le conteneur concu (DC-AC), il s’est avere que ce dernier est capable de 
garder sa temperature interieure entre 2 et 8°C, avec une consommation moyenne de 
1.2kWh/jour et sous les conditions environnementales climatiques sahariennes (Ain- 
Guezzam : Tamanrasset). 

L’objectif de la deuxieme partie, est d’etablir un modele neuronale capable de predire le 
Coefficient de Performance Solaire (COP s ) d’une machine frigorifique a adsorption solaire. 
Cette modelisation s'articule sur l'elaboration d'un programme sur MatLab®, qui comporte 
trois etapes d’ optimisation: 

1. Optimisation des parametres fondamentaux : 

o les fonctions d’ activation pour chaque couche: le choix est porte sur la fonction 
Tansig pour la couche cachee et la fonction Purelin pour la couche de sortie ; vue que 
tous ces deux dernieres fonctions sont des choix tres communs pour ce type de reseau 
(RNMC), (Hagan et al. 1996; MathWorks. 2000). 

o les algorithmes d’apprentissage : l’algorithme Scaled Conjugate Gradient donne le 

resultat le plus performant parmi les six (06) algorithmes testes "Levenberg- 
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Marquardt; Scaled Conjugate Gradient ; Regularisation Bayesienne ; Resilient 
Backpropagation ; BFGS Quassi-Newton ; Variable Learning Rate Backpropagation". 

2. Optimisation de quelques fonctions de normalisation : la fonction normalisation 
MMM donne le meilleur resultat parmi les trois fonctions de normalisation testees : 
"MMM, MSDM et PCA". 

3. Optimisation d’autres parametres : 

o le nombre de neurones dans chaque couche cachee : huit (08) neurones dans la 
couche cachee donne le meilleur resultat parmi le nombre de neurones testes (de 1 
a 30). 

o le nombre de couches cachees : le perceptron monocouche donne une erreur 
minimale parmi les 05 couches cachees testees (0 a 4). 
o le nombre d’iteration (epochs): le test de 1000 iterations represente le meilleur 
resultat parmi le nombre iterations testees (500, 1000, 2000, 3000, 4000). 

Le reseau de neurones statique optimise (RNSO) est obtenu pour une topologie de (8, 8, 
1) apres 1000 iterations et vingt (20) jours d’ optimisation en raison de 08h/j; avec VREQM 
egale a 0.0453% et le R 2 >09985. 

Un programme en langage MatLab® a ete developpe a partir des poids synaptiques et 
des biais du RNSO elabore, permettant a l’utilisateur de calculer le COP solaire a partir de 
huit entrees a savoir : "Surface du generateur, temperature maximale dans le generateur, 
masse du charbon actif employee, irradiation solaire joumaliere, temperature de condensation, 
temperature ambiante journaliere moyenne, latitude et longitude de lieu". 

Afin de mieux controler la performance de cette machine frigorifique a adsorption et 
apres L etude de sensibilite des parametres d’ entrees sur la sortie, le reseau de neurones 
artificiels inverse a ete construit avec une erreur relative absolue moins de 0.08% et un temps 
de calcul de 10 secondes pour le cas du test. 

Comme perspectives a notre etude nous envisageons : 
o une etude statistique du rayonnement solaire en Algerie, afin d’etudier et analyser des 
machines frigorifiques solaires a n’importe quelle cordonnee geographique en Algerie. 
o une etude technico-economique et environnementale des systemes de refrigeration 
solaire photovoltaique et thermique en utilisant l’approche de L analyse de cycle de vie 
(ACV). 
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